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ABSTRACT 
A  new  method  for  design  of  control  systems  in  the  time  domain  is 
presented.   The  adjustable  parameters,  the  roots  and  the  zeros  of  the 
system  are  evaluated  simultaneously  either  by  hand  or  by  computer. 
The  method  is  applicable  to  any  system  configuration  and  is  able  to 
handle  a  large  number  of  adjustable  parameters.   The  method  is  analy- 
tical, any  degree  of  accuracy  can  be  obtained. 
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PART  I  ANALYSIS 
CHAPTER  I 
THEORETICAL  CONSIDERATION  OF  ALGEBRAIC  EQUATIONS 
1-1.  Introduction. 

The  response  of  feedback  control  systems  is  governed  by  the  complete 
set  of  the  closed  loop  pole- zero  configurations.   Since  the  zeros  are 
usually  known,  the  only  problem  is  the  poles.   For  linear  control  systems, 
the  problem  of  poles  is  conveyed  to  the  roots  of  the  characteristic  equa- 
tion which  is  an  algebraic  equation.   Several  methods  have  been  developed 

to  this  end  such  as  Guilemin's  method  ,  root-locus  method  and  Mitrovic's 

2  3  4 
method  '  s  ,   Guillemin8s  method  starts  the  synthesis  by  choosing  a  set 

of  pole-zero  configurations  from  the  consideration  of  specifications  and 
then  finds  out  the  network  needed  to  give  this  chosen  pole-zero  configura- 
tion.  The  root- locus  method  is  essentially  a  graphic  method  for  one  vari- 
able which  has  linear  relation  with  the  coefficients  of  the  characteristic 
equation.   Mitrovic's  method  is  also  a  graphic  method  but  handles  two 
variable  coefficients.   The  method  proposed  in  this  paper  is  rather  an 
analytic  method  and  is  able  to  handle  any  number  of  variables  which  are 
contained  in  the  coefficients  of  the  characteristic  equation.   The  analy- 
sis is  based  on  the  root-coefficient  relation  of  an  algebraic  equation,  but 
a  new  technique  of  partitioning  the  equation  is  employed. 

1-2.  Root-coefficient  Relationship. 

th 
For  a  n   order  algebraic  equation,  there  are  n  relations  among  the 

roots  and  the  coefficients  of  the  equation.   If  all  the  coefficients  are 

actually  known  numerically,  then  all  the  roots  are  fixed  in  such  a  way  that 


they  satisfy  the  root-coefficient  relations.   Conversely,  if  a  set  of 
roots  of  the  characteristic  equation  are  chosen,  then  all  the  coefficients 
of  the  equation  are  fixed  by  those  relations.   However,  if  variables  are 
contained  in  the  coefficients,  then  the  roots  of  the  equation  are  varied 
in  a  definite  relation  to  the  variation  of  the  variables.   For  example, 
if  one  coefficient  of  the  characteristic  equation  is  variable,  then  the 
relation  between  the  roots  and  the  variable  coefficient  can  be  plotted  by 
root-locus  method.   But  if  this  relationship  is  looked  at  from  a  different 
point  of  view,  an  analytical  relation  can  be  obtained.   Consider  a  third 
order  characteristic  equation  with  one  variable  coefficient  as  shown  in 
the  following  equation: 

where  $       and  Q      are  known  numerically,    Ji("'     is  a  function  of  7? 
which  can  be  varied.   Assume  - S.   s  ~S7  and  —  5j  are  the  three  roots  of 
the  equation,  then  the  root-coefficient  relations  are  as  follows: 

SA  +  s,$3  &sss,  =  ftd)  (1-2-2) 

Since  j-  (fij  is  adjustable,  then  equation  (1-2-2)  can  always  be  satisfied 
for  any  choice  of  the  roots  by  varying    j't(n)      accordingly.   However, 
equations  (1-2-1)  and  (1-2-3)  set  constraints  to  the  value  of  the  roots. 
Namely,  the  three  roots  must  satisfy  the  two  equations  (1-2-1)  and  (1-2-3). 


( 


If  the  three  roots  are  treated  as  the  variables  of  a  set  of  simultaneous 
equations,  then  there  are  three  unknowns  but  only  two  equations,  consequent- 
ly one  root  is  arbitrary.   This  arbitrary  root  can  be  chosen  independently 
if  there  is  no  restriction  to  the  value  of    "f("nJ  »  an<*  therefore  this 
root  can  be  treated  as  an  independent  variable,  in  the  root-coefficient 
relations.   In  the  root-locus  method,  the  independent  variable  is  j~  ( n ) 
while  in  this  point  of  view,  the  independent  variable  is  conveyed  to  one 
of  the  roots  of  the  characteristic  equation. 

The  above  reasoning  is  based  on  the  solution  of  simultaneous  alge- 
braic equations,  i.e.,   S    5.  and    3.?   are  treated  as  independent 
variables.   Since  those  variables  are  the  roots  of  an  algebraic  equation, 
a  complex  root  must  be  treated  as  two  variables  because  of  its  conjugate 
root.   If  one  variable  of  a  pair  of  complex  conjugated  roots  is  expected 
to  be  treated  as  the  independent  variable,  this  variable  cannot  be  chosen 
to  be  the  root  itself,  but  can  be  treated  as  follows.   Assume  —  S  and  -5, 
are  a  pair  of  complex  conjugate  roots  of  the  equation  (1-1),  those  two 
roots  can  be  transformed  to  other  two  variables  by  the  following  trans- 
formations : 

S(+  Sz  =  Z^Odn  (1-3) 

iAfc  =  ^n  (1-4) 

Where  \     and  C/Jn  are  known  as  the  damping  ratio  and  the  natural  frequency 
of  the  pair  of  complex  conjugate  roots  —5,  and  -_$,  .   By  this  trans- 
formation, both   y   and   Ujn  can  be  treated  as  independent  variables. 

In  equation  (1-1),  if    j'An)       was  also  a  fixed  number,  then  all 
roots  were  fixed  numbers  by  equations  (1-2)  and  no  root  was  arbitrary. 
The  arbitrary  roots  arise  from  the  introduction  of  variables  to  the  coeffi- 
cients.  If  one  variable  is  introduced  to  the  coefficients,  it  has  been 
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shown  one  toot  is  arbitrary.   For  the  same  reason,  if  two  variables  are 
introduced  to  the  coefficient  of  a  n   order  characteristic  equation, 
then  two  roots  are  arbitrary,,  and  for   JT)   variables  to  the  coefficients, 
jfY)       roots  are  arbitrary.   The  root- locus  method  is  the  case  of  one  vari- 
able while  Mitrovic's  method  is  the  case  of  two  variables,,  but  both  methods 
are  quite  different  ways  of  approach.   It  is  because  of  this  way  of  look- 
ing at  the  root-coefficient  relations^  an  analytic  method  can  be  derived 
and  is  capable  of  handling  any  number  of  variables  contained  in  the  co- 
efficients. 
1-3.  Partition  of  a  Characteristic  Equation. 

If  adjustable  parameters  are  introduced  to  the  coefficients  of  a 
characteristic  equation,  then  this  equation  can  be  partitioned  into  two 
equations.   One  of  them  is  formed  by  the  arbitrary  roots,  the  other  is 
formed  by  the  constrained  roots.   The  constrained  roots  here  are  defined 
as  those  roots  other  than  arbitrary. 

Consider  a  5th  order  characteristic  equation 

5S>  B^S^B^Vb^+^S+B^O   (1-3-1) 

Assume    —5,    —  S  >    —  5,   -  .Sand   _5r  are  t^e  roots.   The  negative 
sign  is  for  convenience.   By  expanding  the  equation 

Set  the  corresponding  coefficients  of  (1-5-1)  and  (1-5-2)  equal,  then  the 
root-coefficient  relations  are  obtained  as  follows; 

B3  =  (5I+5J(3J«^)+5I5J,+5JS^  ta+^+^'V)-^   (1"6"2) 


( 


B,  -  535^(5,+5z;+  ( s,fjy^  i  [<V5J(^+^)  +  V^5V[-£r    (i-6"3) 


B    =    >Xsj\  -f"   [U  +  *J5A   +Sf5z(5Jf^j]55  (1-6-4) 


BD—   5,5*5,5^  d-6-5) 


Assume  the  coefficients   (3,   and   [3  are  adjustable,  then  equation  (1-5-1) 
can  be  written  as  follows: 


Where   -f  (^   f^  \  and   "f  (  $   7?*)  are  functions  of  /;,  and  yi^    which 
are  adjustable  parameters.  g>  ,    g   and  Q      are  assumed  given  fixed 
values.   A  terminology,  "controlled  characteristic  equation"  is  given  to 
equation  (1-7)  on  the  reason  that  its  roots  are  controllable  by  the  ad- 
justable parameters.   Since  there  are  two  adjustable  parameters  intro- 
duced to  the  coefficients,  two  roots  can  be  chosen  arbitrarily.   Let  those 
two  arbitrarily  chosen  roots  be  (  -5,   )  and   (  -5  ) ,   then  the  other  three 
roots  (-5j),  (  --^n  )  and  (  ~S^-  )  are  not  arbitrary  but  constrained  to  the 
fixed  coefficients  (  B  >   B   anc*  £L  )  anc*  tne  arbitrarily  chosen  roots 
as  shown  in  the  last  section.   By  separating  the  arbitrary  roots  and  the 
constrained  roots  equations  (1-6-1),  (1-6-2)  and  (1-6-5)  can  be  arranged 
as  follows: 

From  (1-6-1) 

\+\  +  $S    —    BA-(  S,  +  5J  d-8) 

From  equation    (1-6-2) 
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From  equation  (1-6-5) 

It  can  be  recognized  that  (  $3+\  +  Ss      ),  (  SjS^+S^S^  4 S^Sj    )  and 
(  C>  SoS-   )  are  the  coefficients  of  a  third  order  algebraic  equation 
which  has  roots  of  -5,  ,  - jL  and  -3  .   Define 

r    -±-    s  s  5  d-ii-i) 


£,=§:  S3S^  +  S+5r  +  5^S3  a-"-*> 


C^  =  Sj+S^  +  S^  (i-n-3) 


Then  the  equation  of  the  constrained  roots  is  formed  by 

S?+   CS2+C,SK0   =  O  (M2) 


It  must  be  noticed  that  equations  (1-12)  has  the  roots  —  ix  ,  -5^  and  -5_ 

Substitute  the  definitions  of  (1-11)  into  (1-8),  (1-9)  and  (1-10)  and 
manipulating,  one  obtains 

C>   =  B.-(5,+5i)  (1-12-1) 


C,  =  B3  -Ci(  S.+Aj-SfS, 


(1-12-2) 


c0=  8,  A**  (1-12-3) 

Consider  equations  (1-12),  C0  >  C(  and  C^.  are  functions  of  the  arbitrary 
roots  and  the  fixed  coefficients   g  ,   g   and  &~    of  the  controlled 
characteristic  equation  (1-7).   Therefore  the  roots  of  equation  (1-12) 
are  also  functions  of  the  arbitrary  roots  -5,  and  -5,  and  the  fixed  co- 
efficients.  Assume  there  are  no  restrictions  to  the  variable  coefficients 
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"f  (f)      /;)and  f?(f)      n   )  then  the  controlled  characteristic  equation  (1-7) 
is  partitioned  into  the  following  equation 

because  both  equations  have  the  same  roots.   As  Q      ,  C,  and  Q     are 
functions  of  5  an<^    S2    ,  then  the  five  roots  of  the  original  character- 
istic equation  (1-7)  are  completely  determined  by  equation  (1-12).   Equa- 
tion (1-12)  is  a  third  order  algebraic  equation  and  the  coefficients  are 
expressed  analytically  by  (1-12-1),  (1-12-2)  and  (1-12-3),   Then  the 
original  5th  order  characteristic  equation  is  now  reduced  to  a  third  order 
equation  and  by  this  reason,  a  terminology  "REDUCED  CHARACTERISTIC  EQUATION',' 
is  defined  for  equation  (1-12). 

In  the  above  arguments  it  is  assumed  that  the  variables  ft    and  ft 
are  adjusted  such  that  two  arbitrary  roots  are  chosen.   Therefore,  the 
variables  p),     and  ft      ,  which  are  the  actual  physical  adjustable  parameters 
of  the  system,  must  be  determined.   From  equation  (1-6-3),  and  (1-6-4) 

replace  B^  by  fz  ( ^  ,  $z )   and   B,  by  ft(fitJ   $J  since  they  are 
assumed  adjustable,  and  manipulating,  one  obtains 

f-,(ti,M  =  C(5,  +  3J  +C,S,S,  d-14-2) 

Substitute  £     ,    C,  arid   Ca    from  equations    (l-ljj)    into    (1-14)  s    obtain 

t,W.XJ   =   ftU+y/sf5a  4  £3S,SZ  -  M^(*+V)*  5,5, (^^-(5^/(1-15-2) 

It   can  be   seen   from  equations    (1-15)    that       j    (  fh    "*)    anc*         JjC^/y^j.)  are 
expressed  analytically  as  a    function  of   the  arbitrary   roots   and   the   fixed 


c 
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coefficients  |30  ,   B3   and  84.  .   And  ^  ^  also  can  be  expressed  ex- 
plicitly as  functions  of  the  arbitrary  roots.   Therefore  if  a  set  of 
roots  has  been  chosen  from  the  consideration  of  the  reduced  characteristic 
equation,  then  the  physical  adjustable  parameters  are  readily  determined 
by  equations  (1-15)  or  (1-14). 

By  this  partitioning  process  it  can  be  seen  easily  that  the  original 
fice  root-coefficient  relations  of  (1-6)  are  transformed  into  other  five 
relations  which  are  equations  (1-12)  and  (1-14)  or  (1-13)  and  (1-15). 
For  convenience  (1-12)  and  (1-14)  are  repeated  here: 

f (  =  C0  +  C,(S,+5j  +  ^5,5, 

f2  =  ■  Co(S,  +  5,J  -f  C(S,S2 

In  actual  numerical  computations,  it  is  simpler  to  use  equation  (1-16) 
instead  of  (1-12)  and  (1-15),  because  of  the  progressive  nature  of  the 
root-coefficient  relations  of  algebraic  equations.   Namely,  compute  C0 
and  C^  first,  then  using  these  calculated  data  to  compute  Ct     and  then 

f,   and  f*  • 

In  general,  the  n  root-coefficient  relations  of  a  nth  order  character- 
istic equation  which  has  variables  in  coefficients  can  be  arranged  in  any 
form  as  one  wishes.   The  process  derived  here  has  certain  significances 
in  the  design  of  control  systems.   First  of  all9  the  order  of  the  character- 
istic equations  which  determines  all  the  roots  of  the  system  is  reduced  by 
the  number  of  the  variable  parameters.   Secondary,  the  physical  adjustable 
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parameters  are  expressed  directly  by  the  arbitrarily  chosen  roots.  The 
stability  and  response  of  the  system  is  determined  completely  by  the  re- 
duced characteristic  equation  while  the  adjustable  parameters  are  deter- 
mined directly  by  the  arbitrary  roots.   Since  equations  (1-16)  are  alge- 
braic equations  and  if  no  singularity  is  involved,  then  solutions  of  nt 
and  fa     which  actually  determine  the  physical  adjustable  parameters  al- 
ways exist.   Although  this  process  does  not  guarantee  the  adjustable  para- 
meter obtained  in  this  procedure  to  be  always  realizable,  the  realizable 
region  can  be  found.   In  the  design  of  a  control  system,  this  is  a  logical 
procedure,  namely  the  pole-zero  configuration  of  the  closed  loop  system  is 
chosen  first  from  the  specifications,  and  then  compute  compensator  para- 
meters directly  from  the  chosen  pole-zero  configurations.   For  every 
choice  of  pole-zero  configuration,  a  corresponding  set  of  values  of  the 
compensator  are  obtained.   This  partitioning  process  effectively  trans- 
forms the  desired  pole-zero  configuration  to  be  the  independent  variables 
of  the  design  procedure.   Moreover,  most  of  the  designs  are  based  on  a 
pair  of  dominant  complex  conjugate  roots,  then  two  of  the  arbitrary  roots 
can  be  chosen  (not  necessary)  as  the  dominant  roots.   The  adjustable  para- 
meters of  the  compensator  therefore  can  be  expressed  explicitly  as  func- 
tions of  the  dominant  roots  if  the  dominant  roots  are  defined. 

In  this  process,  there  are  no  limitations.   For  a  control  system  of 
any  order,  any  number  of  adjustable  parameters  and  any  type  of  compensa- 
tor; the  same  technique  can  be  applied. 
1-4.  General  Procedures  of  Partitioning  a  Controlled  Characteristic 

Equation. 

The  example  in  the  last  section  assumed  that  the  adjustable  parameters 
are  contained  in  two  coefficients.   Now  consider  a  cascade  single  section 


<l 


R-C  compensator  as  shown  in  the  block  diagram  of  Fig.  1-1,  in  which  the 
compensator  parameters  enter  all  the  coefficients.  The  characteristic 
equation  of  the  compensated  system  is 


54+  U,  +  PJ53+  (e><+Bxp)s*+(8,t>+*e>0)s-+l>B0 

where  B0  =  K  ,   B,  =  O-b   .       &L  =  CLi  b 


=  0   (1-17) 


Fig.  1-1  Single  section  cascade  compensator. 

(X  =  P/ £  '  ^  ,  a  and  b  are  know  numerically. 
Assume  -5  »  ~S  >  --Sj  an<*  ~~  5^  are  the  roots  of  (1-17).   By  expanding 
the  equation 

(S+S()(S-tSj(S453X5+-%)  =  O  (1-18) 

and  setting  the  corresponding  coefficients  of  (1-17)  and  (1-18)  equal, 
one  obtains  the  following  four  relations: 

B2+p  =  s,  +  5,  +S,  +-S,.  (i-19-i) 


B,+  B2^  =  (*,  +  ^K^  +  5«) +  -*,$*  +5^% 


(1-19=2) 


B,P4o(Ba=  5A(5,+^)  +  W5,  +  4J 


(1-19-3) 


PBo  =^=  ^A^ 


(1-19-4) 


Since  there  are  assumed  two  adjustable  parameters  0(  and  p    of  the  compensa- 
tor, there  are  two  arbitrary  roots.  Assume  -5  s  and  ~5  are  tne 
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arbitrary  roots,  then  (  S,  +  5^.  )  and  (   S  5^.   )  are  the  coefficients 
of  the  reduced  characteristic  equation.   Define 


C     —  5  4  S 


(1-20-1) 


<: 


o 


s3% 


then  the  reduced  characteristic  e 


5  +C.S  +  C 


ion  is 


O 


(1-20-2) 


(1-21) 


Substitute  the  definitions  of  (1-20)  into  (1-19),  one  obtains 


/ 

(5,+sj 
o 


0 

I 


■B, 

6. 


o 
o 

Bo 
o 


C, 

Co 

I3 


B,  -  SA 
O 
O 


(1-22) 


Equations  (1-22)  are  the  partitioned  four  root-coefficients  relations  in 
which  C.  and  C  define  the  reduced  characteristic  equation  while  K>  and  <X 
are  the  compensator  parameters.   These  four  dependent  variables  (C-,  C  s 
P  and  o{     )  can  be  solved  in  terms  of  the  arbitrary  roots  and  the  fixed 
plant  parameters  (  B  »  Bi  an<*  B  )•   For  numerical  computation  purposes, 
only  one  dependent  variable  need  be  expressed  explicitly  as  a  function  of 
the  arbitrary  roots,  the  others  are  computed  readily  by  the  previous  cal- 
culations.  Here  C  is  solved  from  (1-22)  as  shown  in  the  following  equa- 


tion 


where 


C0  =  ^U^-(Vs/-B,-tBz(s,+sj] 


(1-23-1) 


(1-23-2) 
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Now  consider  another  example  which  introduces  three  adjustable 
parameters.   Fig.  1-2  is  a  fourth  order  system  with  first  and  second 
derivative  feedback.   The  parameters  h. ,  h„  and  K  (=  Kii^O  are  assumed 


o 


i  _ 


K, 


O 


K. 


3(i+aJ(S+bJ(S+cJ 


is  4   4s- 


Fig.    1-2     Pure  derivative 


K,  h, 


and  h„,  are  variables  where 
-   ^-j  2 


adjustable.   The  characteristic  equation  is 


^+  B353-f  (B^-f  KzhJsZ+  (  B,+  K2h,)5  +  |C  =0  (1-24) 


where  Bj  =  a  +  b  +  C   ,   B^  =  &-b  +  be  +  CO  ,  and   B,  —  abC 
Assume  the  roots  are  -5,  ,  -£2  %    -  S3   and  -  5^.  then  the  root  coefficients 
relations  are  as  follows: 


Bi=  *|+^  +  *J  +  5« 


(l-25-l) 


3^+  K2h2=  (s,  +  5i)53-f  5,s,  +  (  s, -fia  +  S3)5, 


(1-25-2) 


3,+    fCjl,  =    5,5^3  +S+[5,5i+53C5,  +  5i)J 


K  =   5,  S2  S3  S, 


(1-25-3) 
(1-25-4) 
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Since  there  are  three  adjustable  parameters  (h  ,  h2  and  K) ,  three  roots 
are  arbitrary.   Let  these  arbitrary  roots  be  -  Sf    ,  -  S^  ,  and  -  Sj  . 
The  reduced  characteristic  equation  is  of  order  one  and  its  roots  is  -5^.. 

Define   S.  =  C  s  then  the  reduced  characteristic  equation  is 

5  +  C0^0  d-26) 

Substitute  the  definition  of  S,  =  C   into  (1-25) ,  the  root-coefficients 

4  o 

relations  become 

Co-=-    B3-(S,+  ^+S3j  (1-27-1) 

Kzh±  =  (V5,)53   +  S.S^-f  C0(5,+3a+53)-Ba  (1-27-2) 

Kih,=    ^^  +  C0[5(5i  +  (5,+  5J)Sj]  d-27-3) 

K    =     5>,$2^Co  (1-27-4) 

Equations  (1-27)  are  the  partitioned  root-coefficient  relations  in  which 
(1-27-1)  determine  the  system  response  while  the  other  three  equations 
determine  the  adjustable  parameters  so  that  to  give  the  roots  determined 
by  (1-27-1). 

From  the  three  examples  illustrated  aboves  it  can  be  seen  that  this 
partitioning  process  can  be  applied  to  systems  of  any  order ,  any  number  of 
adjustable  parameters  and  any  system  configurations.   The  important  concept 
of  this  process  is  to  separate  the  roots  of  the  controlled  characteristic 
equation  into  two  parts:  The  arbitrary  roots  and  the  constrained  roots. 
In  the  n  root-coefficient  relations,  the  arbitrary  roots  are  treated  as 
independent  variables  while  the  constrained  roots  and  the  adjustable  para- 
meters are  the  dependent  variables.   The  constrained  roots  are  obtained 
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from  the  solution  of  the  reduced  characteristic  equation  while  the  ad* 
jjustable  parameters  are  calculated  from  the  arbitrary  roots  directly. 
The  general  procedures  are  summarized  as  follows: 

(1)  Formulate  the  controlled  characteristic  equation*.   Set 
the  coefficients  equal  to  the  corresponding  root-coefficient. 

(2)  Define  the  arbitrary  roots+  and  the  coefficients  of  the 
reduced  characteristic  equation  (C  „  C....)«   Substitute  C  s  C. . . .  into 


the  equations  in 

(3)   Solve  for  the  aijustable  parameters  and  C°s  in  terms  of  the 
arbitrary  roots. 
1-5.  Summary. 

This  partitioning  process  effectively  transforms  the  independent  vari- 
ables from  the  adjustable  parameters  to  the  arbitrary  roots  of  the  character- 
istic equation.   As  the  adjustable  parameters  (compensator)  are  determined 
from  the  desired  closed  loop  pole-zero  configurations,  a  design  procedure 
for  control  systems  based  on  the  time  domain  can  be  derived.   The  parti- 
tioning procedures  presented  in  the  last  section  are  applicable  to  all 
system  configurations.   For  specific  cases  of  a  nth  order  systems  general  for- 
mulas can  be  derived.   For  the  cases  that  some  of  the  coefficients  themselves 

*In  discontinuous  operation  (2)s  the  controlled  characteristic  equations 
may  be  different  for  different  modes.  As  far  as  this  step  is  concerned, 
each  mode  must  be  treated  independently. 

+  The  number  of  arbitrary  roots  is  equal  to  the  number  of  the  independent 
adjustable  parameters.  The  it dependent  adjustable  parameters  may  be  less 
than  the  compensator  parameters  because  of  the  specifications  and  other 
restrictions.   For  details,  it  is  discussed  in  Chapter  IV. 
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are  treated  as  variables,  general  expressions  from  the  partitioning 
process  have  been  derived  as  shown  in  Table  1-1  to  Table  1-6.   For 
other  cases,  they  are  shown  in  the  subsequent  Chapters.   In  order  to 
illustrate  the  applicability,  the  following  examples  are  presented. 
Example:   Given  a  characteristic  equation 

5S"  +    !<*S%    |2855+5£oSZ+  \30OS  +  3O00  =  £>  (1-28) 

which  is  unstable.  Assume  a  network  can  be  devised  such  that  the  co- 
efficients of  the  first,  second  and  third  derivative  terms  can  be  ad- 
justed.  Then  the  controlled  characteristic  equation  is 

S*+l6>S*  +     f553  4  fzSljr   f,Si  300O  =  O  (1.29) 

where  j-      ,  -f,  and  f.    are  variables.   From  Table  1-5,  the  coefficients 
of  the  reduced  characteristic  equation  are: 

C,  =   B^-C  5(t5z+55J  (1-30-1) 

C0  =  80/s,S,Si  (1-30-2) 

The  functions  of  the  variable  coefficients  are: 

ft    =    CoCS^  +  S^+ijSj+CS^Sj  (1-30-3) 

fL=   CoiS^St+S^  +  CCS^+^Sj+SjS,)  +  S.**SS  (1-30-4) 

f 3  =    Co  +  G  (5,  +  52  +53)  4  (  s,5J  +  S*%   f  5,5, )  (1-30-5) 

where  (  -£j,  )>  (  ~32  )  and  (  -  Sj   )  are  the  arbitrarily  chosen  roots. 
Since  equation  (1-29)  is  of  the  5th  order,  one  root  must  be  real.   Assume 

15 


this  real  root  is  (  -S,  ).   For  convenience  let  (  -  i>,  )  and  (  -$2  )  be 
the  dominant  complex  roots  and  express  them  by  the  following  notations: 

5(5z  =  CxJn  (1-31-2) 

Substitute  (1-31-1),  (1-31-2),  B,  =  16,  and    gQ  =  3000  into  the  equa- 
tions (1-30),  one  obtains: 

C,  =  /G-UfUJn  +  S^)  (i-32-i) 

£o—    3000/a)j53  (1-32-2) 

f,    =    C0(ixJ^zfu)^3)  4  C,(Jn53  (1-32-3) 

fz    =    Cof^f^n  +  ^+C.^flfOn^)^    k£Sj  (1-32-4) 

f 3    ~     C0  +  C,(2J>u)rtfS5J  f(u>*  -f  ^tOnS^  (1-32-5) 


The  reduced  characteristic  equation  is 


(1-33) 


which  determines  all  the  roots  of  equation  (1-29)  and  consequently  deter- 
mines the  stability  boundary.   By  putting  all  the  three  arbitrary  roots 
on  the  left -half  plane,  then  from  (l-3<  I)     C   is  always  positive,  the 
condition  of  stability  for  the  arbitrary  roots  therefore  from  (1-33)  is 

c,  >o 

or  \(o     h.    2f60n  +^3  (1-33) 
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The  inequality  (1-33)  implies  that  the  three  arbitrary  chosen  roots  are 
confined  in  a  certain  region  if  the  system  is  to  be  stable.   This  region 
is  plotted  in  Fig.  1-3.   Notice   fU/n  is  the  real  part  of  the  complex 
roots.   If  S1  and  S?  are  complex  the  boundary  is  a  line  (  S  =  -8).   If 
S.   and  S?  are  real,  the  boundary  is  a  point,  (S  =  -8  +  j  O   ). 

Assume  the  dominant  roots  (  ~S,    )  and  (  -5^  )  are  desired  such 
that  all  other  roots  must  be  to  the  left  of  the  dominant  roots,  then 
the  following  inequalities  must  be  satisfied 

33  ^  $(Jti  (1-34-1) 

and  C,  >  2.fCOn  (1-34-2) 

Substitute  the  boundary  conditions  of  (1-34-1)  and  (1-34-2)  into  (1-32-1) 
one  obtains 

fCOn   £:  \(o/5  (1_35) 

Equation  (1-35)  defines  the  dominant  root  region  which  is  also  shown  in 
Fig.  1-3.   Within  this  region,  choose   ?  =0.5,  CJn=  4,  S„  =  4,  then 
from  (1-32-1)  and  (1-32-2) 

C,  =  8 
Co  ~*1 

The  roots  of  (1-33)  are:    -4  +   j  5.6  (  CO '  ■  6.85,  y'=  0.59) 
The  root  configuration  for  this  choice  of  the  arbitrary  roots  is  also 
shown  in  Fig.  1-3. 

The  next  step  of  the  design  is  to  calculate  the  variable  coefficients. 
Substitute  €  =  0.5,  6Jn  =4,  S  =  4,  C.  =  8  and  C  =  47  into  equations 
(1-32-3),  (1-32-4)  and  (1-32-5),  one  obtains 

f ,  =  2015 

fL    =  696 

f3  =  1*3 
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For  another  choice  of  the  arbitrary  roots,  the  coefficients  of  the 
reduced  characteristic  equation  and  the  variable  coefficients  are  com- 
puted in  the  same  way. 

The  above  example  assumed  that  three  coefficients  are  variable.   Other 
cases,  such  as  a  cascade  compensator,  feedback  compensator,  etc.,  can  be 
carried  out  in  the  same  procedure,  since  the  expressions  of  the  coeffici- 
ents of  the  reduced  characteristic  equation  and  the  adjustable  parameter 
are  obtained  readily.   The  stable  root  region  an'' dominant  oot  region  are 
important  in  design,  the  general  methods  for  evaluating  those  two  regions 
are  discussed  in  the  next  chapter. 
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boundary  of 
stability 


uouauary 
dominant 
roots 


-10 


5".  ,     M.     L  A, 


Fig.    1^3.      Example    for  5  +  \oS  +  f}S  +  f2S*+  f  S  +3000  =  0 
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order 

coeff.    of   the  reduced 

functions  of   the 

of   the 
equation 

chara.    equation 

variable   coeff,. 

Ct  =  1 

4-th 

<0=l 

C,=  Rf-C^Wj 

5ih 

C,  =  Bj-CiCS.ty-^ 

f,«.M*,tV) +?•*** 

/^    Bo 

fi«c;+Cl($l«J*4*Si 

Q-  1 

C3=  Bj.-.CS.  +  SJ 

b+h 

C«.2  =  1 

CW-3  =  BW.,  -Q.*0^ 

C,h>=:  Bn-i  TC.W<W)"Cm.% 

n+h 

Cr»-t  =  6*^3  —  Cw-4  ty^)~d-j*,S* 

• 

/-     0© 

0  ""?  3,;,. 

Table   1-1     Functions  of  the  variable   coefficients    f      and    -f,  Th* 

func   ions   contained   in  this  table  are  assumed  only  I Md  B^aW 'adjustable. 
Lll  /     2*   fre^e  arbitrarily,  chosen  roots  andVs  arl  the  coefficient; 

of  the  reduced  characteristic  equation. 


n 


S"+  C„.,5 


t-l 


—  -f    C,5  -f  C0  =0 
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order 
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.Functions  of  the 
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4  th 


£+h 


b* 


n+h 


c, 

Co 


Bo  /s,Sx 


c,  — 


-j 

c, 
c, 


=  I 


C3  = 


[  6,  -CoC^)]/^ 


c,  = 


BoAs2 


f2  —  d+A($?^+£V* 


Table  1-2.   Functions  of  the  variable  coefficients  f,  and  f3 
The  functions  contained  in  this  table  are  assumed  only  B  '  and  B  are 
adjustable   (-S^  and  <-s  )  are  the  arbitrarily  chosen  loots  a^d  c'i 
are  the  coefficients  of  the  reduced  characteristic  equation. 
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order 
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equation 

coeff.    of   the   redeced 
cnara.    equation 

Functions  of   tne 
variable  coeff. 

4-+*> 

• 

5"+h 

C,=  (  BZ-C0-C,S,SJ/(S,+S^ 
C0=     B./s,Sz 

6+h 

Ct==    6^  -C3(S,+S4)-  Si 5* 
C  =  6o/5,5i 

n^1 

C„.z  s=    1 

Cn-j  =  B^i  —  Cn-2  ($i^0 

...                      ♦ 

Ci    =    [   B,-  C0  -  Cx  S,S J  /£ 5|f50 
C6   ==    Bo /5,5Z 

Table  1-3   Functions  of  the  variable  coefficients  ~ft       and  f^ . 
The  functions  contained  in  this  table  are  assumed  only  B.  and  B  are 
adjustable.   (-S^  and  (-S2)  are  the  arbitrarily  chosen  roots  and  C's 

are  coefficients  of  the  reduced  characteristic  equation. 
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order 

coeff.    of   the   reduced 

Functions  of   the 

of   the 
equation 

chara.    equation 

variable   coeff. 

3«l 

C,  =    i 

to—   B^-O.+  S,) 

• 

Co=   Bz-C1(5,t5,)~Ci5,Si 

4«i 

C  -   Bi-Cz(VM 

C*«    1 

fo=  C.S.Si 

Co=   B,  ~($,+5,)C,  -S^C^ 

5* 

C,  =i  S}  -(^it)Cx-  s,szc} 

Q  =   ^-(S.+  so 

f{  =  cA+iJ+cs.s, 

Cj=    / 

C0=  B.  —(S.+OC,   -  i.^Ci 

C,  =:   Bj  —  Uffc)  Ci    —   V*Cj 

,                   ... 

,'-■ 

rm 

• 
Co-z-l 

Table  1-4.  Functions  of  the  variable  coefficients  f  and  f 
The  functions  in  this  table  are  assumed  B  and  B„  are  adjustable  ' 
simultaneously.  Where  (-S.)  and  (-S  )  ari  the  arbitrarily  chosen 
roots  and  C's  are  the  coefficients  of  the  reduced  characteristic  equation, 
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■      1 '.  ill 1 

order 

coeff.    of   the   reduced 

fun 

ictions   of   the 

of   the 
equation 

chara.    equation 

variabl 

e   coeff. 

C,=    1 

5ih 

c,=  e*-u+v^j 

■ 

Co=*   So/S^Si 

II 

II 

n 

(1       . 

c>  •■=  | 

d  SS    Sjr-(.S(i$L+$5) 

0 

-+■ 

o 

r-N 

bth 

C,  =  B^-Cats^+^-.CjCWtt*) 

p> 

-V- 

3* 

C0  =    B0 

Co  =    B./5,5tS3 

-4- 
-t- 

3" 

-4- 

Cm 

4- 
0* 

Cn-j=  1 

In 

4- 

4- 

•4" 

C^=:  Bn.,-Cn.3(S,+V^) 

O 

-+- 

£  rw  =•  B*-z  -  C^  ( Jfc«*«j)  -Qjte^ 

A 

n 

nth 

«            .           • 

3* 

-4- 

i 

•>* 

Table  1-5.   Functions  of  the  variable  coefficients  f.  f     and  -f 

The  functions  contained  in  this  table  are  assumed  B  ,  B  and  B  2*are 

adjustable  simultaneously.   Where  (-S  ),  <-s.)  and  \-SJ   are  tlL  arbitrarily 

chosen  roots  and  the  c's  are  the  coefficient!  of  the  rlduced  characteristic 
equation. 
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order 

coeff.    of   the   reduced 

fund 

;ions 

of   tne 

of   the 
equation 

chara.    equation! 

vara 

able 

coeff. 

4-+H 

c,  =  / 

• 

o 

^ 

M-n 

c0=  b,-  c,(5W*)-<  tt+tt+%) 

II 

II 

II 

$* 

Ci  =  B+-(5,+  Sz+5J 

Cz  =  1 

n 

o 

p 

*.                  i 

•-N 

^■N 

Va 

-V* 

i^ 

\r> 

iS* 

*V> 

1_ 

4- 

C0=  Bj  -c,(S,+V*j)  -  c*-(U+;^JfV.) 

£! 

C,  =  fc  -G  ($,♦*♦*)  -  C3  (^  +£&t$5J 

U 

-4- 

Ca  »■  %.-  c,  C s,+vy  -  C4  f5ft*«t5^) 

.         . 

*— s 

Vw 

• 

V/> 

4- 

n+h 

C0-J  =  1 

4- 

4- 
vO 

0* 

Table  1-6.   Functions  of  the  variable  coefficients  f 

unctions    COnf-flind    f«    «-V,,'„    *.«i_ii_ -  ..  '  /<»    s 


The  functions  contained  in  this  table  are  assumed  B  , 
adjustable  simultaneously.   Notations  are  the  same  $s   Tlble  1-5 


f ,  and  f 


and  B„  are 
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CHAPTER  II 
REDUCED  CHARACTERISTIC  EQUATION 
2-1.  Introduction:   The  reduced  characteristic  equation  determines  all  the 
roots  of  the  controlled  characteristic  equation  and  consequently  deter- 
mines the  response  of  the  system.   Usually  the  zeros  of  the  closed  loop 
system  are  known,  then  the  total  performance  of  the  system  depends  upon 
the  choice  of  the  poles  (roots  of  the  characteristic  equation).    In 
order  to  establish  a  guide  for  the  choice  of  the  roots,,  the  properties 
of  the  reduced  characteristic  equation  are  investigated.   For  the  reason 
of  simplicity  and  convenience  most  of  the  examples  in  this  chapter  take  the 
coefficients  themselves  as  variables.   For  the  case  of  the  adjustable  para- 
meters scattered  in  all  the  coefficients,  the  analyses  are  the  same.   In 
fact,  where  the  analysis  of  the  reduced  characteristic  equation  alone  is 
concerned,  it  does  not  make  any  difference  about  what  type  of  compensator 
is  used.   The  analyses  in  this  chapter  also  assumed  that  no  restrictions 
are  imposed  in  the  variables,  consequently  no  restrictions  are  imposed  on 
the  arbitrary  roots. 
2-2.  Coefficients  as  variables. 

Treating  some  of  the  coefficients  of  the  controlled  characteristic 
equation  as  variables  provides  a  convenient  and  relatively  systematic 
means  for  the  purposes  of  analyses.   However,  the  actual  applications 
have  some  difficulties.   Theoretically,  the  variables  can  be  transformed 
analytically  or  graphically  from  one  set  of  variables  to  another  by  the 
root-coefficient  relation  and  the  definition  of  variables.   Thus  the 
variables  which  originally  are  coefficients  of  a  characteristic  equation 
can  be  transformed  to  another  set  of  variables  such  as  the  actual  adjust- 
able parameters.   However,  the  transformation  may  be  very  complicated 
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because  of  the  inherent  nature  of  the  algebraic  equations.   In  actual 
design,  it  is  unnecessary,  because  the  actually  adjustable  parameters 
can  be  treated  as  variables  directly  by  the  partitioning  processes 
described  in  the  last  chapter. 

Mitrovic's  method  treats  two  coefficients  as  variables.   The  applica- 
tion of  Mitrovic's  method  to  the  case  in  which  the  actually  adjustable 
parameters  enter  some  of  the  presumed  fixed  coefficients  requires  that 

transformations  or  approximations  must  be  made  either  graphically  or 

2  4  5. 
analytically.   Some  of  the  work  to  this  end  has  been  done.  8  ' 

2-3.  Effectiveness  of  the  adjustable  coefficients. 

The  effectiveness  of  stabilizing  ability  of  adjustable  coefficients 
are  discussed  in  this  section.   Since  only  a  qualitative  analyses,  i.e., 
the  tendency  of  stabilization  is  required,  the  root-locus  method  is 
applied  to  the  reduced  characteristic  equation. 

Consider  a  fourth  order  characteristic  equation  (2-1]  the  effectiveness 

5*4  6353+  B^+B.S-f-  80— O  (2-i) 

of  adjustable  coefficients  are  analyzed  in  several  cases; 

Case  (1):   B«  adjustable.   Since  only  B«  is  adjustable,  then  the 

roots  are  constrained  by  B„,  B,  and  B  .   Assume  the  reduced  characterise 

3        1  o 


tic  equation   is 


S5+  C25Z-rC,S   4    C0  — -O  (2-2) 


and  the  arbitrarily  chosen  root  is  -S. ,  then 


B  =  Q  +  S, 


»!  -  C0-}-  C,S, 


Bo=  <-S' 
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From  the  above  equations,  solve  for  C, „  C„  and  C  ,    obtain 

1   2      o 

Q  =  6,  -5, 


c,  =  3;(b,-c0)  = 


_B.  _  6c 
J5,     3f 


(2-3-1) 

(2-3-2) 
(2-3-3) 


Substitute  (2-3)  into  equation  (2-2),  one  obtains 

5*  +  (  B,-  5,j5i  +  (■%)  5  -  -|S  (  S  -  SJ  =  O  (2-4) 

Equation  (2-4)  is  the  reduced  characteristic  equation  with  the  coeffici- 
ents expressed  in  terms  of  the  arbitrary  root  -S. .   Assume  S-  has  been 
chosen,  the  tendency  of  the  root  locus  of  equation  (2-4)  by  varying  B 
is  to  be  investigated.   Equation  (2-4)  in  root- locus  form  is 

-%  (  S-5,)        _ 


p. 


(2-4) 


Consider  equation  (2-4)  and  assume  S  ^f  0,  then  there  is  a  Eero  in  the 
right-half  plane  and  the  angle-criteria  of  the  root-locus  is  2-71H     (n  =  0, 
1,  2,)  instead  of  W7L    (m  =  1,  3,  5...).   The  general  root-locus  pattern 
is  shown  in  Fig.  2-1,  in  which  R  and  R„  are  the  roots  of  the  equation 


(A. 


Fig.  2-1  Root-locus  of  equation  (2-4)  Case  of 
B?  adjustable. 
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Case  (2):   B_  adjustable.   Now  consider  the  same  fourth  order 
characteristic  equation  (2-1)  in  which  only  the  coefficient  of  the  3rd  de- 
rivative term  is  adjustable.   By  the  same  manipulation  as  Case  (1),  the 
coefficients  of  the  reduced  characteristic  equations  are  obtained  as 
follows: 


c,  =.  B^St  -  Qo/s2- 


(2=5-1) 


c 


6,/s,  - 


b/5*  •+  &0/s; 


(2-5-2) 


c 


B0/S 


(2-5-3) 


Substitute  (2-5)  into  (2-2) ,  the  reduced  characteristic  equation  be- 


comes 


s5+(4:-  %)s*+&)s+-h(s'-s,s+s?j  =  o 


(2-6) 


In  root-locus  form  for  B  as  variables  (2-6)  becomes 


3.3 


(  sz-  s,s  +  sf  J 


S[^(|-|)s  +  |f] 


-I 


(2-7) 


Assume  S. 


0,  equation  (2-7),  has  two  complex  zeros  in  the 


right  half  plane  and  the  angle-criteria  of  the  root-locus  is  +  180' 
The  root  locus  pattern  is  shown  in  Fig.  2-2. 


o 


cv 


Fig.  2-2.   Root-locus  of  equation  (2-7)s  B„  adjustable 
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Case  (3):   B..  and  B.  adjustable. 


Consider  the  same  fourth  order  equation  but  with  B_  and  B-  adjustable. 


Since  there  are  two  coefficients  adjustable,,  the  reduced  characteristic 
equation  is  of  order  2,   Follow  the  same  derivation,  the  coefficients  of 
the  reduced  characteristic  equation  are: 


C.=  v^K-5A-c0) 


(2-8-1) 


C 


Bo 

5,5^ 


(2=8-2) 


where  -S.  and  -S?  are  arbitrarily  chosen  roots  and  S.  +  S^^^O,  S.S^^p  0 
are  assumed.   The  reduced  characteristic  equation  is 


s"+sfe(*-V--£)s+;§ 


o 


(2-9) 


In  root -locus  form  for  B  as  variable  is 

o 


bc 


WVS^>1 


5-   (5,+  s. 


>j 


(2-10) 


This  equation  again  has  one  zero  in  the  right  half  plane  and  the  angle- 
criteria  of  the  root-locus  is  2.7Lr\      (0  =  0,  1,  2,...).   The  general 
root-locus  pattern  is  shown  in  Fig.  2-3. 


Fig.  2-3.   Root-locus  of  equation  (2=10) 
B~  adjustable. 


B,  and 
1 
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Case  (4):   B?  and  B_  adjustable. 
In  this  case,  the  reduced  characteristic  e 
and  the  coefficients  from  Table  1-2  are  as  follows 

3i        5,  +  ia 


ion  is  of  order  2 


c,= 


3,5, 


(5,5.) 


-Be 


Co  = 


(2-11-1) 


(2-11-2) 


In  root- locus  form,  the  reduced  characteristic  equation  becomes 


'**  +  &) 


(2-12) 


Equation  (2-12)  has  one  zero  in  the  right  half  plane  and  the  angle- 
criteria  of  the  root-locus  is  ZJlf\       (n  =  0S  1,  2J,...).  The  general 

root- locus  pattern  is  shown  in  Fig.  2-4  when  S»  +  S,  ^=  0. 


-CK. 


Fig.  2-4.   Root-locus  of  equation  (2-12) 
adjustable. 


B2  and  B, 
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Case  (5):  B.    and  B?  adjustable  «  From  Table  1-1.   The  coefficient 
of  the  reduced  characteristic  equations  are  as  follows: 


C,  =  B3-  (5,  +  5.) 


(2-13-1) 


C0  = 


Bc 


3,3, 


(2-13-2) 


The  reduced  characteristic  equation  in  root-locus  form  is 


SS,   go 


i[s+  (63-5,-4)] 


=  -l 


(2-14) 


The  general  root -locus  is  shown  in  Fig.  2-5. 


^ 


+cv 


Fig.  2-5  Root-locus  of  equation  (2-14).   B„ 

and  B~  adjustable. 

Case  (6):   B.  adjustable.   In  the  same  way9  the  reduced  characteristic 
equation  in  root- locus  form  is  obtained  as  follows: 

Bo 


5/ 


-S[sN-(B,-SjS-KBa+5;*-  &>) 


=  -I 


(2=15) 


The  general  root-locus  pattern  is  shown  in  Fig.  2-6. 
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Fig.  2-6.   Root-locus  of 
adjustable. 


ion  (2-15)  B. 


Case  (7):   B. ,  B2  and  B~  adjustable: 

In  this  case  three  roots  are  arbitrary,,  only  one  root  is  constrained, 

Assume  S.  is  the  constrained  root,  and  define  C  =  S. «  then 
4  o    4 

The  reduced  characteristic  equation  is 

Bo 

If  no  restrictions  are  imposed  on  the  adjustable  coefficients,,  then  for 
any  choice  of  S„ ,  S„  and  S_  (S,^  0,  S2  ^F  0S  S-  =3p  0),  the  system  can  be 

stabilized  for  any  value  of  B  . 

o 

From  the  above  analyses  and  the  general  pattern  of  the  root- locus 
three  conclusions  can  be  drawn: 

(1)   If  the  adjustable  coefficients  are  not  in  the  lowest  order 
or  not  in  sequence,  as  the  cases  (1),  (2)s  (3)  and  (4)s  right  half  plane 
zeros  are  introduced  in  the  root-locus  plot.   Consequently,,  the  stabiliz- 
ing ability  is  less  than  that  if  no  right  half-plane  zeros  would  be 
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introduced. 

(2)  If  the  adjustable  coefficients  are  in  the  lowest  order 
and  in  sequence,,  as  the  cases  (5),  (6)  and  (7)9  no  right  half-plane 
zeros  are  introduced.   Consequently  the  stabilization  is  more  effective. 

(3)  The  most  effective  scheme  is  the  case  where  all  of  the 
coefficients  are  adjustable  as  case  (7).  The  most  ineffective  scheme 

is  the  case  where  only  the  highest  coefficient  is  adjustable'  as  case  (1). 
In  generals  if  there  are  "V  unadjustable  coefficients  which  have  an 
order  lower  than  the  highest  order  adjustable  coefficients  then  the 
reduced  characteristic  equation  in  root -locus  plot  for  B  as  variable 
has  "r"  zeros  in  the  right  half  plane.   Because  of  this  property ,  it  seems 
that  if  the  adjustable  coefficients  of  a  characteristic  equation  are  not 
the  lower  order  and  not  in  sequence,  the  stabilizing  tendency  of  these  sets 
is  less  than  for  the  case  of  adjustable  coefficients  which  are  of  lower 
order  and  in  sequence.   Consider  case  (5),  in  which  B„  and  B~  are  adjust- 
able, the  reduced  characteristic  equation  is  of  the  second  order  and  there 
are  no  zeros  in  the  right  half  plane,  consequently  the  system  can  always 
be  stabilized  for  any  value  of  K.   However^  for  cases  (1)  to  (4),  it 
cannot  always  stabilize  the  system.   In  general 9  for  any  order  system 
if  the  reduced  characteristic  equation  is  of  second  order  and  all  adjust- 
able coefficients  are  of  lower  order  and  in  sequences 9  then  the  system  can 
be  stabilized  for  any  forward  gain  K.   However9  this  is  not  the  case  if 
the  adjustable  coefficients  are  not  of  lower  order  and  in  sequence.   More- 
over, the  more  zeros  introduced  in  the  right  half  plane  in  the  process  of 
derivation  of  the  reduced  characteristic  equation,  the  less  the  tendency 
to  stabilize  the  system.   This  implies  that  the  more  the  unadjustable  co- 
efficients which  are  of  lower  order  than  the  highest  order  adjustable 
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coefficient ,  the  less  the  effectiveness  of  the  control  parameters.   This 
is  the  inherent  property  of  algebraic  equations.   It  is  because  of  this 
property  that  many  approximations  can  be  made  by  ignoring  those  control 
parameters  which  affect  only  the  coefficients  of  higher  derivatives.   If 
the  design  of  a  control  system  is  based  on  the  characteristic  equation, 
those  coefficients  of  lower  order  derivatives  should  be  considered  first 
if  it  is  possible. 

This  analyses  also  explains  the  well  known  fact  that  a  cascade  filter 
is  more  effective  than  lower  order  derivative  feedback,  especially  for 
high  order  systems.   For  a  cascade  filter,  the  filter  parameters  are 
scattered  in  all  coefficients  of  the  characteristic  equation  in  the  nature 
of  algebraic  equations,  while  derivative  feedback  enters  some  of  the  co- 
efficients only.   This  analyses  also  explains  the  fact  that  tachometer 
feedback  is  good  for  low  order  systems,  but  may  even  hurt  the  performance 
of  high  order  systems. 
2-4.  Root-Region  of  Stability. 

The  reduced  characteristic  equation  is  .a  function  of  the  arbitrary 
roots.   Mathematically  for  any  choice  of  the  arbitrary  roots,  there  ex- 
ists solution  of  the  equation.   However,  in  control  systems,  the  roots 
of  the  characteristic  equations  are  required  to  be  in  the  left  half  plane, 
this  requirement  therefore  confines  the  arbitrary  roots  in  a  certain  region 
in  the  left-half  plane.   This  region  here  is  called  "region  of  stability" 
and  may  serve  as  a  guide  for  the  choice  of  the  arbitrary  roots. 

If  the  number  of  variables  is  one,  this  region  is  a  line  or  lines 
in  the  S-plane  (root-locus).   When  the  number  of  variables  is  more  than  one, 
this  region  is  an  area  or  several  areas  on  S=plane.  The  boundary  of  this 
region  can  be  calculated  by  Routh's  criterion  subject  to  the  reduced 
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characteristic  equation.   The  pattern  of  4th  and  5th  order  equations  for 

f~,   -  fz        j       fi~   T3   and   "f",  -  f$  as  variables  are  shown  in  Fig.  2-7 
and  Fig.  2-8. 
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Fig.  2-7  Root  region  of  a  4th  order  characteristic  equation 
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Fig.  2-8  Root  region  of  a  5th  order  equation 


)<<J 


(c)   f/,f3 

variable 


38 


In  Fig-  2-7  and  Fig.  2-8,,  if  a  pair  of  complex  roots  are  in  the 
stable  region^  then  all  the  roots  are  in  that  region.   If  a  pair  of 
complex  roots  are  in  the  unstable  region^  then  there  is  at  least  one  root 
in  the  right  half  plane.   If  a  pair  of  complex  roots  are  on  the  stability 
boundary  line,  then  there  are  at  least  one  pair  of  complex  roots  on  the 

j  CO  axis  or  one  root  at  the  origin.   The  arbitrarily  chosen  roots 8 
therefore^  must  be  confined  in  the  stable  region.   The  calculation  of 
the  boundary  is  illustrated  in  the  following  example. 
Example  2-1:   Given  a  fifth  order  equation 

<f+    \<oS*+  I23 s3  t  £*osz-h  i loo 3    4  2OO0  =0  (2-16) 

Case  1:  Assume  T"  and  j~     are  variables9  then  the  controlled 
characteristic  equation  becomes 


5*  +  l(~S*  +   f35*  +  "hs2-!-  Iloos   +  Zooo  —O 


From  Table  1-2 


Cz   =  3+-   (5,4  5*  J  (2-18-0 


6, 


■o  S,5Z 


(2-18-3) 


where  -Si  and  -S9  are  the  arbitrarily  chosen  roots.   Assume  these  roots 
are  confined  on  a  constant   ?  line8  then 
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Substitute   into    (2-18) 

C2  =    Bfr-^Wn 


(2-19-1) 


£& 


(2-19-2) 


Co  — 


_6o 


(2-19-3) 


where  B,  =  16,  B  =  2000,  B.  =  1300  from  the  given  equation  (2-17). 
The  Routh's  criterion  for  the  3rd  order  reduced  characteristic 


equation  is 


C,Q  >C 


(2-20) 


Substitute  (2-19)  and  the  values  of  B1 ,  B,  and  B  into  (2-20), 
and  manipulate,  the  Routh's  criterion  in  equality  becomes 

2.4  j^n  +  ~^-  <   18.8  *af  (2.21) 


Let  the  left  side  of  the  (2-12)  be  4-    ,  then 


UJ 


The  +  versus  OOn    plot  for  given  value  of  P  is  shown  in  Fig.  2-9. 


Fig.  2-9  Graphical  solution  of  equation  (2-21) 
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From  the  Inequality  (2-21),  if     CaJ,  £:  0Un  —  ^z  ,  then  all  the  roots 
of  the  reduced  characteristic  equation  are  on  the  left-half  plane.  COt 
and  U)    define  the  stability  boundary  of  the  arbitrarily  chosen  roots  and 
are  the  solution  of  the  following  equation 

fa.  f  2. 

Z.(o  f*>n   "+  -~£  =  \3.S   +  M  (2-22) 

Solutions  of  (2-22)  for      O  ^  J  *^  |   are  tabulated  in  Table  2-1,  and 
the  boundary  for  root  region  are  plotted  in  Fig.  2-10. 

Table  2-1.   Solutions  of  equation  (2-22) 

1  3.7  6.62 

0.9  3.28  7.53 

0.7  2.4  10.5 

0.6  2.05  12.9 

0.5  1.7  14.3 

0.4  1.35  17.9 

0.3  ^05  23.95 

0.2  0.65  36.1 

0.1  0.35  72 


Case  2 :  f ,  f  variable 

Since  the  coefficients  of  the  reduced  characteristic  equation  are 
expressed  by  the  arbitrarily  chosen  roots,  the  mapping  contour  can  be 
chosen  arbitrarily.   In  order  to  demonstrate  this,  let  the  two  arbi- 
trarily chosen  roots  be  on  a  constant  Online,  then 


5(+5^  =  lo^  (2-23-1) 

S)5z   =  ot  +  UU*  (2-23-2) 
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Fig.  2-10  Stable  root  region  of  a  5th  order  equation  with 
variable  of  example  2-1. 
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Substitute  (2-23)  into  the  coefficients  of  the  reduced  characteristic 
equation  in  Table  1-2,  one  obtains 

Cz   =    It,  -  2<?v  (2-24-1) 


i  .    .  .-* 


C,  =    £ZO  -(  Ko-Z^JZ^   -  (CK+LOj  (2-24-2) 

C0  =    £000//(^+  COaJ  (2-24=3) 

The  Routh's  criterion  is    C^C,    3  Co   Substitute  (2=24)  into  the 
above  inequality  and  manipulates  obtain 

(  |^-2^J  J  |28 -(It-^^-j^^  -(a  +U;^)J(^f  CU^J  >  IWO  (2-25) 

The  boundary  value  of  CaJ  for  the  inequality  (2-25)  are  solved  by  the 

same  procedure  as  Case  1  and  tabulated  in  Table  2-2.   The  stable  root  region 

is  plotted  in  Fig.  2-11. 

Table  2-2  Stable  root  region  of  a  5th  order  equation 
with  f,  and  f      variable, 

q^  Stable  region  of  u_) 

0  0.975  <  CJ  <  11.29      -0.975  >  U)  >  -11.29 

2  -8.57  <  uO  <   +8.57 

4  -6.62  ^  ^  <  +6.62 

6  -5.9  <  uJ  <  +5.9 

7  -6.3  <  uj  <  +6.3 
7.5  -5.9  <:  ou  <^  +6,9 

Case  3:  4-       f,  variable. 

In  this  case,  from  Table  1-3,  the  coefficients  of  the  reduced  character- 
istic equation  are: 

Cz   —  a.-C^+V)  (2-26-1) 
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Fig.  2-11. 


Stable  root-region  of  a  5th  order 
equation  with  -f,  and  f  x  variable 
in  example  2-1. 
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C0  =  B0/-S,Si  (2-26-3) 

Assume  S.  and  S?  are  on  a  constant   f  line.   Then 

St  +  Si  =  */<*>n   ,      5  ^°  (2-27"1} 

S(5z  =  CO*  (2-27-2) 

Substitute  (2-26)  and  (2-27)  into  the  Routh  criterion  inequality  and 
manipulate,  obtain 

(|b-^.i)(JWO  -  *2°£  -,bU)l  +  ^f^)U„  ==>  aj  «  looo        (2-28> 

This  inequality  is  a  function  of  OJ  up  to  the  5th  order9  the  solution 
is  more  easily  found  by  a  graphic  method.   The  function  of  the  left  side 
of  (2-28)   for   P   as  a  parameter  are  plotted  as  a  function  of  <aJ n    , 
boundary  values  of  C<Jn  are  evaluated  as  shown  in  Table  2-3.   The  root- 
region  of  stability  is  plotted  in  Fig.  2-12. 

£  5+able    region 

0.1  2.2     <    (jOn  <5.45 

0.3  2.3    <:     U0n  <Z  6 

0.5  2.4    <r    UUn-7.1 

0.7  2.5    <C    UUn<10.7 

0.9  2.6    <    UJn<8.6 

1.0  2.68  <    UJn<:7.3 

Table  2-3.   Stable  region  of  5th  equation  with 

r    ,   .L    variable. 
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Figure  2-12.   Stable  root-region  of  the  5th  order  equation  with 
i     ,  4   variable  in  Example  2-1. 
Ji   JJ 
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In  equation  (2-26-2),    5>  +  5  ifO  was  assumed.   This  implies 
j  ip  0.   For  <?=(),  the  coefficient  C  of  the  reduced  characteristic 
equation  is  undefined.   Rearrange  equation  (2-26-2) 

C,(5,+52J  =■  82~C0  -C,s,52 

For  f  =  0,     5(  +  5^  ■  0   then 

consider  equation  (2-28; ,  the  coefficients  B,s  B«  and  B  are  unadjustable 

coefficients  of  the  given  characteristic  equation,,  therefores  the  solu- 

/  ... 

tion  of  (2-28)  defines  four  points    tjU),   ,  and  Ij^z     on  t*ie 

axis  instead  of  a  region.   Those  four  points  are  convergent  points  for  the 
stability  boundary,  i.e.,  if  there  are  any  roots  on  the  j  UJ   axis  for  this 
type  of  control,  the  roots  must  be  either  ijw,  s  or  +  j  txJ,  regardless 
of  where  the  other  roots  are.   Consider  Fig.  2-13  which  is  the  same  as 
Fig.  2-12,  but  repeated  for  convenience.   By  the  definition  of  stability 
boundary,  if  there  are  any  roots  on  the  boundary,  there  are  at  least  one 
pair  of  roots  on  the  jcJ  axis  or  a  root  at  the  origin.   Assume  a  point 
P  (inside  the  stable  region  or  the  unstable  region)  is  a  root  of  the  equa- 
tion, as  P  approaches  the  boundary  in  any  manner  whatever,,  equation  (2-28) 
implies  there  are  either  a  pair  of  complex  roots  which  approach  the  points 

i  j  6j,  or  Ij^i  j,  or  a  root  which  approaches  the  origin.   In  case  one 
and  case  two  of  this  example,  the  roots  on  the  j UJ   axis  have  an  interval, 
and  the  location  of  the  roots  in  this  interval  depends  upon  where  the  other 
roots  are.   In  this  case,  the  roots  on  the  j  CU  axis  are  independent  of  the 
other  roots.   In  general,  if  the  variable  coefficients  are  not  in  sequence, 
then  the  roots  on  the   i  CO    axis  have  this  property.   In  short,  for  this 
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Fig.    2-13     Convergence  on  the    j  6j   axis 
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u^  =  S.Z 


type  of  control  if  the  variable  coefficients  are  chosen  such  that 

there  are  roots  on  the   j  CO    axis,  then  those  roots  on  the  j  uj    axis  are 

completely  determined  by  the  fixed  coefficients  of  the  given  equation. 

In  this  example,  substitute  B  ,  B  and  B,  into  equation  (2-28) , 
one  obtains 

ou,   -  z.w 

The  discussions  above  are  for  complex  roots.   For  real  roots  the 
analysis  can  be  carried  out  just  by  the  same  procedure. 
2-5.  Dominant  Root  Regions   No  Zeros  Close  to  the  Origin. 

Since  most  of  design  is  based  upon  a  pair  of  dominant  complex  rootss 
a  procedure  of  evaluating  the  "dominant  root  region"  is  to  be  discussed. 
The  definition  of  dominant  roots  depends  upon  the  degree  of  dominance  and 
other  requirements.   It  can  be  defined  as  one  wants.   In  case  of  no  zeros 
close  to  the  origin  and  if  the  pair  of  dominant  complex  roots  are  so  de- 
fined such  that  the  negative  real  part  of  all  other  roots  (complex  and 
real)  of  the  original  characteristic  equation  must  be  greater  than  that 
of  the  presumably  dominant  roots,  then  the  Routh's  criterion  can  be 
applied  to  the  reduced  characteristic  equation  by  a  conformal  mapping 
process, 

/  jV 

/ 


/ 

-7" 

/ 

/ 

/ 
/ 
/ 


^ 


(a)  S-plane 


p-plane 


Fig.  2-14  Conformal  mapping  of  S  =  p- 
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Consider  Fig,  2-14(a),  assume  the  pair  of  arbitrary  complex  roots 
(-S  )  and  (-S«)  are  the  dominant  roots  as  defined  in  the  last  paragraph 
and  also  assume  all  the  other  roots  of  the  original  characteristic  equa- 
tion must  have  the  real  parts  which  are  greater  than  "a"  (a  is  positive 
real  number  and  greater  than  or  equal  to  f(^n   )»  then  all  the  roots  of 
the  reduced  characteristic  equation  must  be  to  the  left  of  the  line  S  =  -a 
in  the  S-plane.   By  the  linear  transformation  S  -  p-a,,  the  region  to  the 
left  of  the  line  S  =  -a  in  the  S-plane  is  transformed  to  the  left -half  plane 
in  the  P-plane.   Therefore  the  Routh's  criterion  applied  to  the  p-plane, 
implies  all  the  roots  of  the  reduced  characteristic  equation  are  to  the 
left  of  the  line  S  =  -a  in  the  S-plane.   In  addition,,  to  put  all  other 
arbitrarily  chosen  roots  (there  may  be  more  than  two  arbitrarily  chosen 
roots)  also  to  the  left  of  the  line  S  -  -a,  then  all  the  roots  of  the 
original  characteristic  equation  except  a  pair  of  complex  dominant  roots 
are  to  the  left  of  the  line  S  ■  -a.   Here  the  line  S  =  -a  was  assumed  to 
be  fixed  if  the  dominant  roots  are  fixed.   In  order  to  find  the  dominant 
root-region  which  satisfy  the  definition,  'a'  must  be  a  variable  as  the 
dominant  roots  are  varying.   Therefore  the  application  of  Routh's  criter- 
ion to  the  p-plane  defines  a  region  in  the  S-plane  if  such  a  region  exists. 

To  obtain  the  transformation  from  S-plane  to  p-plane,  the  process  is 
just  a  substitution  of  S  =  p-a  into  the  reduced  characteristic  equation. 
Consider  a  third  order  reduced  characteristic  equation 

S3+  ClSZ+   C.SIQ^O  (2-29) 

substitute  S  -  p-a  into  the  above  equation,  one  obtains 
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manipulating,  it  becomes 

p5+  U1-3a.)Pz+   (  C,  ~^aq  fia!jp  +LC0-cicL+Cla:-o?)    (2"3°) 

=  0 

Equation  (2-30)  is  the  transformed  reduced  characteristic  equation  in  the 
p-plane.   Now  define: 

J£  ^  Cz  -3a.  (2-3i-i) 

J>t  =  C,-  23Q  -f3o.2  (2-31-2) 

I>o  =  C0  -  C.a.  +  C*af-  a.3  (2-31-3) 

Where  D's  are  the  coefficients  of  the  reduced  characteristic  equation  in 
p-plane,  then  equation  (2-30)  becomes 

33  +  V>SX+   J>,S  +  _D0  =  O  (2-32) 

The  Routh  criterion  of  equation  (2-32)  is 

J>,  D,  =>  J>0  (2'33) 

As  the  coefficients  D's  are  functions  of  the  arbitrarily  chosen  roots  and 
the  definition  of  'a',  the  inequality  from  Routh's  criterion  defines  a 
region  on  the  S-plane  directly. 

The  coefficients  D's  for  a  nth  order  reduced  characteristic  equation 
and  its  Routh' s  criterion  inequality  have  been  formulated  and  tabulated 
in  Table  2-4.   Since  the  reduced  characteristic  equation  of  any  order 
system  and  any  number  of  variables  can  be  deriveds  the  process  of  evalua- 
ting the  dominant  root-region  for  the  definition  described  above  can  be 
applied  to  any  order  system  and  any  number  of  variables. 
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Table  2-4,   Transformed  coefficients  D"s  of  the  reduced  characteristic 
equation  "a"  is  the  mapping  contour,  C's  are  the  coefficient  in  the  S-plane 


order  of   the 
reduced 
equation 

transformedcoeff . fs 

Rough's 
criterion 

Zr6 

J>,  =:   C,  ~2cX 

P>  =  c;-c,a.-fa?" 

3>(>o 

3rd 

p>  =  1 

Pr  =  CA-3a. 

D,  =  c, -2a.cA-J3at 

P0  =  Co-c.a  4  a'c1-^ 

P,>o 
J>A  >J>. 

^th 

p«  =  1 

P,— ;c,-  20.Q,  +3<£Cj~40- 

> 

! 

JVi  =    Cn-i  -    naCn 

Pn-2  =  Cp.2  -Crhija(^|+--— -aGi 

12-3 

-P„-t  -    C  n+  -c  n-3)  a  C-3  t  ^jt^  ^ 

*  ■—rr^rclCn 

• 

i 
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Fig,    2-15     Plot   of   inequality  2-37, 
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Example  2-2.   Consider  the  same  characteristic  equation  given  in 
Example  2-1,  namely  - 

55"  4  3^   +     &jS    ■+  faS*  +  f,S    +  80  =  O  (2-34) 

where  B,  =  16,  B0  =  128s  and  B  =  2000,   The  coefficients  of  the  reduced 
4       3  o 

characteristic  equations  are: 

CL    =    [io-  If^n  (2-35-1) 

C  ,   =     lZ&-3Z.LJn    +  (  4f "  'J  ^n  (2-35-2) 

Co=  *ooo/*Jn*  <2-35~3> 

where  V  and  uJ    define  the  dominant  roots,  since  there  are  only  two 
arbitrarily  chosen  roots.   Assume  a  =  P U-)  ,  i.e.,  all  the  other  roots  must 
have  real  parts  greater  than  or  equal  to  ?^„  »   Substitute  equation  (2-35) 
into  (2-31)  the  transformed  coefficients  of  the  reduced  characteristic 
equation  become: 

P,  =2    l1p-S-fiA)n  (2-36-1) 

Jpi   =    \2&-<pffU)n+(nf-0tOn  (2-36=2) 

J)o  =  looo/u*   -  I2&  fU)n   H6f  o£  -  J»  (  if-  0  ^,t  (2  =  36-3) 

By  substitution  of  (2-36)  into  the  Routh's  criterion  B^  >  Dq  and  manipul- 
ating, one  obtains: 

Let  the  left  side  be  4  (fed),  choose  ?  and  plot  -f  versus  CUn  as  shown 
in  Fig.  2-15,  the  dominant  roots  (defined  by  a  =  fiJj   region  which  satisfy 
the  definition  for  P  =0.3,  0.5,  0.7  are  evaluated  from  Fig.  2-15  as 
follows: 
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f  -  0.3  1.4  ±    6jn  1=  6.3 

f  -  0.5  1.65  £:  U>n  £;  5 

j5  -  0.7  1.9  £  wJn  *=.   4.3 

This  region  is  plotted  in  the  S-plane  as  shown  in  Pig.  2-16;;  in  which 

the  stability  boundary  Is  also  plotted. 

In  the  dominant  root-region,,  assume  the  dominant  roots  are  chosen  to 

be  -S1,  -S2  =1.89  +  j  6;   {   f  =  0.3,  64,-6.3) 

the  secondary  roots  evaluated  from  the  reduced  characteristic  equation 

are : 

-S3,  -S4  =  -2.8   +   j  1.0S    (   J-  0.94s    Cj^  2.82) 

-S_  -  -6.6 
5 

It  shows  that  the  secondary  roots  have  real  parts  which  are  greater 
than  that  of  the  dominant  roots. 

If  the  dominant  roots  are  chosen  as 
-S1,  -S2  ■  2.5   +    j  4.329     (   ^=0.5,  U>n  =  5) 
the  secondary  roots  are: 
-S39  -S4  =   3.56  +    j   2.76,    (  £  -  0.79,   UH-  4.5) 


■S.  *  -3.88 

5 


It  shows  again  the  real  parts  of  the  secondary  roots  are  greater 
that  of  the  dominant  roots. 

In  the  above  two  choices  of  the  dominant  rootsy  all  the  secondary 
roots  also  have  greater  damping  ratio  ?    „   However  the  dominant  root  region 
evaluated  by  this  procedure  doesn't  imply  this  restriction  to  the  secondary 
roots,  but  this  restriction  can  be  obtained  by  choosing  large  value  of  "a" 
relative  to  P  ^n  and  by  inspection  of  the  shape  of  the  region. 

In  Fig.  2-16,  It  can  be  seen  that  the  dominant  root-region  is  inside 
the  stability  root-region.   The  secondary  roots  are  not  necessarily  inside 
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the  stability  boundary.   The  max.  damping  of  the  dominant  roots  can 

also  be  evaluated  from  the  dominant  root -region. 

Example  2-3.   Given  a  6th  order  characteristic  equation 
5%  %5r  +  3^  i    f3  S3  4  f25*  4  f ,  5  4  60  -  O 

Q 

where  f   f    -f _   are  adjustable,  Bc  ■  184,  B.  =  10840  and  B  =  3x10  . 
'  i  j   '  i    *    '  3        J       s   5        4  o 

Here  three  coefficients  are  adjustable,  three  roots  can  be  chosen 
arbitrarily.   Assume  the  requirements  of  the  dominant  roots  is  such  that 
all  the  secondary  roots  must  have  the  real  negative  parts  at  least  twice 
that  of  the  dominant  roots,  then  the  mapping  contour  is  a  =  Zf^n   .   In 
order  to  meet  the  requirements  of  the  secondary  root,  the  third  arbitrarily 
chosen  root  must  be  set  at  least  to  the  min.  value,  namely  S.  ~2f(A-''*  . 
Assume  S«  =  2fto\,   then  the  three  arbitrarily  chosen  roots  are: 


-5,  ,  -32  =  -  fU)n    1  J  UAJ  I-  J* dominant 

-  s3  —    -2fuJn 

From  Table   1-5,    the  coefficients  of  the   reduced  characteristic  equation 
are: 

Ci=     ^-(S,4  52+S3j  (2  =  37-1) 

C,==     5|.-CiCSl+5Jl4S,)-(S,5,+5|SJ+5|Sj  (2-37-2) 


C0   =     &o/s,5zS3 
Substitute   S. ,    S2   and  S      into    (2-37)   and  manipulate,    obtain 


(2-37-3) 


C1  =    B^-4-f^r,  (2-38-1) 

C,=    B4.-4j'BP^»  +('^2-'J^n  (2-38-2) 

C0=    Bl/2fM,3  (2-38-3) 
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Fig.  2-16. 
Shaded  area 


Dominant  root  region  of 
dominant  root  region 


a  5th  order  equation  with  f, ,   f2  variable. 
Dash  line :   Boundary  of  stability. 
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rig.  ^-i/«   Dominant  kooc  Kegions  ox  example  ^-J. 


Substitute  (2-38)  into  the  coefficients  of  the  transformed  reduced  character- 
istic equation  in  Table  2-4,  obtain 

P,   -  %-  SJB^,,  +  (4of-i)Lj* 

The  Routh*s  criterion  D~D,  >  D  after  manipulation  becomes 

^  1    o 

Substitute  the  numerical  values  of  Be»  B,  and  B  .  it  becomes 

5   4      o 

u)tl  ;[5t^-^(l;;Joo\4  mUol$X'0»*-lfW1'-OUn\  >U^f    (2-39) 

For   P  =0.3,  0.5,  0.7  and  0.9,  the  region  of  U)     which  satisfy  the  in- 
equality (2-39)  are  evaluated  and  plotted  as  shown  in  Fig.  2-17, 

<^   =0.3  7.5  fr«JM  ±30 

p   =  0.5  7   ^  ^^19.5 

&     =  0.7  6.8  ^  tj^O 

P   =0.9  None 

Within  the  dominant  root  region,  if  the  dominant  roots  are  chosen  as: 
_Sr  -S2  =  -9.5   +    j   16.4;   (  f   =  0.5,   cJn=  19) 
then:    -S  =  -2  f  oOn   =  -19 

The  other  three  roots  evaluated  from  the  reduced  characteristic 
equation  are : 
-S4,  -S5  ==  -20  +    j25;    (  0   =  0.97,  U^  =  20.6) 

-S6  =     =106 
If  the  dominant  roots  are  chosen  as 
-Sr  -S2  =  -8   +      J13.8,    (   f  =  0.5,  u)n=   16) 
The  secondary  roots  are 
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-S3  =  -16 

-S4,  -S5  =  -23  +    j  13,    (  f   =  0.87,  Cl>„  =  26.4) 

-S,  =  -106 

The  evaluation  of  the  dominant  root  region  is  not  necessarily  carried 
out  point  by  point,  only  those  points  of  interest  are  evaluated.  And 
furthermore,  only  the  boundary  points  which  satisfy  the  Routh's  criterion 
are  of  significance.   The  Routh's  criterion  formed  in  this  process  is  a 
polynomial  of   lAJn  in  any  case,  therefore  for  higher  order  reduced 
characteristic  equations  there  may  be  more  than  one  region  which  satisfies 
the  criteria.   But  only  the  region  which  is  closest  to  the  joJ  axis  is  of 
interest,  because  the  roots  are  assumed  to  be  dominant. 
2-6.  Modified  Dominant  Root  Region  -  One  Zero  Close  to  the  Origin. 

The  assumption  that  the  response  of  a  system  can  be  characterized  by 
a  pair  of  complex  dominant  roots  is  based  on  the  approximations  that  all 
other  roots  are  either  too  far  to  the  presumable  dominant  roots  or  close 
to  a  zero.   Consider  Fig.  2-18,  the  roots  (-S.)  and  (-S~)  are  said  to  be 
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Fig.  2-18   General  Closed  Loop  Pole-zero  Configuration 
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because  the  other  roots  (-Sc),  (-S,)  and  (-S„)  are  far  to  the  left 

DO  J 

and  (-S,)  is  close  to  a  zero  (-Z.)  so  that  they  contribute  negligibly  to 
the  transient  response.   If  no  zeros  lie  to  the  right  of  the  dominant 
roots  (close  to  the  origin) s  then  all  the  secondary  roots  are  required  to 
be  far  to  the  left  if  the  presumed  dominant  roots  are  dominant.   If  there 
is  a  zero  close  to  the  origin,  then  one  root  to  the  right  of  the  domin- 
ant roots  has  negligible  effect  on  the  dominance  if  that  root  is  close  enough 
to  the  zero.   Therefore,  the  dominant-root  region  defined  in  the  last  sec- 
tion must  be  modified  when  there  are  zeros  close  to  the  origin.   Usually, 
the  zeros  of  a  system  are  known.   If  the  known  zero  is  far  to  the  left  of 
the  dominant  roots,  then  the  dominant  root  region  defined  in  the  last  sec- 
tion  is  still  valid.   When  the  zero  is  known  to  be  close  to  the  origin,, 
the  modification  of  the  dominant  root  region  is  required. 

Consider  Fig.  2-19(a)  which  is  a  third  order  system  with  a  cascade 
compensator.   Assume  the  compensator  is  a  lag  network.   Let  (-S..)  and 
(-S^)  be  the  arbitrary  roots  (also  the  dominant  roots  in  this  case),  then 
(Z)  and  (p)  are  functions  of  S  and  S„,  (Equation  1-22).   The  general  pattern 
of  the  root  locus  is  shown  in  Fig.  2-19(b).   Since  S.  and  S~  are  arbitrary 
roots,  then  from  the  tendency  of  the  root  locus,  the  one  only  root  which  is 
possible  close  to  the  origin  is  (-S_).   If  the  arbitrary  roots  are  chosen 
within  the  dominant  root  region  (point  B  in  Fig.  2-19E)  defined  in  the  last 
section  and  the  pole  (p)  and  zero  (Z)  of  the  compensator  are  adjusted  ac- 
cordingly, then  (-S»)  and  (-S.)  are  to  the  left  of  the  dominant  roots  (Fig. 
2- 19b).   If  S.  and  S~  are  chosen  on  the  boundary  of  the  dominant  root  re- 
gion (Point  C  in  Fig.  2-19E),  then  the  root  (-S-)  is  forced  to  be  in  a 
position  as  shown  in  Fig.  2-19C.   As  the  arbitrary  roots  move  toward  the 
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Fig.  2-19a   Cascade  filter. 
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Fig.  2-19b   Root  locus  of  Fig.  2-19(a)  in  case  Z>  p 
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Fig.  2-19C  Dominant  roots  on  the  boundary 
of  the  dominance  for  a  =  fc«JM 


Fig.  2-19D  Dominant  roots 
on  the  boundary 
of  stability 
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Fig.  2-19E  Stable  and  dominant  root-region. 
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Fig.  2-20  Dominant  root  region  (a  =  fU)n  )  and  Modified  Dominant 

Root  Region 
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origin  and  so  does  the  zero.   When  the  arbitrary  roots  reach  the 
boundary  of  stability  (Point  D  in  Fig.  2-19E),  (-S,,)  reaches  the  origin 
and  the  zero  and  pole  of  the  compensator  also  reach  the  origin  as  shown 
in  Fig.  2-19D.   From  this  analysis,  it  can  be  visualized  that  if  the 
dominant  roots  are  chosen  close  to  the  stability  boundary  and  within 
the  stable  region,  this  implies  a  dipole  close  to  the  origin.   Therefore, 
a  region  which  is  close  to  the  stability  boundary  and  between  the  boundaries 
of  stability  and  dominance  (which  are  defined  in  the  last  section  for  a  = 

f  CUn  )  is  also  a  dominant  root  region.   This  region  is  given  the  termin- 
ology "modified  dominant  root  region"  in  order  to  distinguish  it  from  the 
definition  of  the  last  section.  These  two  dominant  regions  are  shown  in 
Fig.  2-20. 

The  pole-zero  configuration  discussed  above  is  often  the  case  for 
systems  of  high  velocity  error  constant  (K  ).   For  very  high  velocity 
error  constant,  a  dipole  near  the  origin  is  necessary. 

The  analysis  above  is  based  on  the  tendency  of  the  root- locus  and 
the  known  zero.   The  same  analysis  can  be  applied  to  other  cases. 

In  this  case  the  dominant  roots  should  be  chosen  either  close  to  the 
stability  boundary  or  within  the  dominant  region  (a  =  £^n  ).   Choosing  the 
dominant  roots  in  the  vicinity  of  the  dominant  boundary  (a  =  V^f)  ) 
should  be  avoided  when  a  zero  is  close  to  the  origin  because  the  root  is 
not  close  enough  to  the  zero  and  it  may  cause  excessive  overshoot.   If  a 
zero  is  known  to  the  right  of  the  dominant  roots,  a  root  close  to  it  is 
necessary. 

Now  consider  the  pole-zero  configuration  of  Fig.  2-19D  which  is  the 
situation  of  stable  limit.   In  this  condition,  the  pole  and  zero  of  the 
compensator  are  at  the  origin,  therefore  the  presumed  dominant  roots  -S 
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and  -S_  are  on  the  original  uncompensated  root-locus  (K  as  variable). 
This  provides  another  means  to  evaluate  the  stable  boundary,  namely  the 
stable  boundary  is  the  uncompensated  root-locus  itself  in  this  case  (lag 
compensator).   This  boundary  line  is  one  of  the  dividing  lines  of  the 
root  relocation  zone  discussed  in  reference  7.   For  high  order  systems 
there  are  several  branches  of  the  stable  boundary  line,  but  for  the  pur- 
pose of  dominant  region  only  the  one  which  is  closest  to  the  jtO   axis  is 
of  interest.   From  the  root -locus  it  can  be  visualized  easily  that  the 
dominant  roots  are  within  the  region  which  is  the  closest  to  the  jUJaxis. 
From  this  analysis,  it  concludes  that  for  lag  compensator  the  region 
closest  to  the  nearest  branch  of  the  uncompensated  root-locus  is  the 
dominant  root  region.   And  the  dominant  boundary  of  the  definition  a  =  fU)n 
provides  a  guide  of  how  close  they  should  be. 

Example  2-4.   In  Fig.  2-21,  K  is  not  to  be  reduced.   Find  the  dom- 
inant and  modified  dominant  root  region: 
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MsN-as-H.*) 


Fig.  2-21  Example  2-4. 

The  expressions  for  the  dependent  variables  Z,  p,  C  and  C.  have 

been  derived  in  (1-22)  and  (1-23).   Substitute  Sl   +   S2  =  £JCJn  ,  S  S 

(jJL      ,  B  =  6,  B,  =  1.2,  B0  =  2  into  (1-22)  and  (1-23),  one  obtains 
n     o      1        £■ 
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Q  =  [(l-^JcjJ-1.2  +4-J^]^/(^<-^j>uJni-^J  (2-40-1) 

p   —     UJ*C0/ (o  (2-40-2) 

C,  =   a  +  p  -afoon  (2-40-3) 

CX    =    [^>nC*   +  £?^«Q  -  U?]/b  (2-40-4) 
From  Table   2-4: 


For  a   =   yuJn 


Po=    Ce-CA+CL* 


The   Routh's   criterion   for  D,    and  D      is  D,     >      0  and  D    ■>  0. 

1      o     1  o 

The  stability  boundary  line  (root-locus  of  the  uncompensated  system) 
is  shown  in  Fig.  2-22a.   In  order  to  evaluate  the  a  =  P  UJn  dominant 
boundary  line,  C  and  C  must  be  calculated.   For  ^  =  0.3  and  €  =  0.5,  the 
results  of  computation  from  equations  (2-40)  are  tabulated  in  Table  2-5a 
and  Table  2-5b  respectively.   The  dominant  boundary  line  for  a  = f"Jn   is 

Table  2-5a.   Example  2-4  for  J5  =  0.3 

u)n        Co  C,  f>                          o<  E  j>    J>0 

0.4      0.65  1.777  0.01735  0.0698  0.0249  +    + 

0.5      0.473  1.72  0.0197  0.0914  0.216  +    + 

0.6      0.278  1.657  0.01662  0.118  0.141  +    + 

0.7       0.0539  1.584  0.0044  0.1322  0.0335  + 

0.8  -0.085 
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Table  2-5b  Example  2-4  for  ^  =  0.5 


U)n 

Co 

c, 

P 

<X 

Z 

J>, 

£o 

0.2 

0.81 

1.8 

0.0054 

0.0543 

0.0995 

+ 

+ 

0.3 

0.615 

1.7 

0.00923 

0.0545 

0  1695 

+ 

+ 

0.4 

0.418 

1.611 

0.01115 

0.0685 

0.1625 

+ 

+ 

0.5 

0.213 

1.509 

0.00886 

0.079 

0.112 

+ 

- 

0.6 

0 

Shown  in  pig.  2-22a.   If  the  dominant  roots  ?  =  0.5S  cj  =  0.55  are  chosen 
which  is  close  to  the  stability  boundarys  the  closed  loop  pole-configura- 
tion is  shown  in  Fig.  2-22b.   For  this  choice  of  the  dominant  roots,  the 
compensator  is 

c  S  +  o.oo54-3 


68 


Fig. 

2-22a 

Dominant   root   region  of  Example   2-4  #            V 

Shaded  area   is  the  modified  dominant              * 
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Fig.  2-22b  Closed  loop  pole-zero  configuration 
of  Example  2-4. 
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CHAPTER  III 
FUNCTION  OF  VARIABLE  PARAMETERS 
3-1  Realization  Problem: 

The  result  of  the  partition  process  employed  in  this  method  is  to 
separate  the  dependent  variables  into  two  parts;   the  coefficient  of  the 
reduced  characteristic  equation  and  the  variable  parameters  (parameters 
of  the  compensator).   The  design  is  carried  out  essentially  by  two  in- 
dependent procedures.   One  is  to  choose  appropriate  closed  loop  pole- 
zero  configuration  from  the  reduced  characteristic  equation  and  the  other 
is  to  realize  the  compensator  according  the  chosen  pole-zero  configuration, 
Consider  a  cascade  compensator  of  the  following  transfer  function. 


C  sN   b,5  •+  b0) 


If  there  are  no  restrictions  to  the  type  of  compensator,  the  parameters 

a, ,  a  ,  b,  and  b  are  linearly  dependent  in  the  root-coefficient  relations, 
1   o   1      o 

Their  solutions  from  the  partition  functions  are  always  real.   For  the 
function  of  (3-1)  with  all  real  coefficients  always  can  be  realized  either 
by  a  passive  network  or  an  active  network.   However9  if  the  poles  of  the 
compensator  are  required  to  be  on  the  negative  real  axiss  then  (3-1)  be- 


comes 


c"         (s-rp.Ks.+  p^) 


(3-2) 


From  (3-2)  and  (3-1) 

?^?^-=  b,  ,  Br,  =  b0  (3-3) 
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In  equation  (3-3),  P  and  P_  are  quadratic  functions.   For  certain  choices 
of  the  roots  the  solution  of  P  and  P.  may  be  complex  numbers,,  although 

fc>  and  Jq0    are  always  real.   Therefore  the  choice  of  the  arbitrary  roots 
are  not  only  necessary  in  the  dominant  root  region  but  also  necessary  in 
the  region  which  gives  a  physically  realizable  compensator.   For  certain 
kinds  of  specifications  and  restrictions^  the  compensator  may  be  forced 

o 

to  be  an  active  network  ,  but  for  most  of  the  systems,,  the  designer  should 
try  a  passive  network  first.   The  same  realization  problem  may  arise  from 
feedback  compensators  and  even  pure  derivative  feedback. 

3-2  Linear  Dependence. 

If  the  variable  parameters  are  linearly  dependent  in  the  coefficients 
of  the  controlled  characteristic  equation,,  then  the  compensators  are 
physically  realizable  passive  or  active  for  every  choice  of  the  arbitrary 
roots,  because  solutions  of  them  are  always  real  numbers.   For  this  reasons 
the  variable  parameters  in  the  coefficients  of  the  controlled  characteris- 
tic equation  should  be  arranged  (as  possible)  such  that  they  are  linearly 
dependent.   In  many  cases,  the  linear  dependent  relations  can  be  obtained 
by  transformations.   Consider  a  single  section  cascade  compensator 


Q    =  — (3-4) 

If  fC  ,  E.     and  p  are  variable  parameters  then  the  term  |(z  is  a 
quadratic  term.   However,  by  defining    |<  £  —  <X  then  (3=4)  becomes 

Ci.  —       c.  ^  -       where  K"c  =  <* /k   (3~5> 

For  the  transfer  function  of  (3-5)  alone  the  parameters  are  linearly  de- 
pendent.  For  the  same  reason  a  single  section  compensator  in  the  feedback 
path  the  relations  can  be  linearized  in  the  same  way. 

The  linear  dependence  of  the  variable  parameter  is  a  sufficient  condi- 
tion for  the  compensator  to  be  realizable,,  but  it  is  not  necessary.   The 
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necessary  and  sufficient  condition  can  be  found.   For  example >  a  single 
section  cascade  compensator;,  the  roots  within  the  single  section  relocation 
zones  formed  by  the  rules  of  reference  (7)  are  the  necessary  and  sufficient 
condition. 

Realization  problems  may  even  arise  for  pure  derivative  feedback. 
Consider  Fig.  3-1  in  which  the  point  "X"  is  assumed  to  be  able  to  insert 
a  summer  and  the  signal  at  which  can  be  measured.   Assumed  pure  deriva- 
tive feedback  is  devised  as  shown  in  the  Fig.  3-1.   From  rules  of  signal 
flow  graph  the  determinate  of  the  system  contains  a  term  which  is  the  pro- 
duct of  the  two  nontouching  loop  transmissions.   This  term  contains  the 
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Fig.  3-1.   Pure  Derivative  Feedback 
product  of  h  and  h..   For  a  certain  choice  of  the  roots. 


h_  may  be 


complex  numbers.   For  some  purposes  in  control  systems ,  the  system  may  be 
expected  to  be  overdamped  and  to  operate  in  a  given  set  of  roots  such  as 
discontinuous  operation  in  reference  9,   For  that  purpose  the  scheme  of 
Fig.  3-1  should  be  avoided.   For  pure  derivative  feedback  to  achieve  a 
given  set  of  roots  certain  rules  and  methods  have  been  worked  out.   They  are 
discussed  in  detail  in  Appendix  I. 
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3-3  Choice  of  Arbitrary  Roots. 

Any  root  may  be  chosen  as  the  arbitrary  root.   For  convenience 
the  dominant  roots  are  usually  chosen  as  the  arbitrary  roots.   If  there 
is  only  one  arbitrary  root  the  damping  ratio   P  is  the  convenient  in- 
dependent variable.   If  the  arbitrary  roots  are  more  than  two,  two  of 
them  are  chosen  as  the  dominant  roots,  others  may  be  either  fixed  at  a 
certain  value  or  be  defined  to  have  a  certain  relation  with  the  dominant 
roots.   Suppose  the  original  system  has  a  zero  which  is  to  the  right  of 
the  expected  dominant  roots.   Then  in  the  case  of  three  arbitrary  roots 
one  of  them  may  be  chosen  to  cancel  that  zero  or  be  close  to  it.   In 
general,  the  more  the  arbitrary  roots,  the  more  the  design  freedom  one 
has.   The  forward  gain  K  may  be  always  treated  as  a  variable  parameter, 
so  one  additional  arbitrary  root  may  be  introduced  even  though  K  is  us- 
usally  specified  within  some  limits.   By  treating  K  as  a  variable  para- 
meter, sometimes  it  may  simplify  the  relations,  but  sometimes  may  make 
the  problem  more  complicated.   It  depends  on  the  specifications  and  the 
nature  of  the  compensator,  but  the  increase  of  design  freedom  is  always 
true. 

For  the  case  when  a  large  number  of  arbitrary  roots  i.e.,  a  large 
number  of  variable  parameters  is  involved,  a  clear  picture  of  the  closed 
loop  pole-zero  configuration  must  be  assumed.   With  that  pole-zero  configura- 
tion the  arbitrary  roots  can  be  properly  related  to  the  dominant  roots.   In 
the  mapping  process  from  the  contour  of  the  arbitrary  roots  to  the  vari- 
able parameters,  only  one  independent  variable  is  varied,  i.e.,  the  para- 
metrical  line  is  obtained  in  the  parametrical  space  (or  coefficient  space 
of  the  reduced  characteristic  equation).   There  are  infinitely  many  para- 
metrical lines  in  the  parametrical  space,  but  only  those  which  are  within 
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the  specification  limits  are  of  interest. 

3-4  Coefficient  as  Variables:   Mitrovic's  Method. 

In  case  of  coefficients  as  variables,  this  method  provides  a  relatively 
simple  relation.   For  lower  order  systems  and  for  derivative  feedback  the 
application  is  direct.   For  the  case  of  any  two  coefficients  taken  as 
variables,  the  function  of  the  variable  coefficients  are  the  same  as  Mitrovic's 
coefficient  functior  which  has  a  certain  interpretation  from  the  mapping  pro- 
cess.  In  order  to  show  the  coefficient  functions  are  the  same  as  that  of 

Mitrovic's  method,  take  a  5th  order  equation  with  B  and  B,  as  variables 

o  1 

from  Table   1-4,    and  make   substitution  of  S^+S^   =    £?  UJn         and  S-S?   =  U)n 
from  which  one  obtains: 

C,  —   B3-  2j^n  [\-  ^h)  -   W* 

=  B3  -  ^  B^n  +  (  <+f-  ij  u>* 


fo=  Hi  Sz-2f%U)n+  (V-liW'  -  ^(«"0^«  ]  (3-6) 


=  ^j  ba  +  ( i-tj*;  eX  +  ij  (<+f-0  <V°n  +■  c<tf -iK*  (3.7) 


Equations  (3-6)  and  (3-7)  are  the  same  expressions  as  obtained  from  Mitrovic's 

2 
method  .   By  the  same  manipulation,  functions  from  Table  1-1  to  Table  1-3 
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3 
are  the  same  expressions  as  obtained  from  the  modified  Nitric" s  method  . 

Mitrovic's  original  coefficient  function  is  derived  by  mapping  a  constant 

P   contour  through  the  characteristic  equation.   By  the  same  mapping 

process,  other  contours  also  can  be  used.   In  this  analytical  method  the 

functions  are  expressed  in  terms  of  complex  variables,  the  functions  for 

any  contour  can  be  obtained  readily.   Because  of  this  extension  of  Mitrovics 

Method  to  the  contour  other  than  constant    P    line,  some  of  the  features 

are  discussed  briefly  in  the  next  three  sections. 

3-5  Cauchy's  Theorem  Applied  to  the  Function  of  Variable  Coefficients. 

In  the  case  of  two  variable  coefficients  of  a  characteristic  equation, 
the  interpretation  of  the  parametrical  line  of  f^  vs  j-  -   is  discussed 
in  detail  in  reference  (2)  and  reference  (4).   Here  only  the  cases  of 
variable  coefficients  which  are  not  in  sequence  are  discussed.   Take  the 
case  of  -f-,  and  f3      as  variable  coefficients.   Assume  the  original  character- 
istic equation  is 

F  =  5% +  63W  BaSa+  6,S  +  B0       (3-8) 

and  assume  the  controlled  characteristic  equation   is 

f  =    5n4    ----  -t  f35*+   &/+  f.S-f  B0  (3-9) 

If    -f,     and     f  3     are   so  chosen  such  that   -f    -  0S    then 

F   =  (8i-h)s*+   (6,-fJS  (3-10) 

Assume  the  arbitrary  roots  are  on  a  constant   P  line  in  the  left  half  of 
the  complex  plane,  then  the  angle  of  the  S  term  of  (3-10)  is  /*_        and 
the  angle  of  the  S  term  is  /&   if  B_  >  f"3    and  B,  >  ff     respectively. 
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If  B^  <  J ^  and  B1  ■<.  -f ,   ,  then  the  angle  of  the  corresponding  term 

is  izi@      an<i  /-Q    respectively  as  shown  in  Fig.  3-2 (b).   The  functions 
of  -f,  and   f  3   have  negative  values  as  cJ<<  I  from  Table  1-3.  Assume  all 
the  roots  of  equation  (3-9)  are  to  the  left  of  the  f     line8  then  the  mapping 
of  U)n   from   0  -*  Co  line,  must  give  counterclockwise  loops.   At  the  4  points 
on  ptS)     curve  as  shown  by  (1),  (2)  (3)  and  (4),  f,    and  -f      must  have 
the  following  values: 

at  point  (1)     f3  —  B3      f,  <  B, 

at  point  (2)     f3>B3     i ,=  B,  (3-11) 

at  point  (3)     f3z=  B»3     f,  ^  B, 

at  point  (4)     f3  <  ^     -f(  _  B, 

By  tracing  the  curve  -£   vs  -f  by  (3-11) ,  the  f  vs  f  curve  must  have 
a  counterclockwise  loop  and  enclose  the  point  M.   In  general^  the  curve  f. 
vs  f.  in  which  j  >^+  1,  then  a  counterclockwise  loop  of  the  curve  f\  vs  -f  ■ 
implies  all  the  roots  are  to  the  left  of  the   f   line.  \ :  point  M  is 
inside  the  loop.   In  addition  to  the  conclusions  from  the  reference  (2)  and 
reference  (4),  the  interpretations  of  -f.   vs  f-     curve  are  summarized  as 
follows : 

(1)  If  J  ■  i  +  1  (  j  =  0S  1,  2J,.,.),  clockwise  loop  of  f,  vs  f- 
curve  implies  all  the  roots  are  to  the  left  of  the  ?  line  providing  point 
M  is  inside  the  loop. 

(2)  If  j  >  i+1  (j=0,  1,  2,...),  counterclockwise  loop  of  -f^ 
vs  f«  curve  implies  all  the  roots  are  to  the  left  of  the  f*  line  pro- 
viding point  M  is  inside  the  loop. 

The  above  rules  can  also  be  applied  when  the  mapping  contour  is  a 
constant  vertical  line  (  cx_  constant). 
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3-6  Evaluation  of  roots  from  the  f".  vs  f:  plot. 

The  TX  vsfjplot  is  a  mapping  process  from  the  complex  plane.   The 
mapping  functions  of  this  process  are  the  functions  of  the  variable  co- 
efficients, namely  -f .    =   -+-.  (s).   In  this  mapping  process  regions  map 
into  region,  curves  map  into  curves.   An  unbounded  region  on  the  S-plane 
maps  an  unbounded  region  unto  the  coefficient  plane,  while  bounded  regions 
map  unto  bounded  regions.   Therefore,  the  stable  root  regions  discussed  in 
the  last  Chapter  map  a  stable  region  in  the  coefficient  plane.   The  map- 
ping process  from  S-plane  to  coefficient  plane  is  unique ,  but  the  inverse 
(from  coefficient  plane  to  S-plane)  mapping  is  not  necessarily  unique,  it 
may  have  multi-values.   By  using  this  mult i- value  property  of  inverse  map- 
ping, more  than  one  pair  of  complex  roots  can  be  evaluated  from  the  coef- 
ficient plot.   Figure  3-3  is  a  4th  order  equation  with  J~,   and  "F?  as  vari- 
ables.  Fig.  3-3(a)  shows  the  stable  root-region  and  Fig.  3-3(b)  shows  the 

~f,  vs.  f  plot.   It  can  be  seen  from  Fig.  3-3  that  the  unbounded  stable 
root  region  on  S-plane  map  an  unbounded  region  unto  f  'v  f  plane. 

The   f  =0.5  lines  on  S-plane  map  a  clockwise  loop  unto  the  plane. 
Every  point  on  S-plane  map  a  point  on  f^f.,  plane.   However,  the  inverse 
mapping  is  double-valued,  such  as  points  (2)  and  (3)  on  f   ^  j-       plane 
map  two  points  (2)'  and  two  points  (3)"  on  S-plane.   Therefore,  if  suf- 
ficient curves  have  been  drawn  on  f-  ^  Tj  plane,  by  using  different  con- 
tours (such  as  constant  $    and  constant  c\.  lines),  two  pairs  of  complex- 
roots  can  be  evaluated  from  the  T,  ^  Tj   plane  directly.   For  the  case  of  a 
4th  order  equation,  all  the  roots  can  be  seen  on  the  'ff'^f-   plane  for 
every  choice  of  f;  and   "f~  •   .   Real  roots  are  not  shown  on  the  coeffici- 
ent plot,  but  they  can  be  evaluated  easily  from  the  reduced  characteristic 

equation. 
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Since  for  every  region  on  the  S-plane,  there  is  a  corresponding 

region  on  the  coefficient  plane,  the  required  sys;-m  root-region  can  be 

mapped  on  the  coefficient  plane  if  such  a  regv'a  exists,   For  example  in 

Fig.  3-3,  if  the  roots  are  required  such  t*»t  all  roots  have   J—  0,5 

and  time  constant  less  than  1,  then  tb.«  r°°t  region  of  the  equation  has 

a  corresponding  region  on  the         Plot  which  is  the  intersection  re- 

gion  of  the  clockwise  loops  p ,*M  =  1  and  )    =   0.5„   This  region  is  shown 

by  the  shaded  area.   For  ^er   requirements,  the  region  on  the  coefficient- 

plane  which  satisfied  the  requirement  can  be  located  by  the  same  procedure. 

/ 

Figure  3-4  &A   Figure  3-5  are  numerical  examples  with  variables  f", 
z 
// 
and  t  for  4t¥^rder  an<*  5th  order  equations  respectively.   In  Figure  3-4 
/ 

only  one  constant  C\-  (=  1)  contour  is  drawn  and  in  Figure  3~5S  only  one 

If' 
constant^    (  =  0.5)  is  drawn.   From  those  figures,  it  can  be  seen  that  the 

constant     and  constant  CK.    contour  constitute  curvilinear  coordinates 

on  the  t"  '  -  plane.   From  the  point  of  view  of  curvilinear  coordinates, 

two  pairs  of  complex- roots  can  be  read  for  every  point  on  the  f-  plane. 

3-7  Stability  boundary  on  the  coefficient-plane  when  the  variable  coefficients 
are  not  in  sequence: 

As  discussed  in  Section  Three  of  Chapter  Two,  when  the  variable 
coefficients  are  not  in  sequence,  the  functions  of  the  variable  coefficients 
are  undefined  on  the  ( ^U  axis.   Consequently,  the  stability  of  the  bound- 
ary line  on  the  coefficient-plane  cannot  be  obtained  by  mapping  the  ^axis 
of  the  S-plane.   However,  the  boundary  on  the  S-plane  is  well  defined  as 
has  been  shown  in  Chapter  Two,  therefore,  the  boundary  on  the  coefficient- 
plane  can  be  obtained  by  mapping  those  boundary  lines  from  the  S-plane  to 
the  coefficient  plane.   This  process  can  be  carried  out  with  point  by  p oint 
calculations  or  by  approximating  the  boundary  on  the  S-plane  by  a  function. 
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The  boundary  in  Fig.  2-12  has  been  calculated  point  by  point  and  plotted 
in  Fig.  3-6.   In  Fig.  3-6,  a  constant  P  =0.1,  line  is  also  mapped  onto 
the  ~f  /v  f     plane.   The  P  *  0.1  line  shows  a  counterclockwise  loop  which 
implies  that  there  exists  a  region  in  which  all  roots  have)>0,l. 
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PART  II  DESIGN 
CHAPTER  IV  DESIGN  PROCEDURE 
4-1   Specifications. 

The  transient  response  of  a  control  system  is  determined  completely 
if  the  closed-loop  pole-zero  (root-zero)  configuration  is  known.   The 
contributions  to  the  transient  response  by  each  pole  and  zero  can  be 
calculated  precisely.   If  a  pair  of  complex  roots  is  dominant s    then 
approximations  can  be  made.   Based  on  these  dominant  roots2  the  effect  of 
all  the  zeros  and  all  the  secondary  roots  to  the  transient  response  can 
be  visualized  from  the  complete  root-zero  configuration.  s   With  this 
relation  the  specifications  from  the  desired  response  can  be  conveyed  to 
root-zero  configuration.   As  in  this  method,  the  design  is  based  on  root- 
zero  configuration,  all  specifications  must  be  interpreted  in  a  set  of 
appropriate  configurations  of  roots  and  zeros.   Frequency  response  is 
readily  interpreted  if  dominant  roots  are  assumed.   The  steady  state 
specifications  (stiffness  and  error  constant)  are  contained  implicitly  in 
the  equations  and  formulas  to  work  with,  then  they  are  satisfied 
cally  if  the  final  results  satisfy  all  other  specificati 
meet  all  the  specifications  may  not  be  possible  until  some  adjustments  are 
made.   For  example,  to  design  a  system  of  high  velocity  cons« 
dipole  near  the  origin  may  be  required.   If  the  original 
are  defined  so  that  all  secondary  roots  are  far  away  to  the  left  of 
dominant  roots,  it  may  result  that  no  root-zero  configurations  for  this 
definition  of  dominant  roots  are  able  to  meet  the  high  K  requirement. 
In  this  case,  the  definitions  of  dominant  roots  may  necessarily  be  adjust- 
ed so  that  a  dipole  to  the  right  of  the  dominant  roots  is  allowed. 
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In  other  words,  a  term  of  long  tail  of  relatively  small  magnitude  in 
the  transient  response  is  allowed.   Usually,  this  long  tail  term  in 
transient  response  is  avoided.   But  in  order  to  achieve  high  velocity 
error  constant  by  a  relatively  simple  network,  this  adjustment  may  be 
acceptable  in  consideration  of  all  the  factors  in  design. 

4-2  Design  Procedure  --  General. 

In  the  design  of  control  systems,  many  methods  and  tools  are  avail- 
able such  as  the  root -locus  method,  frequency  domain  method,  pole  zero 
configurations,  etc.   Problems  frequently  rise  about  what  procedures  and 
patterns  are  more  logical.   In  pure  synthesis,  the  designer  is  able  to 
take  the  specifications  and  proceeds  logically  in  a  step  by  step  manner 
to  a  system  which  will  meet  these  specifications.   In  these  senses  many  of 
the  design  techniques  are  not  true  synthesis  procedures.   On  the  other  hand, 
the  design  procedure  proposed  by  Guillemin  which  takes  the  form  analogue 
to  modern  filter  theory  is  a  synthesis  technique.   In  Guillemin' s  procedures 
the  closed  loop  transfer  function  is  first  determined  f rom t he  specifications 
and  the  compensation  is  determined  subsequently.   The  design  procedure  by 
the  analytic  method  developed  in  this  paper  is  more  or  less  close  to  Guillemin's 
procedure.   Briefly,  the  procedures  are  as  follows: 

(1)  The  regions  of  closed  loop  pole-zero  configurations  are  deter- 
mined from  the  specifications. 

(2)  The  type  of  compensator  (cascade,  feedback  pure  derivative, 
or  mixed  network)  is  chosen  from  the  consideration  of  specifications  and 
other  restrictions. 

(3)  Analytic  expressions  of  all  the  roots  and  the  parameters  of 
compensation  in  terms  of  the  arbitrary  roots  are  obtained. 

(4)  Within  the  pole-zero  region  defined  in  (1),  choose  a  set  of 
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appropriate  pole-zero  configurations.   Then  the  parameters  of  the  compen- 
sator are  calculated  from  the  expressions  in  (3). 

Procedure  (1)  is  an  interpretation  or  transition  from  the  specifica- 
tions to  the  closed  loop  pole-zero  configurations.   Procedure  (2)  depends 
upon  the  specifications,  the  nature  of  the  plant  and  also  on  engineering 
judgment  and  experience.   These  two  procedures  are  general  to  all  methods, 
if  the  design  is  based  on  transient  response  or  time  domain.   They  are  the 
prerequisites.   Procedures  (3)  and  (4)  are  concerned  with  the  method  it- 
self.  For  this  analytic  method,  they  are  explained  in  detail  in  the 
following  sections. 

4-3  Design  Procedure  --  Detail. 

When  the  type  of  compensation  has  been  chosen,  the  design  procedure 
for  this  methid  is  as  follows: 

(1)  Formulation  of  the  Characteristic  Equation. 

The  characteristic  equation  with  the  compensator  is  formu- 
lated in  the  way  such  that  the  coefficients  of  this  equation  are  functions 
of  the  free  parameters  of  the  compensator.   In  other  words,  all  parameters 
are  independent  from  one  another. 

(2)  Formulation  of  the  Root-Coefficient  Relations. 

Set  the  coefficients  in  (1)  equal  to  the  corresponding  root 
coefficient.   If  the  order  of  the  equation  is  "n"  there  are  "n"  relations. 
If  there  are  "r"  free  parameters  in  the  compensator,  then  "r"  roots  are 
arbitrary  and  the  reduced  characteristic  equation  is  of  the  order  (n- r) 
and  has  (n-r)  coefficients.   The  root-coefficients  relations  are  formulat- 
ed in  such  a  way  that  the  coefficients  are  functions  of  the  arbitrary  roots 
and  the  coefficients  of  the  reduced  characteristic  equation. 
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(3)  Formulation  of  the  Functions  of  the  Dependent  Variables. 

In  the  "n"  root  coefficients  relations  formulated  in  Step  2, 
there  are  "n"  dependent  variables  and  r  independent  variables.   The  'V 
dependent  variables  are  the  "r"  parameters  of  the  compensator  and  the  (n-r) 
coefficients  of  the  reduced  characteristic  equation.   The  independent  vari- 
ables are  the  "r"  arbitrary  chosen  roots.   Since  there  are  "n"  equations  and 
"n"  dependent  variables,  the  dependent  variables  can  be  solved  and  expressed 
in  functions  of  the  arbitrary  roots. 

(4)  Determination  of  Stability  and  Dominance. 

From  the  reduced  characteristic  equation,  determine  the 
stability  root-region  and  the  dominant  root-region.   The  dominant  root 
region  depends  upon  the  degree  and  the  definition  of  the  dominance  which 
is  defined  by  the  designer  from  the  consideration  of  specifications. 

(5)  Evaluation  of  the  Compensator  Parameters. 

If  the  desired  dominant  roots  are  within  the  dominant  root- 
region  which  has  been  determined  in  step  (4),  then  choose  the  dominant 
roots  (usually  independent  variables)  and  substitute  them  into  the  formula 
of  step  (3),  then  the  parameters  of  the  compensator  are  evaluated. 

(6)  Realization  of  the  Compensator  Network. 

If  the  parameters  of  the  compensator  evaluated  in  the  last 
step  are  within  the  specified  range  and  can  be  realized,  then  the  design  has 
been  completed. 

(7)  Design  Adjustment. 

If  the  value  of  the  compensator  parameters  are  not  in  the 
specified  range  or  the  desired  dominant  roots  are  not  in  the  dominant -root 
region,  then  it  indicates  either  this  type  of  compensation  cannot  meet  the 
requirements  or  the  definition  of  the  dominant  roots  must  be  adjusted. 

89 


other  types  of  compensation  have  been  chosen  or  the  degree  of  dominance 
has  been  adjusted,  the  design  procedure  is  just  a  repeat  of  the  above 
steps. 

The  procedures  described  above  are  general.  The  first  three  steps 
are  the  partitioning  process  as  discussed  in  Chapter  I.   For  some  specific 
cases,  formula  have  been  derived  as  shown  in  the  subsequent  sections,  then 
the  design  is  just  a  computation  and  is  started  from  Step  (4).   The  compu- 
tation is  best  carried  out  by  tabulating  all  previous  data  and  by  pro- 
gressing computation.   Once  those  data  are  tabulated,  all  information 
can  be  found  from  the  table  and  if  some  changes  of  the  type  of  the  com- 
pensation are  made,  most  of  the  data  can  be  utilized. 

4-4  Free  Parameters  of  the  Compensator. 

In  step  (2)  of  the  last  section,  a  terminology  "Free  Parameters"  was 
used.   It  may  need  further  explanation.   Consider  a  cascade  single  section 
R-C  compensator.   Let  the  transfer  function  be 

r  Kc(S  +  l) 

C<c-  -JTp- 

The  apparent  parameters  are  K  ,  Z  and  b     .   Assume  the  system  is  type  I 
and  the  velocity  error  constant  K  is  specified.   If  the  forward  gain  K 
of  the  uncompensated  system  was  determined  by  the  K  requirement  and  in 
the  partitioning  processes,  K  was  treated  as  a  given  numerical  constant 

then  K  is  not  arbitrary  in  order  to  meet  the  requirement  of  K  .   In  other 

c  J  v 

words  K  must  satisfy  the  following  relation: 

Kc=  P/a 

Thus  the  three  apparent  compensator  parameters  K  ,  fc>  ,  and  Z  are  related 
by  the  requirement  of  K  .   Only  two  of  them  are  independent  variables,  the 
third  one  is  dependent.   The  independent  variables  of  the  compensator  in 
this  sense  are  given  the  terminology  "free  parameters". 
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For  other  cases,  the  same  treatment  can  be  applied  to  the  compen- 
sator parameters.   In  the  above  example,  if  the  forward  gain  K  was  also 
treated  as  an  adjustable  parameter,  then  all  the  three  apparent  para- 
meters are  independent  because  K  K  is  also  independent.   However,  the  re- 
sult of  the  partitioning  processes  are  different  for  the  two  cases.   But 
the  design  can  be  carried  out  by  either  case.   In  the  former  case  (K  is 
constant)  there  are  two  arbitrary  roots,  while  in  the  latter  case  (K  is  ad- 
justable) there  are  three  arbitrary  roots.   The  expressions  of  the  para- 
meters are  different  for  different  cases,  but  the  design  can  be  carried 
out  in  both  cases. 


91 


CHAPTER  V  DESIGN  OF  PURE  DERIVATIVE  FEEDBACK 
5-1  Tachometer  Feedback  -  Second  Order. 

For  compensating  lightly  damped  instrument  servos ,  tachometer  feed- 
back  is  commonly  used.   For  second  order  cases,  as  all  information  can 
be  determined  precisely,  only  the  formulas  are  presented. 


Fig.  5-1  Second  Order  System  with  Tachometer 

Feedback 

The  block  diagram  of  a  Type  1  second  order  system  is  shown  in  Fig. 

5-1.   The  separation  of  the  forward  gain  by  K.  and  K~  is  for  the  purpose 

of  proper  adjustment  of  h.   Assume  the  specifications  are  given  K, K«  and 

\   .   The  design  procedures  are  as  follows: 

(1)  Formulation  of  the  Characteristic  Equation 
The  Characteristic  Equation  is 

5^+  (cl4  Kih)S  4  K,K,  =  0  (5-1) 

(2)  Formulation  of  the  Root-Coefficient  Relation. 

Since  there  is  only  one  parameter  of  compensations  there  is  only  one 
independent  variable.   The  independent  variable  is  not  necessarily  the 
arbitrary  root  itself,  it  can  be  a  function  of  the  roots.   In  the  case 
of  one  arbitrary  root  and  that  root  has  to  be  complex,  the  independent 
variable  can  be  chosen  as  either   f   or  UOn    of  the  pair  of  complex 
conjugate  roots.   For  second  order  case,  the  only  choice  is   ?   , 
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because  GOn  *-s  fixed  by  the  gain.   By  doing  this  the  root-coefficient 
relations  are: 

a  +-  Kth   =■  I  f^>n  C5=2-i) 

K,K,=  0Jn  (5-2-2) 

(3)   Formulation  of  the  Functions  of  the  Dependent  Variables, 
Since  the  root  of  the  reduced  characteristic  equation  is  the  domin- 
ant root  itself,  there  are  no  stability  and  dominant  region  problems. 
Here  the  dependent  variables  are  cOn  and  h(  ?  is  chosen  as  independent 
variable).   Their  expressions  as  functions  of  the  independent  variable  are 
as  follows: 

From  (5-2-2): 

From  (5-2-1):     j<dh  =  Z$U)„  -  O.  =  af/K^-CL 

^fTMC-^-  (5-3-2) 

For  the  given  value  of   P  ,  h  is  computed. 

5-2  Third  Order:   Tachometer  Feedback. 

The  block  diagram  of  a  third  order  type  one  system  with  tachometer 
feedback  all  around  the  plant  is  shown  in  Fig.  5-2. 
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Fig.  5-2  Third  Order  System  with  Tachometer 
Feedback 
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The  characteristic  equation  of  the  uncompensated  system  is 

53+  dzS2-t   B,S  -f  £0  =0  (5-3) 


where  B2  =  a  +  b,   B.  «  ab,    B  =  K.IC. 

Assume  the  forward  gain  K  and  the  damping  ratio  )   based  on  the 
dominant  roots  are  given.   The  design  procedures  are; 

(1)  Formulation  of  the  characteristic  equations   It  is  easily 
formulated  as  follows: 

i>3+   B2s'  +  (  B,  +  KA)S  -T-    B0~0  C5~4> 

(2)  Formulation  of  the  root-coefficient  relation:   Since 
there  is  only  one  control  parameter  h,  only  one  root  is  arbitrary. 

This  root  can  be  chosen  either  the  real  root  or  ^  of  the  complex  roots. 
Assume  the  real  root  is  chosen  and  designated  as  -S.  (S.  is  positive  on 


the  negative  real  axis).   Then  S.  is  the  independent  variable. 

duced  characteristic  equation  is  of  order  2  and  therefore  has  the  form 

5Z  +  C,S   -f  C0  =  O 

Since  the  roots  of  this  equation  are  actually  the  dominant  roots9  it  is 
expressed  by  P  and   U-)n  for  convenience.   Here  ?     and  U)    are  not  in- 
dependent variables  but  dependent  variables.   By  the  choice  of  the 
above  definition,  the  root-coefficients  relations  are: 

B^  =  Z/UJ^  t-S,  (5=5=1) 

B0=  ^nS>  C5-5  =  3) 

(3)  Determination  of  the  functions  of  the  dependent  variables. 
Here  the  three  dependent  variables  are  P      ,  u)    and  h.   In  the 
actual  computation,  only  one  dependent  variable  must  be  found  explicitly 
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as  a  function  of  the  independent  variable,  the  others  can  be  computed 
by  using  the  implicitly  functions  of  formulas  (5-5).   But  for  the  pur° 
pose  of  illustration,  they  are  solved  and  shown  in  equations  (5-6) 


5  ={B,-S.,)Js,/ZJb-0  (5-6-1) 

UJh  =  J  B./S,  (5-6-2) 

K2h  =  B0/5,  +  6,S,  -5,'-  S,  <5"6-3> 

Note  the  three  variables  f  ,  C<Jn  and  h  are  expressed  as  functions  of  the 
independent  variable. 

(4)  Determination  of  stability  and  dominant  root  region: 

In  this  one  control  parameter  case,  the  stability  and  dominant  root 
regions  are  lines.   These  lines  are  the  same  as  the  root  locuss  by  treat- 
ing h  as  variable  and  K  as  a  fixed  number.   The  only  difference  between 
them  is  the  independent  variable.   In  the  root- locus  method9  the  indepen- 
dent variable  is  h  while  in  this  method  it  is  the  arbitrary  root  S.o   A 
numerical  example  is  shown  in  Fig.  5-4.   The  degree  of  dominance  can  be 
visualized  easily  from  those  plots. 

For  this  particular  scheme,  it  indicates  a  max.  damping.   This  max. 
damping  can  be  found  from  equation  (5-6-1)  analytically  because  f*  is  ex- 
pressed explicitly  as  a  function  of  S. .   Differentiate  €   with  respect  to 
S.  and  set  it  to  zero,  one  obtains: 

Substitute  into  equation  (5-6-1),  obtain 

?*   =  &/z.lb0  (5-8) 

Substitute  this  value  of  S.  into  ?  and  C<Jn  ,  the  real  part  of  the 
complex  root  becomes 
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f^l  =  ®J5    =  s, 


m 


(5-9) 


Equation  (4-7)  shows  that  if  the  tachometer  loop  gain  K~h  is  so 


adjusted  such  that  the  real  root  is  (   — 


3 


)  then  the  complex 


roots  have  max.  damping  for  the  given  value  of  B  (  =  K)  and  this  max. 

o 

damping  is  shown  in  equation  (5-8).   Equation  (5-9)  shows  that  for  max. 
damping  the  real  part  of  the  complex  roots  is  equal  to  the  real  root. 
Fig.  5-3  shows  this  root  configuration. 


Fig.  5-3  Root  Configuration  for  max.  f   of  3rd  Order 
System  with  Tachometer  Feedback. 

For  this  root  configuration,  the  tachometer  loop  gain  and  oJ  are 

obtained  from  equations  (5-6-2)  and  (5-6-3)  as  follows; 


Kth  =  -^  +K-8, 


(5-10) 


(5-11) 

Equation  (5-8)  gives  a  quick  look  of  the  max.  possible  damping  for  given 
K  or  the  max.  K  for  a  given  9    in  a  third  order  system  with  tachometer 
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Fig.  5-4(a)   Plot  of  Equations  (5-6-1)  and  (5-6-2)  for  numerical  values 

given  in  example'  5-1. 


Fig.  5-4(b)   Root  locus  of  the  characteristic 
Equation 
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feedback. 

Example  5-1:   For  the  block  diagram  shown  in  Fig,  5-2,  if  b  a  109 
a  =  5,  K^o  =  750;  the  uncompensated  system  has  a  pair  of  complex  roots 
on  i CO  axis.   For  this  pole-zero  configuration:   B0  =  15,  B,  =  50,  B  =  750. 

^  JL  1  O 

From  equation  (5-8),  the  max.  possible  damping 


IS' 


xiso 


From  (5-7)   S.   =  5. 
1m 


Assume  this  value  of  f    is  taken,  then  from  (5-6-2)  and  (5=6-3) 

<^n=  12.3 
K2h  =  18.5 
For  other  root  configurations,  ?   and  C0n  as  a  function  of  S.  are  plotted 
in  Fig.  (5-4) (a)  and  the  root-locus  for  h  as  variable  is  shown  in  Fig.  (5=4) 
(b).   If  the  dominant  root  is  defined  so  that  the  secondary  root  has  greater 
real  part  than  that  of  the  dominant  roots,  then  the  portion  of  the  root 
locus  when  S.  ^  5  is  the  dominant  root  region. 

5-3  Third  Order:   Tachometer  Feedback  to  Interior  Node. 

Consider  Fig.  5-5  in  which  point  "X"  is  assumed  to  be  available  to 
insert  a  summer.   The  analytic  design  procedures  are: 


O 
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Fig.  5-5  Tachometer  Feedback:   Minor  Loop 
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(1)  The   characteristic  equation   is: 

53+   {  62+  KA)S*  +   (B,+  bKzh)S  +  &0—0  (5-12) 

where       B0   =  a+b,        B.    =  ab,      B     =  K.K0 
2  '1  o  1   2 

(2)  Root-coefficient   relationship.      For  the   same   definition 
of  dependent  and   independent  variables  as  described   in  the   last   section^ 
the  relations  are 

6z+    K2h  =    ZJUJ^  +  S,  (5-13-1) 

B,  +  bk,h  =  U)U   2;«0„S,  (5-13-2) 

B0  =  U);  5,  <5-"-3> 

(3)  Function  of  dependent  variables  (  ?  ,  aJ  s   K,  h   )• 
From  equation  (5-13-1)  and  (5-13-2),  eliminate  K~h 

Substitute  C»Jn  from  (5-13-3)  into  the  above  expression  obtain 

e   —  B'b~6'  +"^r  ~b5'  (5-i4) 

The  other  two  variables  can  be  obtained  in  the  same  way  but  are  not 
necessary. 

(4)  Stability  and  dominant  root  region:   For  the  numerical 
values  given  in  Example  5-1,  the  dominant  roots  P  and  oU  are  plotted 
in  Fig.  5-6  as  functions  of  S.  from  equation  (5-14).   The  degree  of  dom- 
inance can  be  seen  easily. 

(5)  Take  f  =  0.5  which  corresponding  to  S  =  20.   From  (5-13- 

3) 

#  CO0=  6.14 

From  either    (5-13-1)   or    (5-13-2),    obtain 

Kth  -  11.14 
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Fig.    5-6       %       and     CJn  as  a   function  of  the  real 
roDt   of  the   systems   shown   in  Fig.    5-5. 
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5-4  First  and  Second  Derivative  Feedback. 

In  case  of  first  and  second  derivative  feedback  around  the  plant 
for  a  given  value  of  forward  gain,  only  the  coefficients  of  the  first 
and  second  derivative  term  are  variable.   The  design  may  be  carried  out 
using  two  variables.   But  for  more  freedom  it  is  better  to  take  the  con- 
stant term  also  as  a  variable  then  three  arbitrary  roots  are  available. 
The  coefficients  of  the  reduced  characteristic  equation  and  functions  of 
the  variable  coefficients  have  been  derived  and  tabulated  in  Table  1-6. 

Example  5-2.   Given  the  original  characteristic  equation, 

S*+  >0S*  +  3IS*  +  30S   +  <  =  O  (5-15) 

The  forward  gain  K  must  be  greater  than  1000.   The  gain  at  the  stability 

limit  for  the  uncompensated  system  is  84,  first  and  second  derivative  feed- 
back are  chosen.   When  three  variable  coefficients  are  used,  equation  (5- 
15)  becomes 

S%  |OS%  f^-f-f.5  +    f0=O  (5=16) 

where  -f  =  K 

From  Table  1-6,  the  coefficient  of  the  reduced  characteristic  equa- 
tion and  the  functions  of  the  variable  coefficients  are; 

C0=  BJ-(-Sl  +  S2  +  S3j  (5  =  17-1) 

f0  =  C0S,S>ZS3  (5=17-2) 

f,    =    CCJA-tSA+SjSj  (5-17-3) 

fa=  C0(S1  +  S1+53jt(i^+Vi+^(J      <5"17-4> 
where  (-S..),  (-S2)  and  (-S.)  are  the  arbitrary  roots  and  C  =  S,.   Here 
the  order  of  the  equation  is  four,  all  the  roots  may  be  complex.   There 
is  only  one  constrained  root  so  it  is  convenient  to  replace  C  by  S,  it- 
self.  The  order  of  the  reduced  characteristic  equation  is  one;  equation 
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(5-17-1)  alone  determines  all  the  roots.  Assume  (-S-)  and  (_S2)  are  the 
dominant  roots,  then  (5-17-1)  becomes 

10  —   Z$U)n   +  S,  ti+  (5-18-1) 

fo=^S1S4  (5-18-2) 

There  are  no  zeros  in  this  system,  all  secondary  roots  (-S-  and  -S.)  are 

required  to  be  to  the  left  of  the  dominant  roots.   Assume  (-S„)  and  (-S,) 

are  real  roots,  then  for  the  definition  of  dominance,  the  max.  K  occurs 

when  S-  ~z  S,  =  f^Or\       •   From  (5-18-1)  for  this  value  of  the  secondary 
roots,  the  conditions  are: 


fu)n=  1.5 


(5-19-1) 


*  .  ± 


mQX.    K  =  f^n  (5-19-2) 

For  P  =0.5,  the  max.  K  =  156.  This  is  far  below  the  required  K,  then 
the  secondary  roots  must  be  complex.  Assume  P  and  u^  denote  the  second 
roots,  then  equations  (5-18)  become 

fiO   4  f '(*>„'  =  5  (5-20-1) 

k  =  u>:^a  (5-2o-2) 

Notice  in  equation  (5-20-1)  there  are  three  arbitrary  variables  and  this 
equation  alone  determines  all  the  roots.   Assume  the  dominance  is  defined 
as    <?  UJn  >■  P^n  then  from  (5-20-1),  the  dominant  root  region  is 

Pu3  <■  2.5,  which  is  shown  in  Fig.  5-7a.   Within  this  region  choose 
C  =  0.5,    -.Cr   =  2,  then  froir  (5-20-1)  the  condition  for  secondary 
roots  is  P  <4,  =  4 

Substitute  into  (5-20-2),  obtain     «  =.  fo^"/  J/i 
For  K  =  1000,  obtain   ^  =  0.253.   For  this  choice  of  the  dominant  roots 
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f  and  f~  are  readily  calculated  from  (5-17-3)  and  (5-17-4)  as  follows 


fx   =  532 
f2  =  270 
The  root  configuration  is  shown  in  Fig.  5-7b.   For  another  choice  of  the 
dominant  roots,  the  computation  can  be  carried  out  in  just  the  same  way. 
Example  5-3:   Given  the  original  characteristic  equation 
5s"  4  i<?JS%  U8S*+  £zos*i  I5oo5  -i-  K   —  O 
The  forward  gain  K  must  be  greater  than  2000.   For  this  value  of  K  the 
system  is  badly  damped,  first  and  second  derivative  feedback  are  chosen. 
By  using 
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Fig.  5-7a 
Dominant  Root  Region  of 
Example  5-2 


Fig.  5-7b 
Root  -configuration  of  Example  5-2 
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three  variable  coefficients,  the  working  formulas  from  Table  1-6  are  as 
follows : 

C,=  %~(  S^S^S5)  (5-21-1. 

C0  =:  B3  -  C,  (  3,  +  Sa  4-%]  _  (  ^  +  JA  +  5,  S,  )  (5-21-2) 

fo  -  S,^Co  (5~21"3) 

f,  =  Co(5,-S,  +  ^sj  +  -SjSj  +  C,  5,5^3 


(5-21-4) 
(5-21-5) 


Assume  (-S.)  and  (-S„)  are  the  dominant  roots.  The  order  of  the  equation  is 
five,  one  root  must  be  real.  Let  this  real  root  be  (-S~).  Here  three  roots 
are  arbitrary.   S~  may  assume  any  value.   For  dominance,  assume 

5j  =  aj^oOn  (5-22) 

Substitute  (5-22),  B,  =  16  and  B^  =  128  into  equation  (5-21),  one  obtains: 

C,  =   Id,-  4^n  (5-23-1) 

C0=  U8  -4-yu)nC,-(|  +  4fJu;na  (5-23-2) 

ff    =   C0UJ(;(l+«lJ  ^JCdnC,  (5-23-4) 

fx  =  4f(0nC„  +(l+  rj^nC,  ^fU^  (5-23-5) 

Notice  (5-23-1)  and  (5-23-2)  determine  all  the  roots,  (5-23-3)  determines 
the  forward  given  K  while  (5-23-4)  and  (5-23-5)  determine  the  derivative 
gain.   Since  there  are  no  zeros  in  the  system,  the  secondary  roots  are 
required  to  be  to  the  reft  of  the  dominant  roots.   Let  a  =  f  ^r\   »  the 
transformed  coefficients  from  Table  2-4  are  as  follows: 
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£,  =  C,~   ^J^n  (5-24-1) 

J>a  =   C0  ~;cJnC,  +  f*60*  (5-24-2) 

The  stability  boundary  lines  are  determined  by  C.>°and  Q    >  O  ° 

The  dominant  boundary  lines  are  determined  by  J)>o  and   I)  >  O 

i  © 

They  are  shown  in  Fig.  5-8a.   For  P  =0.5,  the  results  of  computation  of 

C  .  C,.  f  .  D  and  D,  are  tabulated  in  Table  5-1.   Examine  Table  5-1. 
o   1   o   o      1 

when  COn  is  between  4  and  5  the  gain  K  is  greater  than  the  specified 
value  2000.   This  range  of  C0n  is  within  the  dominant  region  since  both 
D.  and  D  are  positive.   The  root-configuration  for  (0  =4,    ^  =  0.5s, 
is  shown  in  Fig.  5-8b.   For  this  choice  of  the  dominant  rootsa  the  vari- 
able coefficients  are  computed  readily  from  (5-23-3),  (5-23-4)  and  (5-23-5) 

f  =  K  =  2040 
o 

f  =  1534  (original  1300) 
f2  =  576  (original  520) 
For  other  values  of  S-  and  other  choices  of  the  dominant  root  the 
computation  is  the  same. 


U)n 

c, 

Co 

fo 

J>, 

^o 

2 

12 

72 

576 

+ 

+ 

3 

10 

50 

1350 

+ 

+ 

4 

8 

32 

2040 

+ 

+ 

5 

6 

18 

2250 

+ 

+ 

6 

4 

8 

1728 

- 

+ 

7 

2 

2 

686 

+ 

8 

0 

0 

0 

Table  5-1  Example  5-3  for  €  =   0.5 
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Fig.  5-8a  Stability  and  Dominant  Root  Regions  of  Example  5-3 
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Fig.  5-8b  Root  Configuration  of  Example  5-3  for  R    =  0.5 
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CHAPTER  VI 
FEEDBACK  COMPENSATOR  DESIGN 
6-1   Introduction. 

The  compensation  by  pure  derivative  feedback  discussed  in  the  last 
chapter  usually  adjusts  only  some  of  the  coefficients  of  the  character- 
istic equation.   The  effectiveness  of  this  type  of  compensation  is  less 
than  that  if  the  compensator  parameters  enter  all  the  coefficients.   This 
has  been  shown  by  the  reduced  characteristic  equation  in  Chapter  II. 

Moreover,  the  tachometer  feedback  reduced  K  and  acceleration  feedback  re- 

v 

duced  K  and  so  on;  and  a  system  can  seldom  be  stabilized  without  adjust- 
ing  the  lower  order  coefficients  of  the  characteristic  equation.   There- 
fore, if  the  compensator  is  chosen  in  the  feedback  paths  the  compensator 
often  consists  of  a  tachometer  followed  by  a  network.   The  tachometer 
effectively  changes  the  coefficient  of  the  first  derivative  term,  while 
the  poles  and  zeros  of  the  network  enter  all  the  coefficients  in  the  way 
of  the  coherent  nature  of  algebraic  equations.   This  chapter  is  a  discus- 
sion about  the  design  of  such  a  compensator  by  the  analytic  method. 

6-2  Second  Order  Compensator. 

Consider  Fig.  6-1  which  is  a  second  order  system  with  feedback  compen= 
sator  consisting  of  a  tachometer  followed  by  a  lead  network. 


K, 
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5(-S  +  <\) 


KaS' 


s+  p 


Fig.  6-1  Second  Order  System  with  Feedback  Compensator, 
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The  characteristic  equation  of  the  compensated  system  is 

_5%(  e.+p+KA)^-*-  (  60+  B,p;s  4ft,jo  —  o   (6-i) 

where  B,  =  a,  B  =  K,  K„.   In  this  characteristic  equation,,  there  are 
1       o    1  2 

two  variables  in  the  coefficients,  therefore  two  roots  are  arbitrary 

if  K  and  P  have  no  restrictions.   The  reduced  characteristic  equation 

is  of  order  one  and  has  the  form  of  S  +  C  =0.   The  root  coefficient  re- 

o 


lations  are 


e>0p  =  cu)n 


(6-2-1) 

(6-2-2) 

6-2-3) 


In  equation  (6-2),  the  dependent  variables  are  to  ,  k  and  C  ;  the  in- 
dependent variables  are  ?     and   ^n  whi-ch  are  assumed  to  define  the 

dominant  roots.   C  determines  the  stability  and  dominant  root  region,, 

o 

while  p  and  jb  must  assume  values  according  to  the  dominant  roots. 
Since  equations  (6-2)  are  linear  functions  of  the  dependent  variables^ 
it  is  convenient  to  write  in  the  matrix  form  as  follows: 


/ 

0 
/ 


I 


6,   -2 

a.    - 


- 1 

K^Ro. 

f*>« 

p 



^ 

Co 

— 

0 


»-3) 


The  three  dependent  variables  can  be  solved  readily,  but  only  C  must  be 
expressed  as  an  explicit  function  of  the  independent  variables.   By  Cramer's 
rule,  one  obtains 


(6-4) 
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C   in  this  case  is  the  actual  real  root.   The  stability  and  dominant 
o 

root  region  are  determined  by  equation  (6-4)  by  varying  \    and  ^n    „ 

For  design  purposes,  only  those  values  of  ?  and   OO^  which  are  within 

the  specifications  are  computed.   When  ?  and  uJ    have  been  chosens 

K«ka  and  p  are  calculated  and  the  design  is  completed; 

Example  6-1:   In  the  block  diagram  Fig.  6-1,  K..K,,  =  10  ,  a  =  23„ 

Therefore  B  =  10  and  B,  =  23.   From  equation  (6-4): 
o  1  c 

The  stability  criterion  is     Cq^O        i.e, 


2.X|0^tJn-23u£ 

or 

These  two  inequalities  define  the  stability  region  and  is  shown  in  Fig, 
6-2.   If  the  dominant  region  is  defined  by     C  ^>  <2  f  ^On    then  the 
criterion  for  this  definition  is 

10  "^    ->a?tOn 


This  region  is  also  shown  in  Fig.  6-2.   For  <?  =  0.5,  the  range  of  Cdn   is 
0<C60n<220.   The  values  of  CQ,  p  and  K_ka  for  the  dominant  roots  in 
this  region  (  ?  =  0.5)  are  computed  as  shown  in  Table  6-i.   The  final 
decision  of  the  design  is  to  pick  up  a  value  which  is  favorable  to  all 
compensator  parameters  and  other  specifications.   From  Table  6-1  assume: 


^n 

Co 

F> 

M 

10 

io4 

10 

9770 

20 

5000 

20 

4960 

40 

2480 

38.4 

2460 

60 

1600 

56 

1580 

80 

1200 

77 

1180 

100 

922 

92 

900 

200 

312 

125 

365 

Table  6-1  Values  of  p  ,  f^ka  and  C 

of  Example  5-1  for  ?  =0.5 


109 


Fig.  6-2  Stable  root-region  and  Dominant  root-region 

of  Example  6-1 
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The  dominant  roots  are  chosen  to  be  ^  =  0.5,  (xj    -   1008  then  the  third 
real  root  is  -922,  p   =  92  and  K~ka  =  900.  The  closed  loop  transfer  fane- 

tl0n  ^  ©S  =     K.KA^P) 

The  closed  loop  pole-zero  configuration  of  the  above  choice  of  the  domin- 
ant roots  is  shown  in  Fig.  6-3. 


1** 

-X- 


*  fr 


-12. 


^tf.5 


JW 


+  <^v 


Fig.  6-3  Closed  Loop  Pole-zero  Configuration  of 
Example  6-1  for  the  dominant  roots 
<J   =0.5,   uJn  =  100 

From  Fig.  6-3,  it  can  be  seen  that  the  system  is  essentially  second 

order  because  the  real  root  is  far  away  from  the  origin.   The  effect  of 

P  to  the  transient  response  is  also  easily  visualized  from  the  pole  zero 

configuration. 

6-3  Third  Order. Compensator. 

Fig.  6-4  is  a  block  diagram  of  a  third  order  system  with  feedback 
compensator.   Assume  the  forward  gain  K.K„  is  fixed  from  the  steady  state 
accuracy  and  the  velocity  constant  must  be  greater  than  a  given  value  K  . 

Ill 


•^ 

K, 

+.o 

K2 

J 

- 

•v. 

— 

5(S+aJ(S  +  b; 

Fig.  6-4  Third  Order  System  with  Feedback  Compensator. 
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The  characteristic  equation  of  the  compensated  system  is 


.4 


where   B  —  0.+  b^    8(  =  <3b,  B0=K,K^ 

In  the  above  equation,  there  are  three  variabless  namely  fc     „  jo   and  £ 

However,  since  K  is  specified,  those  three  variables  are  not  independent 

one  to  another,  they  are  constrained  by  the  specification  K  .  K  of  the 

v  v 


compensated  system  is 
Rearrange  equation  (6-6) 


B,P  (6-6) 

PB,+  K^2 


Substitute  this  expression  into  the  coefficient  of  the  first  derivative 
term  of  equation  (6-5),  it  becomes 

B,p  t  K^ti  t  B0  =  6,p  +  ■£&  -PB,  *  Bo  =  J^?  +  B0 


Then  equation  (6-5)  becomes 


V' 


8. 


S4  +  (  62+  p)  S3+  (6,  t  BiP  +  IC^t)5i+  (T°P  t  8„)  S  -t  6pp  =  O    (6"7) 

The  coefficients  of  characteristic  equation  now  have  only  two  variables » 
namely  i>   and  p  .   The  specification  of  K  eliminates  one  of  the  vari= 
ables  of  the  compensator.   This  is  always  true.   The  more  specifications 
given,  the  less  the  freedom  of  the  compensator.   Since  there  are  two  vari- 
ables in  the  coefficients,  two  roots  are  arbitrary.   The  reduced  character- 
istic  equation  is  of  order  two  and  has  the  form: 
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B  *f>  =  2j^n-t-C 


[6-8-1) 


6,+  B,p+  |C2£t  =  ui*  4  Co   ■*  2JH.C 


6 


(6-8-2) 

(6-8-3) 

(6-8-4) 


where   §   and  CcL  are  the  arbitrary  roots  while  C  and  C,  are  the  co- 

*  n  J  o      1 

efficient  of  the  reduced  characteristic  equation.   As  the  four  dependent 
variables  C. ,  C  , p and  ft>.  in  the  above  equations  are  linear  functions8 
it  is  convenient  to  express  by  matrix  form  as  follows: 


-2fU)n 


B 


-w. 


-cdf 


r<v 


o 

"C|~ 

1 

a 

0 

P 

0 

Mt 

Bo 

0 


By  Cramer's  rule,   one  obtains 


c  = 


B  -  &2&0<*>*   +  2joUn 


M^n-fX2+  ^ 


(6-10) 


The  stability  and  dominant  region  are  determined  by  C  and  C, .   Choose 

o      1 

dominant  roots  the  parameters  K,Kt   »  b  anc*  Z.  °f  tne  compensator  are 
calculated  from  equations  (6-8)  and  (6-6) 

Example  6-2.   For  the  system  of  Fig.  6-4  given  a  =  10,  b  ■  30,  K„K?  = 

21,000  K  —  60,  damping  ratio  P      of  the  dominant  roots  be  approximately 
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0.7  and    U)     =  10. 

For  the  given  values,   B~  =  40 

B.  =  300 

B  =  21,000 
o 

Substitute  into  equation  (6-10)  for    €  =  0.7 

f  44X/08-  8.+  *IO*U)Z   4  i^^n 


For  (Jn  =  10,       Cq  =  1321 

From  (6-8-4)  p  =  6.26 

From  (6-8-1)  C,  =  \    +  p  -  2^n  =  32.3 

C.  and  C_  determine  the  secondary  roots,  they  must  be  checked  before 

going  further.   The  reduced  characteristic  equation  is 

5*+  32  2  S  +  1320  =0 

The  roots  of  this  equation  are  -16.1  +  j  32.5.   Compare  the  real  part  of 
these  roots  to  that  of  the  dominant  roots  (  ?LOiy   =7),  it  indicates  the 
presumed  dominant  roots  are  dominant. 
From  (6-8-2) 

Mt  =  ^n  +  C„  +  2j>  UJnCt   -  B,  -  6,p  =  /3^3 
and  finally  from  (6-6) 

The  compensator  becomes 


The  calculations  above  are  for  one  point  only.   The  choice  of 
dominant  roots  does  not  insure  a  satisfactory  result.   For  more  freedom 
and  adjustment,  the  stability  and  dominant  root  region  are  calculated 
first  as  shown  in  Fig.  6-5.   From  Fig.  6-5,  adjustments  of  the  dominant 
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Fig.  6-5.   Stability  and  Dominant  Root-regions  of  Example  6-2 


Nfc.i  +  i3^r 


t<7^ 


~(0 


Fig.  6-6.   Closed  loop  pole-zero  Configuration  of  Example  6-2  for 


3 


=  0.7,   O>0=  10. 
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roots  can  be  made. 

The  closed  loop  transfer  function  is  readily  formulated  as 

0R  "~  cs+s.jcs+s^cs-fs^cs+s^ 

The  closed  loop  pole-zero  configuration  for  ?      -   0.7a    uj  =  10 
(dominant  roots)  is  shown  in  Fig.  (6-6). 

Example:  6-3.   Given  the  same  plant  and  the  compensator  as  example 
6-2  except  the  specifications  of  K  is  relaxed.   Here  are  three  free  ad- 
justable parameters,  three  arbitrary  roots  are  available.   The  control- 
led characteristic  equation  is  the  same  as  (6-5) ,  the  root-coefficient 
relation  are  as  follows: 

Batp    =   2?oJn+±2  +  S±  (6-12-1) 

B,+  %p+K^t  =  U>»+ifU„(^^)  +  %s*  (6-12-2) 

B,p  +  K2*t  £  t  K  =  H(V5J  +  ^f^A       <6"12-3> 

The  four  dependent  variables  are  K«K  ,  p  ,  jz  and  S,9  the  three  independent 
variables  are   C   ,   <-oL  and  S~.   Assume  S  and  S,  are  complex^  and  are 
denoted  by  j  and  CJn  .   If  (^jy  is  chosen  as  the  dependent  variables  then 
from  (6-12-1)  and  (6-12-4)  one  obtains 

This  equation  alone  determines  all  the  roots  since  it  contains  P   ,  cd      s 

P       ,  and  (jj     .  Assume  f  (dominant  root)  =  0.59  ^    (  secondary  root) 
=  0.7,  and  substitute  the  numerical  values  of  K  =  21s0008  B  -  300 
B~  =  40  into  (6-13)  and  solve  for  U)n   ,  one  obtains 

'    W1XIO*   f     , . " 
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i-14) 


Assume  only  10   :fr  CiJ   —  20  are  of  interest,  the  results  of  compu- 
tations from  Equations  (6-12)  and  (6-14)  are  tabulated  in  Table  6-2. 

Notice  in  Table  6-3  when    U)n   =  10  and  Lj     =  23.5  Z  is  negative. 

270      348  62660         10.25         34 

23.5      2.63  575  -7.25        -3.62 


10 


15 


20 


109.2     128.5         9085  17.2  7.48 

21.4 

54.4      56.4  5324  22.3  2.53 

19.6 

Table  6-2  Example  6-3  for  f=   0.5,   ?'*  0.7 
For  simple  realization  of  the  compensator,  the  roots  of  cJn  =  15, 
C/J  =  109.2  are  chosen.   For  this  choice  of  the  roots,  the  compensator 

is  CT  =     $tS(s  +  l1±L 

where    $   =     L  ■-  •   If  K  =  1,  k  =  0.433.  fie  four  roots  of 

t  Kz 

the  system  are  f  =  0.5,  u>  =  15,  j>'  =  0.7,  U)^-  109.2.  There 
is  a  zero  at  -128.5.  From  the  closed  loop  pole-zero  configuration,  the 
roots    ^  =0.5,   UJn  =  15  have  high  degree  of  dominance. 
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CHAPTER  VII 
CASCADE  COMPENSATION  DESIGN 
7-1  Single  Section: 

As  cascade  compensation  is  commonly  used,  the  design  procedure  dis- 
cussed in  Chapter  IV,  is  applied.  In  Chapter  I  a  third  order  system  was 
analyzed,  here  a  nth  order  system  in  the  general  case  is  considered. 


a 


i  _ 


k(S+2) 


U  +  R) 


K 


Sn+/Lls"-+- 


+  4. 


Fig.  7-1  Single  Section  Cascade  Compensation  for  a 

nth  Order  System 

Consider  the  system  of  Fig.  7-1.   Assume  the  characteristic  equa- 
tion of  the  uncompensated  system  is 

.2 


S%Bn_,5"-'+B„.2S"-% 


+  as  +  6.s  +  &—o    o-i) 


For  specified  forward  gain  K  and  error  constant  (K  ,  K  ,  etc.)  the 
three  compensator  parameters  Z,    Jp  and  <X  are  not  independent  one 
to  another,  but  related,  namely 

:  -£—  (7-2) 


The  characteristic  equation  of  the  compensated  system  is 


sTihws% 


+  (8^4  B^pJ-S^  --+(B«*$f>)s        (7-3) 
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In  equation  (7-3),  there  are  two  variables  t?  and  o{    in  the  coef- 
ficients, then  two  roots  are  arbitrary.   The  reduced  characteristic 
equation  is  of  order  (n-1),  and  has  the  following  form. 

5""'+   CmS""**    ---     +c,s  +  C„=0 

If  the  two  arbitrary  roots  are  presumed  to  be  the  dominant  roots  of  the 
compensated  system  and  are  expressed  by    P   and    CO      then  the  root- 
coefficient  relations  of  equation  (7-3)  are: 


6op  =-  C0u)n 
B,   +  E^p  =  C2^n    +  C,ifU„  4  Co 


B, 


(7-4) 


Be.,  t    &*P  =C-u)H%   C,zW  +    C2 


The  N+l  dependent  variables  (    p    ,   o(     ,  C  ,  C.....C  2)  are  linear 
functions,  then  equation  (7-4)  can  be  expressed  in  matrix  form 


3n  0 

B„-.  0 

B„-,   0 


0 

B,  3C 
B„  o 


O 

o 


o 


0 

o 
- 1 


o 
-I 

-2fK 
o 


0 


-  I 


KJi. 


-2f 


o 


-  I 


-2(H, 


-w 


Co 

c. 


c 


r 


0 

o 


_^n-2 


0 


ih 


B, 

O 


(7-5) 
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By  Cramer's   rule  one  obtains 

-B0[^B,  +  ?&<*>«  + V&Ldt*  4      ---  *?&(*)„"  (7.6) 

where  7,    —  -  I 

(7-7) 

■ 

I 
t 

The  other  dependent  variables  can  be  expressed  in  the  same  way  but  this 
is  not  necessary  as  far  as  computational  purposes  are  concerned.   When  C 
has  been  calculated,  the  other  variables  are  readily  evaluated  from  equa- 
tions (7-4).   Choose   f   ,  evaluate  the  stability  and  dominant  root 
interval  of  cJn  .   Within  the  dominant  root  region,  choose  u)n  and  evalu- 
ate p  and  0{      .   If  the  network  is  realizable,  the  design  wrok  is 
completed. 

Example  7-1.   Given   Q  =  — r£  and  K  is  not  to  be  reduced.   Then 
the  characteristic  equation  of  the  uncompensated  system  is 

S   +  400  =  0 

B„  =  1,    B,  =  0,   B  =  400,   BJ  =  0  where  i  :>  2 
2         1        o  i 

The  order  of  the  compensated  system  is  three  and  that  of  the  reduced 

characteristic  equation  is  one.   Then 

C,  ■  1.   C.  =  0   when  i  Zf    1 
1         1 
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From  Equation  7-4: 


B,+  B^p  =  C.ijto,,  +  Co 
Substitute  the  values  of  B's  into  the  above  equations,  obtains 

400  p  =  Co^o* 

p  —    I^UJ0  +Co 
From  equation    (7-6)   and    (7-7) 


(7-10) 


Substitute  the  values  of  B's,  obtain 


c    _       6co  f  o^»  (7.U) 

For  stability  the  criterion  is   C0--^  0»  namely 

■ — -H —  >o 

(*j£  -400 
The  numerator  never  becomes  negative,  then  the  inequality  becomes 

Define  the  dominant  roots  such  that  the  real  part  of  the  secondary  roots 

must  be  at  least  twice  that  of  the  dominant  roots,  then 

_  goo  ?<*>» 

C°  >  2  J40"  '     CJi-^O      >  2^" 
Then  the  criterion  for  dominant  roots  is 

Those  regions  are  shown  in  Fig.  7-2. 
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Within  the  dominant  root  region,  for   f  =0.5,  and   P  =  0.69  the 
compensator  parameters  are  computed  from  equation  (7-8)  for  different 
values  of  u)      and  tabulated  in  Table  7-1.   Table  7-1  includes  not  only 
the  compensator  parameters,  but  also  all  the  roots  of  the  system  for  dif- 
ferent choice  of  the  compensator.  Assume  the  dominant  roots  are  chosen 

f  =  0.6,  CO  =  30,  then  from  Table  7-1,  C  =  28.8,  p   =  64.8, 

p  r        Wi(S+if.8) 

n(      =4.41,        £  =  •  =  14.8,  The  compensator  is  0"  =  77TT 

^  cX,  c    (546f:8J 


as 


f  —   0.5T 

f-'ot(o 

^n 

Co 

P 

<* 

Co 

P 

<*v 

22 

220 

242 

13.2 

264 

291 

14 

24 

54.6 

78.5 

4.72 

67.7 

94.3 

5.34 

26 

37.7 

63.4 

4.12 

45.3 

76.1 

4.6 

28 

29.2 

57.2 

4 

35.1 

68.5 

4.41 

30 

24 

54 

4.05 

28.8 

64.8 

4.41 

32 

20.5 

54.6 

4.2 

24.6 

65.6 

4.53 

34 

19.3 

55.7 

4.53 

23.2 

66.8 

4.87 

Table  7-1  Computation  Result  of  Example  7-1 
The  closed  loop-pole- zero  configuration  for  this  choice  of  the  dominant 
roots  is  shown  in  Fig.  7-3(a).   In  Fig.  7-3(a)  it  can  be  visualized  that 

the  peaking  time  is  small  but  the  overshoot  may  be  too  large  because  the 

c. 
zero  is  too  close  to  the  origin  .   If  P  =  0,6,  <^n  =  24,  are  chosen  as 

the  dominant  roots,  then  from  Table  7-1,  C  =  67.7,  f>    =   94.3  t  <U    ■  5.34 

and  £  =17.7.   The  compensator  is 

__   SMC    S-HMj^ 

The  pole-zero  configuration  for  this  choice  of  the  dominant  roots  is  shown 
in  Fig.  7-3 (b).  For  this  choice  of  the  dominant  roots  the  peaking  time  is 
still  small,  but  the  overshoot  is  reduced. 
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$tcJ»tHy 


-30        -zo       -10 
Fig.   7-2 

Stable   and  Dominant  Region  of  Example  7-1 


*<x 


-t<^ 


Fig.   7-3(a) 

Closed  Loop  Pole-zero  Configuration  for  Dominar*- 

Roots     P  =  0.6,  iJn  =  30   for  Example   7-1. 


Fig.  7-3(b) 


*<K 


Closed  Loop  Pole-zero  Configuration  for  Dominant 
Roots  €   =  0.6,  uj '  =  24  of  Example  7-1. 
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Example  7-2:   Given  (r    =  ~77Z ^   TT    and  K  is  not  to  be 

reduced.   The  characteristic  equation  of  the  uncompensated  system  is 

S3  +  l£   +  10^  +  300  J  0 
From  Equation  7-4: 

7-f-  p  =  c,  +  2^oj0 

lO-fTP=    <^m   +c"   +^5'UJ"C.  (7~12) 

(Op  +  3ooot  =  C,^  +  zfoJ„Co  t      loop  =  C0^„ 
From  Equation    (7-6)   and    (7-7) 

For    P    =  0.5 

°~      76J^-CJ^-3O0 
For  =0.3 

To  determine  the  stability  region,  the  criterion  is 


Co  ^o  and         C,>0  (7"14) 

Define  the  dominant  roots  by  a  =  fU)n   ,  from  Table  2-7  of  Chapter  II, 
the  transformed  coefficients  of  the  reduced  characteristic  equation  are: 

P,  =  c/-^Q_ 

The  criterion  for  dominant  roots  by  this  definition  of  dominant  roots  is: 

P,  ^>o  and  J>0^°  <7"15) 

Computations  for  P  =0.3  and  f   =  0.5  from  Equations  (7-13)  and  (7-12) 
are  shown  in  Table  7-2  and  Table  7-3  respectively. 
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0Jn 

Co 

P 

Cj 

Oi 

£, 

"*0 

0.2 

9.13 

0.001217 

7 

0.0046 

+ 

+ 

0.4 

8.22 

0.0044 

6.88 

0.0101 

+ 

+ 

0.6 

7.25 

0.0084 

6.81 

0.0166 

+ 

+ 

0.8 

6.34 

0.0135 

6.49 

0.0235 

+ 

+ 

1.0 

5.28 

0.0176 

6.41 

0.0313 

+ 

+ 

1.2 

4.16 

0.02 

6.27 

0.0397 

+ 

+ 

1.4 

2.97 

0.0193 

6.21 

0.049 

+ 

+ 

1.6 

1.71 

0.0146 

6.07 

0.057 

+ 

+ 

1.8 

0.397 

0.0043 

8.86 

0.097 

+ 

+ 

1.9 

-0.312 

-0.0038 

4.73 

+ 

- 

Table  7-2   P  =  0.3  for  Example  7-2 


u)n 

C0 

P 

c, 

<X 

», 

Po 

0.6 

5.84 

0.007 

6.4 

0.019 

+ 

+ 

0.8 

4.46 

0.0095 

6.2 

0.0191 

+ 

+ 

1.0 

3.07 

0.01022 

6.01 

0.0216 

+ 

+ 

1.1 

2.24 

0.00904 

6.02 

0.0300 

+ 

- 

1.2 

1.56 

0.00748 

6.05 

0.0322 

+ 

- 

1.3 

0.87 

0.0049 

5.75 

0.036 

+ 

- 

1.4 

0.1925 

0.00125 

5.61 

0.0375 

+ 

- 

1.5 

-0.48 

-0.0036 

5.46 

0.0383 

+ 

- 

Table  7-3  £>=  0.5  for  Example  7-2 
From  the  criteria  (7-14)  and  (7-15),  the  stability  and  dominant  region 
for  P  =  0.3  and  f  =  0.5  are  determined  from  the  computation  tables. 

Stability  Dominant  (  a~  f^J 

f=0.3  UJn<l.8  LJn<l.8 

^=0.5  u)„<\.4r  uJn<\.oS 

Since  this  computation  is  a  progressing  processs  when  the  stability 
and  dominant  root  region  have  been  located,  the  compensator  parameters 
also  have  been  calculated  already.   The  computation  tables  not  only  in- 
clude the  compensator  parameters,  but  also  the  secondary  roots.   The 
coefficients  (D's)  of  the  transformed  equation  are  shown  signs  only  in 
the  compensation  tables,  because  in  this  case  only  signs  are  needed  to 
determing  the  dominant  region. 
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Assume  the  dominant  roots  are  chosen  as  <?   =  0.3,  LUn   =  1.0,  then 

from  Table  7-2,  C  =  5.28,  C  =  6.41,  p    =  0.0176,  OK    =  0.0313  and 

p 
/  —  — —-    =   0.563.   The  reduced  characteristic  equation  is 

*-    o<, 

S2  +  6.41S  +  5.28  =  0 
Then  the  secondary  roots  are  S  ■  -5.44,  -0.97.   The  closed  loop  pole- 
zero  configuration  is  shown  in  Fig.  7-4(a).   In  Fig.  7-4(a)s  all  the 
secondary  roots  are  to  the  left  of  the  dominant  roots. 

Consider  Table  7-2  and  7-3,  they  indicate  C{  ■<-  1,  this  implies  the 
compensator  is  a  lag  network.   For  a  lag  network  the  zero  is  close  to 
the  origin.   Then  a  root  to  the  right  of  the  dominant  roots  is  allowed. 
This  dipole  near  the  origin  has  negligible  effect  on  the  transient  response 
if  the  root  is  close  enough  to  the  zero.   In  this  sense,  the  definition  of 
the  dominant  roots  may  be  adjusted  as  such  that  all  the  secondary  roots 
must  be  to  the  left  of  the  dominant  roots  except  the  root  which  is  close 
to  the  zero. 

In  this  example:   if  $  -   0.5,  U)n   =  1.2  is  chosen  as  the  dominant 

roots,  from  Table  7-3,  Cq  =  1.56,   (^  =  7.25  |b  =  0.00748,  <*  =  0.035 

p 
and  i   —  — -  =  0.214.   The  reduced  characteristic  equation  is 
Ov 

S2  +  7.25S  +  1.56  =  0 
Then  the  secondary  roots  are  -7.03  and  -0.225.   The  closed  loop  pole- 
zero  configuration  is  shown  in  Fig.  7-4(b).   Compare  the  pole-zero  con- 
figurations of  Fig.  7-4(a)  and  Fig.  7-4(b),  it  can  be  easily  seen  that 
the  scheme  of  Fig.  7-4(b)  has  better  performance  than  that  of  scheme 
Fig.  7-4(a). 
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Fig.  7-4 (a)  Closed  Loop  Pole-zero  Configuration  for  the  Dominant 

Roots    f    =  0.3,  u)n  ~   1  of 
Example  7-2 
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Fig.  7-4(b)   Closed  Loop  Pole-zero  Configuration  fr>r  a  Dipole 
Near  the  Origin  of  Example  7-2.   Dominant  Roots 
are   P  =  0.5, 


LJn  =  1.2. 


7-2  Stability  Boundary  Lines. 

It  has  been  pointed  out  that  the  stability  boundary  lines  for  a 
cascade  compensator  are  the  same  boundary  lines  of  relocation  zones  in 
Reference  7.   In  the  single  section  zone  if  the  compensator  is  to  be 
realized  by  an  R-C  circuit  with  pole-zero  on  negative  real  axis,  then 
the  necessary  and  sufficient  condition  in  all  the  roots  are  within  that 

zone.   That  is  in  the  single  section  zone  the  compensator^  Q    ~ 

c   U  +  pj 
always  can  be  realized  by  an  R-C  circuit  for  any  choice  of  the  roots  pro- 
viding those  roots  are  within  that  zone.   In  the  multiple  zone9  the  roots 
being  in  that  zone  is  a  necessary  condition  but  not  a  sufficient  condi- 
tion. 

For  higher  order  systems  there  may  be  several  zones  for  a  single 
section  compensator.   But  for  dominance,  only  the  zone  which  is  the 
closest  to  the  j6U  axis  is  of  interest.   It  has  been  proved  that  putting 
a  pair  of  complex  roots  within  the  lag  area  and  close  to  the  stability 
boundary  lines  (uncompensated  root- loci)  implies  a  closed  loop  dipole 
close  to  the  origin.   By  the  same  reason,  putting  a  pair  of  complex  roots 
within  the  lead  area  and  close  to  the  boundary  line  (pseudo-root  loci) 
implies  all  secondary  roots  are  far  away  from  the  dominant  roots.   The 
question  of  how  far  away  and  how  close  they  are,   the  dominant  root  bound- 
ary lines  of  the  definition  a  =  6"  0Jn       serves  a  guide  for  both  condi- 
tions. 

7-3  Double  Section  Cascade  Compensator. 

If  the  dominant  roots  are  expected  to  be  within  the  multiple  sec- 
tion zone,  at  least  two  sections  are  necessary.   In  general,  the  design 
is  carried  out  by  two  steps.   First  choose  an  intermediate  dominant  root 
which  is  close  to  the  stability  boundary  lines.   Form  the  transfer  function 
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obtained  by  the  intermediate  dominant  roots,  then  repeat  the  process. 
The  working  formulas  for  the  open  loop  transfer  function  with  one  zero 
in  it  have  been  derived  in  Appendix  2  for  an  original  nth  order  system. 
The  design  by  intermediate  dominant  roots  has  two  advantages.   First  the 
quadratic  term  as  discussed  in  Chapter  3,  equations  (3-3)  is  avoided. 
Secondly  the  compensator  can  be  designed  as  a  double  lead  sections,  a 
double  lag  section  or  a  lead-lag  network.   The  intermediate  dominant 
roots  should  be  chosen  close  to  the  final  dominant  roots  because  the 
regions  which  are  close  to  the  dividing  lines  (both  uncompensated  root- 
loci  and  the  pseudo  root-loci)  are  dominant  root  regions. 

4-0  o 
Example  7-3:   Given    Cx   ==  ~    ,  the  system  is  expected  to  operate 

on  a  pair  of  dominant  roots  €    =0.7,   00^  =  20,  and  K  is  not  to  be  re- 
duced. 

From  Example  7-1  and  Fig.  7-2,  U)^   =  20  is  on  the  stability  boundary 
line,  then  double  section  is  needed.   From  Table  7-1,  choose  the  inter- 
mediate dominant  roots  P  =  0.5,  u)    =  22,  which  is  close  to  the  final 
dominant  roots  P  =0.7,  U)   =  20.   The  transfer  function  becomes 

5Z30(jS  +  \8S)  (7-16) 

Sx(  S  +  *4aJ 

From  Appendix  2  equation  A-2-4,  the  four  equations  are 
p-C,  =  Zf^h  ~  Bi 

B,p+  |a*  -Co  -^HC,  =  ^-Z,-*-^ 
60p  -  w,tc0  =  o 
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(7-17-1) 

(7-17-2) 
(7-17-3) 
(7-17-4) 


From  A-2-6: 

s~     ^ — .  — — (7-ioj 

where  B  ,  B,  and  B0  are  the  coefficients  of  the  characteristic  equation 
o   1      2 

of  the  open  loop  transfer  function  (7-16),  (-C)  is  the  open  loop  zero. 
The  stability  boundary  lines  are  shown  in  Fig.  7-5a  by  the  rules  of 

reference  7.   From  Fig.  7-5,  it  indicates  that  the  roots  of   f  =  0.7, 

LOy.    =  20,  are  within  the  single  lead  section  zone  and  close  to  the  pseudo 
root-loci.   Then  they  are  expected  to  be  dominant.   Substitute  C  =  185s  B 

97500,  B.  =  5280,  B2  =  242,  P   =0.7,  (J„  =  20  into  equation  (7-18) ,  one 
obtains 


From  (7-17-4) 
From  (7-17-1) 
From  (7-17-3) 
Then 

The  secondary  roots  are  determined  from  the  equation 


c 

o 

=  6340 

P  = 

26 

V 

=  240 

(X  = 

1.4 

z  = 

18.6 

(7-19) 


The  solution  of  (7-19)  is  -210,  -30. 

The  compensator  for  the  second  stage  is    Cr  —  - 


The  overall  compensator  is 


i8.s  a  +  fz.&f 


H:  (s  +  24-2J  Cs+^&J 

The  closed  loop  pole-zero  configuration  is  shown  in  Fig.  7-5b.   From 
Fig.  7-5b  it  can  be  visualized  that  the  response  is  speeded  up  by  the 
two  zeros  while  the  overshoot  is  small  because  of  the  approximately  equal 
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Fig.  7-5 (b)   Final  Closed  Loop  Pole-zero  Configuration  of 

Example  7-3 
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distances  of  the  zeros  and  the  dominant  roots  to  the  origin. 

7-4  Mixed  Compensators. 

For  a  single  section  cascade  or  feedback  compensator  it  has  been 
shown  that  the  adjustable  parameters  are  linearly  dependent  (Matrices  6-9 
and  7-5),   In  general  a  compensator  has  the  transfer  function 

r     _   apSP  +  CLr-/^-'  +  --•  -irQ-o  (7.20) 

The  coefficients,  Q    ,  (X,    ,  •-  Q„  and  b0  >    bt     --  b*     are  linearly 
dependent  for  either  a  cascade  or  a  feedback  path.   The  solutions  given 
a  set  of  roots  are  always  real  numbers.   If  there  are  more  than  one  compen- 
sator and  in  different  loops,  the  adjustable  parameters  are  not  linearly 
dependent.   The  solutions  for  a  given  set  of  roots  may  be  complex  numbers 
and  consequently  the  compensators  are  not  physically  realizable.   However, 
the  subset  of  roots  which  gives  the  solution  of  real  numbers  can  be  found 
if  such  a  subset  exists  in  the  stable  root  region.   Note  the  coefficients 
of  the  reduced  characteristic  equations  are  always  real  numbers  for  any 
choice  of  the  roots  and  for  any  type  of  compensators.   The  reason  for  this 
fact  is  due  to  the  complex  conjugate  roots.   Because  of  this  fact  the  stable 
root  region  and  dominant  root  region  are  always  defined  if  they  exist. 
The  physically  adjustable  parameters  also  appear  as  real  numbers  in  some 
transformations  for  any  choice  of  the  roots,  but  the  individual  parameters 
may  not  be  real  numbers.   Consider  Fig.  7-6  which  consists  of  a  single 
section  compensator  in  cascade  and  another  in  the  feedback  path,  Cq     is 
the  given  plant  and  2,  ,  £t  ,   p;  ,  f^  and  £,  are  adjustable  parameters. 
From  a  signal  flow  graph,  the  characteristic  equation  is 

z  (7-21) 
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Fig.  7-6    Mixed  Compensators 


Fig.7-6(a)  Mixed  Compensation  of  Example  7-4 
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Consider  the    first   term  of    (7-21),    it   can  be   expressed  as    follows: 

(i+P,)(S  +  P*)Q  =    [V*  Cfi+ft)*   *  P.fjtf  (7-22) 

Let 

P.-f  P2    =    b  (7-23) 

P,  ft  "  O.  (7-24) 

where  o  and  CL  are  always  real  numbers  for  any  choice  of  roots  of  (7-21). 
However,  p  and  p  are  not  necessarily  real.  To  find  the  real  solution 
of  p  and  p  ,  the  roots  must  be  subject  to  some  restrictions.  In  this 
case  it  is 

b     >^0<_  (7-25) 

The  inequality  (7-25)  defines  the  realizable  region.   For  a  passive  net- 
work, additional  constraints  are  necessary,  that  is 

b  >o  a  >o  (7-26) 

The  same  constraints  of  the  roots  can  be  applied  to  other  cases. 

Example  7-4.   Given        G-  —  ,-^~- ,   The  velocity  error 

H  iCS+s-jCit/oj2- 

constant  K  is  not  to  be  reduced.   The  band  width  is  less  than  ^0  rad/sec< 
v 

From  the  specifications  dominant  roots  of  f    =   0.5  Cj  =  10  are  chosen. 
For  cascade  compensation  alone,  more  than  two  lead  sections  are  required 
by  the  Ross  Warren  method  .   A  scheme  of  Fig.  7-6 (a)  is  attempted.   Be- 
cause of  the  gain  requirement,  only  second  derivative  feedback  is  used 
and  the  parameters  o(    ,  z  and  p  are  related  by  the  following  relation 

d  —  -£-  (7-27) 

There  are  three  free  parameters  (  <s(  ,  p  ,  j£  ),  three  roots  are  art  it - 
rary.   The  constant  controlled  characteristic  equation  is 

Sr+  (2L-  +  PJS*  +  {zoo+soop  +  id" 'tU.)Sl V  (too  +  *>of>  fto'X^s2 

4.    (i-OOp  +I0yocj  5    -p     ,</p    =0  (7-28) 
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Let  (-S-),  (-S-)  and  (-S  )  be  the  arbitrary  roots  and  define  the 
following  terms  for  simplification. 

F,  =.    S.+S^S,  (7-29-1) 

P2  —  _S,St  4-  S25,  +  5^,  (7-29-2) 

F3   =    S.S^Sj  (7-29-3) 
The   root-coefficient  relations  are: 

ZOO    tp    -      F,  4/  C,  (7-30-1) 

a00  4   Jr00p+I0*£c=    C*+    CF,     +   F4  (7-30-2) 

25    +  aoop  Wgj0    =   CoF,    +  C^  +  F3  (7-30-3) 

ZiTp    +  ^10*  =   CoF^   +CB  (7-30-4) 
\Q*?    =   C0F3 


(7-30-5) 


Where  C,  and  C  are  the  coefficients  of  the  reduced  characteristic  equation 
1      o 

which  defines  the  constrained  roots  (-S,)  and  (-S_)  by  the  equation 

5*+  C,S  +  C0  — O  (7-31) 

From  equations  (7-30-1),  (7-30-2),  (7-30-3)  and  (7-30-5),  C   is 

o 

obtained: 

F3(flPJ  -isFj+iQ^Co    -  iOL~[^^  -tSrP.Fj  f/O^JCo 

-f  10'°  [  F,  Fa  -  25"  F,  -  F3  +  50oJ  =  O  <7"32> 

Assume  (-S.)  and  (-S~)  are  the  dominant  roots  and  denoted  by  P     and  UJn   m 
(-S-)  is  arbitrary  and  must  be  greater  than  f^n    •   Assume  S-  =  ZftOn 
For  f  =  0.5,  U)n=   10  from  (7-29) 

F,  =  20 

Fa  =  200 

P  =  1000 
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Substitute  into  (7-32)  one  obtains, 


Take 

C  =  2080 
o 

From 

(7-30-5) 

From 

(7-30-1) 

From 

(7-30-2) 

From 

(7-30-4) 

C  =  +  2080, 

o         ' 


p  =  20.8 


Cx  =  25.8 


31.6 


/or4.=  2°76 

<X  =  4.32 
The  reduced  characteristic  equation  is 

52  +  zs.&s  +  zoao  =  o 

The  constrained  roots  are  -S, ,   -S_  =  -12.9  +  j  43.5 

4     5         — 


The  five  roots  are: 


The  compensator  is 


-S1,    -S2     =  -£    +   j   8.66 

-s3  =  -10 

-S4,    -S        =  -12.9     +     j   43.5 


-     gj3g  (S-K83J 
Ctc—       s  +  20.8 


i^=    0  0ZO1Q 


The  closed   loop  pole-zero  configuration   is   shovm   in  Fig.    7-7.      For  other 
choices  of  the  arbitrary  roots,    better  results   can  be  obtained. 

X 
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-  IO 


-30 
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Fig.  7-7  Closed  Loop  Pole-zero  Configuration 
of  Example  7-4. 
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VIII.   CONCLUSIONS 
8-1  Summary  of  Results. 

A  new  method  for  design  of  control  systems  in  the  time  domain  has 
been  described.   The  method  is  applicable  to  any  system  configuration  and 
is  able  to  handle  any  number  of  variable  parameters.   The  design  procedure 
is  based  upon  a  set  of  simultaneous  algebraic  equations,  and  is  well  suited 
to  implementation  by  digital  computer.   The  method  is  flexible  and  has  more 
freedom  in  design  than  other  methods.   The  values  of  adjustable  parameters 
are  computed  from  analytical  expressions,  so  any  desired  degree  of  accuracy 
can  be  obtained. 

The  main  advantage  of  this  method  over  those  methods  available  in  the 
literature  lies  in  the  fact  that  the  method  is  a  direct  syntheses  from  the 
closed  loop  pole-zero  configuration.   No  trial  and  error  is  needed.   The 
specifications  and  restrictions  of  the  compensator  confine  the  arbitrary 
roots  within  a  region  on  the  S-plane,  once  this  region  or  regions  have  been 
determined,  the  rest  of  the  design  work  is  just  a  computation.   The  compen- 
sator obtained  by  this  computation  guarantee  that  the  specifications  and  re- 
strictions are  satisfied. 

The  main  difference  between  this  method  and  other  available  methods 
is  that  the  domain  of  the  system  can  be  any  number  of  variables.   Theoreti- 
cally, the  more  variables  available,  the  more  design  freedom  one  has.   In 
the  root-locus  method,  only  one  parameter  is  varied  others  must  be  fixed. 
In  Mitrovic's  method,  only  two  parameters  are  varied.   This  method  is  able 
to  handle  any  number  of  variable  parameters  simultaneously. 
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8-2  General  Description  of  the  Method. 

The  design  procedure  can  be  summarized  as  shown  in  Fig.  8-1.   The 
analytical  expressions  of  the  dependent  variables  of  commonly  used  com- 
pensators have  been  derived  in  the  previous  chapters.   In  general  for  the 
case  in  which  coefficients  of  the  characteristic  equations  as  variables,, 
the  dependent  variables  are  always  in  linear  relations  one  to  another. 
A  compensator  of  the  transfer  function  as  (8-1)  is  in  the  forward  path  or 
in  the  feedback  path  alone,  the  coefficients  a's  and  b's  are  always  linear 

5m  +  b„Sn-  +  --  -  +    b0 

relations  in  the  partitioned  formulas.   Those  linear  algebraic  equations 
are  always  consistant  and  therefore  always  have  a  solution.   For  the  case 
in  which  compensators  are  in  both  forward  and  feedback  paths,  a  quadratic 
expression  is  always  obtained.   But  this  quadratic  expression  does  not 
imply  the  dominant  and  realizable  regions  are  smaller.   It  shows  only  the 
natural  algebraic  relations.   For  example,  the  coefficients  of  the  reduc- 
ed characteristic  are  always  real  no  matter  how  complicated  their  functions 
are. 

8-  3  Extensions. 

The  method  can  be  extended  to  sampled  data  systems  and  adaptive  sys- 
tems.  The  modifications  are  the  different  interpretations  of  the  variables, 
For  example,  in  the  case  of  sampled  data  systems,  if  £  transform  is  used 
then  the  criteria  of  stability,  dominance  and  realizability  are  different. 
However,  the  general  idea  of  approach  can  be  applied. 

In  the  previous  chapters,  the  arbitrary  variables  are  chosen  as  the 
roots  of  the  system.   This  choice  is  not  necessary.   Theoretically,  any 
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Fig.   8-1  General  Description,  of  the  Method 
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combination  of  "r"  variables  can  be  cbosen  as  the  arbitrary  variables  such 
as  a  combination  of  roots  and  variable  parameters. 

This  method  can  be  applied  to  any  mode  of  linear  operation  in  dis- 
continuous operation.   By  proper  choice  of  the  compensators  the  dominant 
roots  of  two  modes  can  be  chosen  in  a  relation  as  the  designer  wants. 
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APPENDIX  I 
CONSISTENCY  &  LINEAR  DEPENDENCE  OF  DERIVATIVE  FEEDBACK 
1.    Introduction. 

The  coefficients  of  the  characteristic  equation  of  a  linear  control 
system  are  functions  of  the  system  parameters.   If  all  the  coefficients 
can  be  adjusted  by  a  set  of  the  parameters,  then  the  roots  of  the  system 
can  be  controlled.   Conversely,  if  the  roots  are  chosen,  then  the  system 
parameters  are  forced  to  have  certain  values.   The  conditions  of  using 
system  parameters  to  satisfy  the  coefficients  of  a  characteristic  equation 
for  a  given  set  of  roots  depend  upon  the  number  of  adjustable  parameters 
and  the  consistency  of  the  simultaneous  equations.   Theoretically,  if  all 
the  derivatives  of  a  system  can  be  obtained  from  the  output  of  the  system, 
then  by  feedback  as  shown  in  Fig.  A-l-1,  the  system  roots  can  be  chosen 
arbitrarily  providing  no  restrictions  are  placed  on  the  forward  gain  K 
and  the  control  parameters  h's.   This  can  be  shown  as  follows.   The  charact- 
eristic equation  of  the  system  with  feedback  is: 

3n+UA_l+a()_jSn~'-f   4(  Kzh,i«,)S  +    K,KZ=D    (A-l-1) 

Assume  the  roots  of  the  system  are  chosen  from  the  considerations  of 
the  specifications,  and  the  required  characteristic  equation  is  formulated 
as  equation  (A-l-2) 


3%  /l^.S0"'^  +  fl2S  f  A,S  4  f\c   =0    (A-l-2) 


For  the  system  to  have  the  chosen  roots,  then  the  corresponding  coeffi- 
cients of  equations  (A-l-1)  and  (A-l-2)  must  be  equal.   Then 
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Fig.  A-l-1   Block  Diagram  by  Feedback  of  Derivatives  .  £r  is 
Transfer  Function  of  the  Plant >  K  and  K. 
are  the  System  Gains. 
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k,Ka  =  H0 

M,  +  &,  —  ^i  (A-l-3) 


Equation  (A- 1-3)  always  has  a  unique  solution  if  there  are  no  re- 
strictions on  K' s  and  h's.   And  in  the  particular  cases  each  parameter 
h  controls  one  coefficient  independently  and  the  forward  gain  controls 
the  constant  co-efficient.   If  the  control  parameters  have  certain  restric- 
tions, then  the  system  roots  cannot  be  arbitrarily  chosen^  but  have  to  be 
subject  to  these  restrictions.   For  example,  if  K.K^,  is  specified  from 
the  consideration  of  steady  state  accuracy,  then  the  product  of  the  roots 

must  equal  to  K, K_.   If  h  ,  cannot  be  obtained  (h  ,  =  0),  then  the  si 
12       n-1  n-1 


mm 


of  roots  must  be  equal  to  a   ,.   The  restrictions  on  the  control  para- 

n-1 

meters  restrict  the  roots  in  a  certain  constraint.   The  more  restrictions9 
the  higher  is  the  degree  of  the  constraints  of  the  roots.   In  a  practical 
system,  the  higher  derivatives  are  very  noisy  and  very  difficult  to  obtain, 
This  is  a  common  restriction  for  derivative  feedback. 

In  Fig.  A-l-1,  it  was  assumed  that  there  were  no  zeros  in  the  forward 
path.   If  one  zero  is  in  the  forward  path,  then  the  allowed  highest  order 
of  the  feedback  derivative  is  (n-2).   In  general,  if  Vn"  zeros  are  in  the 
forward  path,  only  the  derivatives  of  the  order  up  to  n-(M+l)  are  allowed 
to  feedback.   This  can  be  proved  as  follows 

Refer  to  Fig.  A-l-1,  assume  m  zeros  in  Gf      ,  and  the  highest  deriva- 
tive has  an  order  of  f       .   Then  the  characteristic  equation  is 


5n  +  (an.,+hr+m)  sn~'  +  -  --   =o  (A-i-4) 
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The  subscript  (  P+m  )  of  H  r+in  indicates  the  highest  order  of  the  feed- 
back derivative  and  the  number  of  zeros.   From  equation  (A-l-4) 

r  +  m    <.  n  -  i 
r  <    n-  (  m  +  ij 

Therefore  if  in    zeros  are  in  the  forward  path,  then  only  n-(  ,n+  I  ) 
control  parameters  can  be  introduced,  if  the  feedback  path  is  as  shown  in 
Fig.  A- 1-1.   However,  other  schemes  of  lower  order  derivative  feedback  are 
possible  to  introduce  (n-1)  control  parameters  and  consequently  control  all 
the  coefficients  of  the  characteristic  equation.   The  following  sections 
are  discussions  about  these  possibilities.   In  the  following  sections,  two 
terminologies  are  defined  as  follows: 

Control  variable  is  the  variable  of  the  system  which  can  be  measured. 
Control  parameter  is  the  adjustable  parameter  in  the  system. 
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2.    Feedback  of  the  First  and  Second  Derivative  f rom  t he  Output  to 
Several  Points  in  the  Forward  Path. 

Since  the  first  and  the  second  derivative  are  usually  allowed  to 
feedback  for  high  order  systems  and  can  be  obtained,  only  these  two  deriva- 
tives are  considered.   The  object  of  this  discussion  is  to  adjust  all  the 
coefficients  of  a  nth  order  characteristic  equation  to  desired  values.   In 
order  to  accomplish  this,  n  control  parameters  are  needed.   As  the  forward 
gain  K  is  considered  as  one  control  parameter,  (n-1)  control  parameters 
from  the  feedback  path  are  necessary  for  this  particular  scheme.   Using 
first  and  second  derivative  feedback  around  the  plant  intorduces  only  two 
control  parameters.   However,  if  there  exist  junction  points  between  the 
energy-storage  elements  in  the  forward  path,  such  that  signals  can  be  fed 
in  or  picked  up,  then  additional  control  parameters  can  be  introduced  by 
feeding  derivatives  to  those  points.   As  two  control  parameters  have  been 
introduced  by  feeding  derivatives  to  the  error  channel,  then  the  neces- 
sary number  of  points  between  the  energy  storage  elements  in  the  for- 

_  _  o  n  —  ? 

ward  path  is  —z —   (if  n  is  cdd  and  ^ —   if  n  is  even)   since  two 
derivatives  can  be  fed  into  each  node.   Briefly,  this  can  be  summarized  in 

the  following  two  statements. 

n-3 
(1)*  If  a  nth  order  system  has  at  least  — r—  available  points 

(if  n  is  odd,  and  — r—  if  n  is  even)  between  the  energy  storage  elements 
in  the  forward  path,  then  all  the  coefficients  of  the  characteristic  equa- 
tion can  be  adjusted  to  arbitrary  values  by  introducing  first  and  second 
derivatives  from  the  output  to  each  point  providing  the  two  derivatives 
are  defined  and  the  feedback  paths  are  proper. 

*all  statements  in  this  discussion  may  be  considered  as  rules.   All  of  them 
are  necessary  conditions,  but  not  sufficient  conditions. 
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(2)   If  a  nth  order  system  has  at  least  n-2  available  points 
between  the  energy-storage  elements  in  the  forward  path,  then  all  the  co- 
efficients of  the  characteristic  equation  can  be  adjusted  to  arbitrary 
values  by  introducing  only  first  derivative  feedback  from  the  output  to 
each  point,  providing  the  feedback  paths  are  proper. 
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3.    Formulation  of  Characteristic  Equation  for  Pure  Derivative  Feed- 
back. 

From  the  signal  flow  graph,  the  determinant  of  a  system  can  be 
obtained  by  the  following  rule. 

where  z=.  product  of  the  kth  possible  combination  of  A.   non- 

t—  K 

touching  loops 

/  ,  =  loop  transmission 
For  this  particular  derivative  feedback  systems  two  terminologies 
are  defined  as  follows: 

Loop  gain  =  numerator  of  each  loop  transmission. 
In  Fig.  A-l-2.   Loop  gain  of  L.  =  h.P 

Loop  gain  of  L2  =  h2-2P3 
Loop  gain  of  L-  =  h  P  P2P 
Loop  gain  of  L,  =  PjP2P3 
Characteristic  function  =  denominator  of  each  of  the  forward  transfer 
functions,  or  the  product  of  them. 

Consider  Fig.  A-l-2:   from  the  graph  rule 

+  P.fiPjhj  +  ^?S,  (A-l-5) 

The  first  term    Q   Q  Q,   4    last  term  P  P_P  =  characteristic  equation 
of  the  original  system. 
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Fig.  A- 1-2  Formulation  of  the  Determinant 

P;(s),  Q,-(s),  b;(sj  ore  polynomial  ,*p  5 
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The  second  term  =  loop  gain  L.  X  nontouching  chara.  function. 
The  third  term  =  loop  gain  L_  X  nontouching  chara.  function. 
The  fourth  term  =  loop  gain  L_  X  nontouching  chara.  function. 
From  equation  (A-l-5),  the  formulation  of  the  characteristic  equa- 
tion for  pure  derivative  feedback  may  be  stated  in  the  following  rule. 
The  characteristic  equation  of  a  pure  derivative  feedback  system  = 

characteristic  equation  of  the  original  system  +     ^~      loop  gain  X 

K 

(*) 
characteristic  function  of  nontouching  loop  +      *^~  loop  gain  X 

characteristic  function  on  nontouching  loop. 

Example  A-l-1:   Devise  a  derivative  feedback  scheme  of  Fig.  A-l-3, 
such  that  all  the  coefficients  of  the  characteristic  equation  can  be  ad- 
justed.  Assume  the  required  characteristic  equation  from  the  specifica- 
tions is 

5*"  +  A4s*+  aj5*4  &^+  A«s  +  ^o  —  o 

find  the  relations  of  the  control  parameters  and  the  coefficients  of  the 
given  characteristic  equation. 

Since  intermediate  points  are  available,  then  n-2  =  5-2  =  3.   From 
rule  (2),  it  is  not  possible  to  use  first  derivative  only. 

By  rule  (1)  — t—     =     -z —  =  1.   Then  a  combination  of  first  and  second 
derivative  can  meet  the  requirement.  There  are  many  possible  ways  to 
accomplish  this  since  only  4  control  parameters  are  needed. 

Two  schemes  are  listed  as  follows: 

(a)  Hx  =  t^s       H2  =  h2s2  +  h3s     H3  =  t^s2 

(b)  Hj   =  h^  H2   =  h2s2  H3  =  1y2   +  h4s 

Assume  scheme  (a)  is  chosen.   The  chara.  equation  is  formulated  in 
a  tabulated  form  for  convenience. 
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(5+P2)(s+8) 


U 


L  H,  =  1S; 


Fig.  A-l-3.   Example  A-l-1  -  Mark  "X"  denotes  signals  are  allowed 

to  feed  in  or  pick  up 
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Required  Chara.  Equation 
Original  Chara.  Equation 
Gain  [_.   x  Chara.  Nontouching 

Gain   [_j   x  Chara.  Nontouching 


Gain  L^,   x  Chara.  Func.  Nontov 

h,s  Kj  *(s+ftj(*+fcK**ft) 


A*S 


vl 


4*5        A,5      flo 
Z£Pj    |ZPC     PK      77  (^    *  W*J 


k,h, 


K.Ww 


K2K3h; 


*A(p.*p2 


K,w  H+ 


Wh, 


MP.6  + 


iy   gp^p.k 


iC^ 


The  equations  for  the  coefficients  are 

2  Pift  +  K.fcfc  h«.  f    K^  h^  +   ^O.  +  Pz*  ftj^i,  =  ^3 
TT  P^  f  fclCKiKj  f   kRj  p,£  +  KsP.^&^i   -  /»i 

The  first  four  equations  of  (A-l-6)  can  be  written  in  matrix  form: 


K3P.P.P3 


0 
KZJC3 

K*  «3  p, 


0 
o 


KZK3 

Wifi 

\5L 


0 

f 

k",K2Kj 

h2 

_-. 

K.K.K^ 

Hi 

0 

— . 

>_ 

V 

-x  P. 

A3 

-*&&■ 

A, 

-*w« 

fl,- 

77(^-  n. 

In  (A- 1-7)  whenever  the  determinant  is  not  zero,  there  always  exists 
a  solution  of  the  control  parameters  h's  for  arbitrary  values  of  A. 

In  the  above  example,  the  table  of  the  characteristic  equations  gives 
a  clear  picture  of  how  the  coefficients  are  controlled  by  h-parameters. 
Also  equation  (A-l-7)  shows  that  the  h-parameters  are  linearly  dependent 
to  the  coefficients  of  the  given  characteristic  equation.   It  is  this 
linear  dependence  which  makes  the  arbitrary  adjustment  possible.   In 
equation  (A- 1-6),  the  constant  term  is  controlled  by  the  forward  gain  only. 
This  is  the  nature  of  derivative  feedback,  because  the  derivatives  control 
the  coefficient  of  the  terms  of  a  differential  equation.   Any  derivative 
feedback  from  the  output  has  the  scheme  of  properties  of  linear  dependence 
of  h-parameters  and  the  independence  of  forward  gain  K. 

Example  (A-l-2).   Fig.  A-l-4  is  a  third  order  system,  one  point  between 
the  energy-storage  elements  in  the  forward  path  is  allowed  to  feed  in  signal, 
then  the  system  roots  may  be  arbitrarily  chosen.   Assume  the  roots  to  be 
-3  and  -0.5  +  J  1,   K. ,  h  ,  h_  can  be  determined.   The  characteristic  equa- 
tion for  the  given  roots  is; 

S3   +    4S2  +    4.25S   +  3.75  =  0 

Main  loop  1         1        0         K 

L.  :  loop  gain  =  Lh.s 
chara.  function 
not  enclosed  =  1  ^in9 

L~ :  loop  gain  =  h  s 
Chara.  function 
not  enclosed  =  s  +  1  h        h 
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Fig.    A-I-5     Example   a-I-3 
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Set  the  corresponding  coefficients  equal,  variable  parameters  are  solved. 
KT  =  3.75,         h  =3,  h2  =  0.334 

Example  A-l-3.   In  Fig.  A-l-5,  it  is  assumed  that  signal  can  be  fed  in  at 
the  point  "X"  in  the  forward  path,  and  the  forward  gain  of  the  system  is 
1500  based  on  the  steady  state  specification.   According  to  rule  1,  — r—  =  1, 
then  by  first  and  second  derivative  feedback,  the  system  roots  can  be  chosen 
arbitrarily  providing  the  product  of  the  roots  is  equal  to  1500. 

The  uncompensated  system  is  unstable.   The  feedback  paths  are  chosen  as 
shown,  and  the  roots  are  assumed  to  be  S  =  (-3  +  j5)9  -6,  and  -7.36. 
The  characteristic  equation  for  the  chosen  roots  is 

S4  +  19.36S3  +  158. 52  +  720S  +  1500  -  0 
The  system  coefficients 

Uncompensated  system  =1  10        31        30        1500 

L,:  (  h,S+h^)      1500  1500h2   1500h1 

L2:  h3S   30(5  +  *Xi  +  $-J  30h3      240h3     450h3 

The  coefficients  relations  are: 
30h3  +  10  =  19.36 
1500h2  +  240h3  +  31  =  158 
450h3  +  1500h1  +  30  =  720 
Solve  equation  (A-l-8),  get  (A-l-8) 


h1  =  0.367 
h2  =  0.0347 
h3  =  0.312 
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5.    Feedback  from  Intermediate  Points. 

It  can  be  shown  that  if  the  output  of  each  block  in  the  diagram  (A-l-6) 
is  chosen  as  the  state  variable  of  the  system,  then  there  exists  a  linear 
transforamtion  between  the  state  space  and  the  phase  space  (output  and  all 
its  derivatives)  of  a  system.   In  Fig.  A-l-6,  the  state  variables  are 
chosen  as  shown,  the  relationship  between  state  space  and  phase  space  is 


<r~ 


3tCL 


Kb 


5+b 


3 


X,  =  6>«, 


derived  as  follows: 


Fig.  A-l-6  State  Variables 


x,=  ec 


(A-l-9-1) 


PC,  __  _L 

Xi    5  j 


(A-l-9-2) 


-rr-M^  ^ +  *>*<  -  ^ ,  «,  =  -i&j  ♦  t&) 


(A-l-9-3) 


Equations  (A- 1-9)  are  the  required  relations.   It  shows  that  the  out 
put  of  each  block  is  a  linear  combination  of  the  output  0      and  its  de- 
rivatives. 
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With  this  concept,  therefore,  the  feedback  from  intermediate  points 
is  equivalent  to  derivative  feedback  from  the  output.   Fig,  A-i-7(a)  is 
equivalent  to  Fig.  A-l-7(b). 

This  equivalence  also  can  be  shown  by  block  diagram  manipulations. 
Because  of  this  equivalence,  if  the  signals  flowing  in  the  intermediate 
points  can  be  measured,  then  those  signals  can  be  taken  as  control  vari- 
ables and  the  analysis  is  the  same  as  for  derivative  feedback  from  the 
output.   Moreover,  if  the  derivatives  of  those  signals  are  defined  and 
can  be  obtained,  they  also  can  be  taken  as  control  variables. 

However,  if  consistency  and  linear  dependence  are  considered,  the 
feedback  paths  and  the  number  of  control  parameters  have  many  restrictions. 
(1)  Consider  the  scheme  of  Fig.  A-l-8,  it  introduces  sufficient 
control  parameters  but  the  coefficients  of  the  characteristic  equation  are 
not  linear  combinations  of  the  h-parameters.   It  consists  of  two  indepen- 
dent nontouching  loops,  the  contribution  of  these  two  loops  to  the  whole 
system  is  the  product  of  the  two  loop  transmissions.   So  the  coefficients 
are  quadratic  functions  of  the  h-  parameters.   This  can  be  shown  by  formu- 
lating the  characteristic  equation  as  follows: 


s3 

s2 

S1 

s° 

Original   system 

1 

ci-fb 

c\b 

KaK 

lx  :    kbh,5  CS+*J 

Kbh, 

KP^, 

4      K^MsVbsJ 

KA 

fCi)ha 

L.    x  L2 

KK^x 

If  the  required  characteristic  equation  is 

3         2 

S    +   A0S  +  A,S   +  A   =  0 
2        1       o 
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Fig.  A-l-7(a)   Feedback  from  Intermediate  points, 
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Fig.   A-l-7(b)     Equivalenence  of  Fig.   A-l-7(a) 
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Fig.  A-l-8  Quardrtic  Function  of  h  parameters 
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Fig.  A-l-9  Example  A-l-4 
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then 

fCaQH,  4  J^bh,  +  1^^  =  4,-afc  (A-i-io) 

-  K  K  =  A  0 

The  solution  of  equation  (A- 1-10)  for  real  values  of  h„  and  H~  implies 
the  restriction  of  the  roots  in  a  certain  region.   In  generals  if  there 
are  "  /\    "  independent  feedback  loops  in  a  system^  the  coefficients  of 
the  characteristic  equation  consist  of  the  terms  h . ,  (h  h. ) , . . .  (h. . . .h  ). 
This  leads  to  the  following  statement. 

If  all  the  coefficients  of  the  characteristic  equation  of  a  derivative 
controlled  system  are  a  linear  function  of  the  control  parameter,  then  all 
feedback  paths  must  be  interconnected,  (or  touching).   The  above  state- 
ment may  be  stated  in  another  way:   If  all  the  coefficients  of  the  character- 
istic equation  of  a  derivative  controlled  system  be  a  linear  function  of 
the  control  parameters,  then  those  points  from  which  the  signals  are 
taken  as  control  variables  cannot  be  used  as  feed-in  nodes. 

Example  A-l-4.   In  Fig.  A-l-9,  assume  the  signal  flowing  at  the  point 
marked  "X"  can  be  measured,  then  it  can  be  picked  up  as  a  control  variable. 
Assume  the  feedback  paths  are  as  shown,  they  are  interconnected.   The 

characteristic  equations  are: 

3        2 
Required  characteristic  equation  =  S   +  A„S   +  A, S   +  A  =0 

2        1      o 

Characteristic  equation  J      a-*-k>      ab     KIC 

Li  ■  Ks  K°«>  K,K,h, 

The  coefficient  relations  are: 
K0hL=  A,-  (a+bj 

Kx  Kb  ~  ^  ° 
h.  and  h.  can  be  solved  for  arbitrary  value  of  A's. 
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(A-l-11) 


(2)   If  zeros  are  Involved  in  the  forward  path,,  the  feedback 
from  the  intermediate  nodes  may  introduce  the  problem  of  consistency. 
Consider  Fig.  A-l-10,  the  signal  at  "X"  can  be  measured  and  then  can  be 

taken  as  control  variable.   The  characteristic  equations  are: 

3         2 
Required  Chara.  Equation      S    +   A~S     +   A, S   +  A  =  0 

2  1        o 

Original  loop  1        a  +  b     ab+k  k,     K  K 

The  coefficient  relations  are: 

hcA  =1 


(A-l-12) 


Equation  (A-l-12)  are  linear,  but  are  not  consistent  for  arbitrary 
A's,  since  the  number  of  equations  is  one  more  than  the  number  of  unknowns 
(h  and  h«).   Substitute  h  from  the  first  equation  into  the  second  and 
third  equation,  obtain: 

KafCtA  +  C^  +  bJ  =  fl2-(a+  bj 

Manipulate,  get 

cx/\2-/q,  —  <x  (<*  +  a  - 3b)  -  ab  -  I^Kt,  (A-1-13) 

Equation  (A-1-13)  is  the  restriction  of  the  roots  for  this  systm,  i.e.s  if 
equation  (A-l-12)  has  a  solution,  then  equation  (A-1-13)  must  be  satisfied. 
Therefore,  when  a  system  has  the  required  number  of  control  parameters  and 
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Fig.  A-l-10.   Effect  of  Zeros  to  Feedback 


Fig.  A- 1-11.   Feedback  of  the  output  to  more  than  one  point. 
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the  linear  relationship  of  the  coefficients9  the  choice  of  roots  may  still 
be  subject  to  some  restriction  if  the  feedback  is  not  proper.   In  general  it 
can  be  shown  that  if  the  loop  transmission  has  the  same  order(s)  for  both 
numerator  and  denominator,  then  it  introduces  the  restriction  of  the  roots. 
This  leads  to  the  following  statement. 

If  the  roots  of  derivative  controlled  systems  are  to  be  chosen  arbi- 
trarily, then  for  each  loop  transmission,  the  order  of  the  numerator  must 
be  at  least  one  less  than  the  denominator. 

(3)  There  is  another  obvious  inconsistency  which  arises  from 
feedback  of  the  output  to  more  than  one  point.   Therefore,  if  an  inter- 
mediate signal  which  is  not  followed  by  pure  integrator,  be  fed  to  an- 
other point,  the  equations  of  the  coefficients  is  inconsistent  for  arbi- 
trary choice  of  roots.   Fig.  A-l-11  is  an  example. 


•  .  i    *  „3        „2        1        o 

original  loop  __^       S        S        S        S 


The  characteristic  equation  of  Fig.  A-l-11  is 

3        2        1 
S        S        S 

/  Q+b  +  c        ab+bC  +  ba        K^il+nbC 

^^     Kh,  (<\+b)K,h,  ObK,h, 

There  are  three  equations  but  two  unknowns.   This  leads  to  a  con- 
straint of  the  roots. 

(4)   Since  the  feedback  from  an  intermediate  node  is  equivalent 
to  derivative  feedback  from  the  output,  then  it  is  possible  that  the  para- 
meters introduced  by  both  control  variables  may  have  the  same  effect  on 
the  system;  i.e.,  the  parameters  may  not  be  independent.   This  can  be 
easily  detected  by  just  looking  at  the  loops.   In -Fig.  A-l-12^,  h_ 
has  the  same  effect  as  h„  and  h  ,  since  h  is  a  combination  of  the  first 
and  second  derivative  of  the  output.   That  iss  h  enters  the  coefficient 
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of  the  first  and  second  derivative  terms  as  h  and  h_.   Therefore  in 
Fig.  A-l-12,  the  system  actually  has  only  two  control  parameters. 

However,  if  the  derivative  of  the  signal  at  "X"  can  be  obtained 
and  the  scheme  changes  to  Fig.  A-l-13s  then  each  parameter  plays  its  own 
role. 
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Fig.  A-l-12.  Two  parameters  have  the  same  effect. 
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5CS+cj 
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Fig.  A-l-13.   Each  parameter  has  it's  own  effect 
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5.    Summary. 

Derivative  control  is  an  effective  method  to  control  the  coefficients 
of  the  characteristic  equation  of  a  system,  since  it  uses  derivatives  of 
the  system  variables  to  control  the  coefficients  of  the  characteristic 
equation.   Since  an  intermediate  signal  feedback  has  the  same  effect  as 
derivatives  from  the  output,  they  are  considered  as  in  the  category  of 
derivatives.   The  criteria  in  Section  I  also  applies  to  the  intermediate 
signal  feedback  because  of  the  condition  (1)  in  Section  4.  The  inter- 
mediate signals  do  not  introduce  more  control  parameters  but  increase  the 
possible  ways  of  obtaining  derivatives  and  the  feedback  path. 

In  general,  one  constraint  of  the  roots  is  almost  always  required  in 
a  system, that  is,  the  product  of  the  roots  equal  to  the  system  forward 
gain,  because  the  gain  is  usually  specified  from  the  stiffness  of  the 
system.   This  is  only  one  degree  of  constraint,  it  does  not  introduce 
much  difficulty.   But  if  additional  constraints  are  involved  such  as 
quadratic  functions  and  inconsistent  equations,  then  the  solution  is  not 
so  obvious.  Therefore,  it  is  logical  to  look  for  a  set  of  linear  and 
consistent  control  parameters.   From  the  discussion  in  the  previous  sec- 
tions, it  may  be  summarized  as  follows; 

(1)  The  number  of  h  parameters  must  be  equal  to  (w-1) 

(2)  h  -  parameters  must  appear  only  in  the  coefficients  from  the 
first  derivatives  to  the  (n-1)  derivative. 

(3)  h-parameters  must  be  linear  combinations  of  the  given  co- 
efficients. 

Example  A-l-5.   Fig.  2-1-14  is  a  6th  order  system.   Assume  two  nodes 
between  the  energy-storage  element  are  allowed  to  feed  in  or  pickup  signals, 
and  the  derivative  of  the  signal  at  "X".  can  be  obtained.   From  the 
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6-2 
criteria  of  intermediate  nodes,  (  — ~—  =  2)  the  system  coefficients  can 

be  adjusted  to  arbitrary  value  by  first  and  second  derivative  feedback. 

Fig.  A-l-14  (a)  (b)  and  (c)  are  the  possible  proper  feedback  paths.   Fig. 

A-1-I4(c)  again  has  5  possible  combinations  of  H°s.   They  are 

Hl=h,5+hz  Hx  -  h,s(  or  ht)  ^  =  h,5  +  h2 

H2   =    h3S+h1.  H2   =    hzS  +  h3  H2   =    h35 

H3   =    hfS  (  <n,M*J  H3   =     h4S*  4  h^  H3   =    h4*2  +  M 

Therefore,  there  are  8  possible  ways  to  control  the  coefficients  of  the 
characteristic  equation  of  this  system.   The  choice  of  the  scheme  depends 
upon  the  specifications  and  perhaps  the  engineering  judgments. 
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Fig.   A-l-  14.      Example  of  A-l-5. 
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APPENDIX  II 
PARTITION  OF  A  nTH  ORDER  SYSTEM  WITH  ONE  ZERO  FOR  SINGLE  SECTION  CASCADE 

COMPENSATOR 
Take  a  4th  order  system  as  shown  in  Fig.  A-2-1. 


(S  +  p) 


K(S+C) 


5(sV/+fis  +  PJ 


Fig.  A-2-1.   Cascade  Compensator 
The  controlled  characteristic  equation  is 

The  characteristic  equation  of  the  uncompensated  system  is 


(A-2-1) 


5^+  fts*4  fti'-f  ((3,  +  xji  4-  /CC  =  O 


A-2-2) 


In  order  to  express  the  coefficients  of  (A-2-1)  in  terms  of  the  coefficients 
of  the  uncompensated  system  in  a  simpler  form9  define: 
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ft  =  B* 


P,+  K 


B, 


(A-2-3) 


ICC  =  6, 

Substitute  (A-2-3)  into  (A-2-1)  and  manipulate,  one  obtains 

+  (  pB,  +  oCB0JS  +  B0p  =  0 


(A-2-4) 


Set  the  coefficients  of  (A-2-4)  equal  to  the  corresponding  root-coeffici- 
ents and  denote  the  arbitrary  roots  by   9      and  ud    s    one  obtains 


I       0  0  0-| 

6,    Bo   '^jH    -<*£     ° 
8a    0    -u>*  O     O 


p 

c„ 

■^T 

c, 

O 

c, 

O 

L 


(A-2-5) 


J 


Where  C  ,  C  and  C«  are  the  coefficients  of  the  reduced  characteristic  equa- 
tion. 

By  Cramer's   Rule. 

c  _  gjgrffj^  -<'- *0 t>A -h? +»K  ]  ~  'th-B2<^J<ii-nX]  (A_,  6) 


Define 
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: 


1^-1^%^^)  ^r    K>1 


For  an  uncompensated  nth  order  system  which  has  one  open  loop  zero  at 
(-C)  and  has  the  characteristic  equation  with  coefficients 


n-ls  n-2 


....B  ,  then 
o 


0     [a0m^r^n%%^--]^ 


(A-2-7) 
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